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Abstract 

This  document  presents  a  theoretical  study  di-muonic  hydrogen  and  helium 
molecules  that  have  the  potential  of  enhancing  the  muon-catalyzed  fusion  reaction  rate. 
In  order  to  study  these  di-muonic  molecules  a  method  of  non-adiabatic  quantum 
mechanics  referred  to  as  a  General  Particle  Orbital  (GPO)  method  was  developed.  Three 
mechanisms  that  have  the  possibility  of  enhancing  the  muon-catalyzed  fusion  rate  were 
discovered.  Two  involve  the  formation  of  di-muonic  hydrogen  molecules,  and  the  other 
uses  di-muonic  molecules  to  liberate  muons  stuck  to  He  nuclei.  The  effects  of  muon 
spin  on  di-muonic  hydrogen  molecules  was  studied.  The  nuclear  separation  in  di-muonic 
hydrogen  molecules  with  parallel  muon  spin  is  too  great  for  the  molecules  to  have  a 
fusion  rate  which  can  enhance  the  fusion  yield.  The  possibility  of  these  molecules 
transitioning  to  single  muon  molecules  or  triatomic  oblate  symmetric  top  molecules 
which  may  fuse  faster  is  examined.  Using  two  muons  to  catalyze  He-  He  fusion  is 

■j 

shown  to  be  impractical;  however,  using  two  muons  to  catalyze  He-d  fusion  is  possible. 
While  studying  the  physical  properties  of  di-muonic  hydrogen  and  helium  molecules 
some  unique  properties  were  discovered.  Correlation  interactions  in  these  molecules 
result  in  an  increase  in  the  calculated  nuclear  bond  length. 
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THEORETICAL  STUDY  OF  THE  EFFECTS  OF  DI-MUONIC 


MOLECULES  ON  MUON-CATALYZED  FUSION 

I.  Introduction 

Since  its  discovery  in  1947  researchers  have  hoped  that  muon-catalyzed  fusion 
could  provide  a  means  of  providing  an  almost  endless  supply  of  fusion  energy.  While 
great  strides  have  been  made  in  understanding  the  muon-catalyzed  fusion  process,  so  far 
the  yields  obtained  fall  short  of  what  has  been  hoped  for.  Some  proposed  mechanisms 
for  producing  energy  via  hybrid  fusion-fission  reactors  hold  some  promise,  however,  no 
practical  path  to  useful  quantities  of  pure  fusion  energy  has  yet  been  found. 

This  document  outlines  a  unique  method  of  using  non-adiabatic  quantum 
mechanics  to  study  muonic  molecules.  The  method  is  used  to  study  two-muon  catalytic 
processes  that  are  predicted  to  have  significant  effects  on  the  overall  fusion  yield  when 
the  thennal  muon  flux  is  large  (see  Chapter  8).  Reactions  which  have  the  potential  of 
increasing  the  muon-catalyzed  fusion  rate  and  reactions  that  could  free  muons  stuck  to 
helium  nuclei  are  presented. 

The  second  chapter  in  this  document  outlines  many  of  the  properties  of  muons 
which  are  important  to  muon-catalyzed  fusion  and  gives  a  history  of  muon-catalyzed 
fusion  research.  Included  in  this  chapter  are  sections  which  discuss  how  single  muon 
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catalyzed  fusion  occurs  and  many  of  the  factors  that  affect  fusion  yields.  The  chapter  is 
concluded  with  a  short  discussion  on  uses  of  muon-catalyzed  fusion  with  currently 
obtainable  yields. 

Chapter  3  gives  an  introduction  to  the  two  muon  processes  that  are  discussed  in 
detail  throughout  the  remainder  of  this  document.  Changes  to  the  muon-catalyzed  fusion 
cycle,  that  are  expected  to  occur  when  the  thermal  muon  flux  is  high,  are  presented. 

Chapter  4  presents  a  general  multi-configurational  quantum  mechanical  method 
of  modeling  molecules  that  allows  any  number  of  any  type  ( i.e .,  mass  and  charge)  of 
quantum  particles  in  the  presence  of  fixed  (i.e.,  classical)  particles  to  be  studied.  The 
methods  presented  in  Chapter  4  are  particularly  useful  for  studying  exotic  particles  that 
cannot  be  accurately  modeled  using  the  Bom-Oppenheimer  approximation.  Of  particular 
importance  to  this  chapter  are  the  methods  of  calculating  correlation  interactions  between 
particles.  These  interactions  are  shown  to  contribute  significantly  to  the  physical 
properties  of  the  muonic  molecules  studied. 

In  Chapter  5  the  effects  of  correlation  interactions  on  the  binding  energy  and  bond 
length  of  di-muonic  hydrogen  molecules  are  presented.  It  is  shown  that  the  relative 
impact  of  these  interactions  is  much  greater  in  muonic  molecules  than  it  is  with  similar 
conventional  molecules  which  contain  only  protons,  neutrons,  and  electrons.  These 
interactions  are  responsible  for  much  of  the  binding  energy  di-muonic  hydrogen 
molecules  contain  and  they  have  a  significant  impact  on  the  nuclear  bond  lengths.  It  is 
shown  that  these  correlation  interactions  result  in  the  nuclear  bond  lengths  being  larger 
than  what  would  occur  if  these  interactions  did  not  exist.  These  results  are  surprising,  in 
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that  they  are  significantly  different  than  what  has  been  observed  in  conventional 
hydrogen  molecules. 

Of  all  of  the  factors  which  limit  the  maximum  obtainable  muon-catalyzed  fusion 
yield,  none  is  more  important  than  the  fact  that  muons  sometimes  stick  to  helium  nuclei. 
Chapter  6  analyzes  the  possibility  of  using  di-muonic  fusion  processes  to  free  muons 
stuck  to  helium-3  nuclei.  Quasi-classical  molecular  dynamics  is  used  to  calculate  the 
vibrational  energy  levels  and  vibrational  spectra  of  the  di-muonic  molecules  presented  in 
this  section. 

Many  of  the  di-muonic  hydrogen  molecules  which  may  fonn  in  a  muon-catalyzed 
fusion  reaction  chamber,  with  a  high  thermal  muon  flux,  have  bond  lengths  that  are  too 
large  for  the  molecules  to  have  a  significant  fusion  rate  (see  Chapter  5).  If  these 
molecules  are  to  contribute  to  an  increase  in  the  muon-catalyzed  fusion  yield  a 
mechanism  for  them  to  transfer  energy  and  form  molecules  that  do  have  high  fusion  rates 
must  exist.  Chapter  7  analyzes  some  reaction  paths  that  could  transfonn  di-muonic 
hydrogen  molecules  with  relatively  large  bond  lengths  to  molecular  species  that  have 
much  more  closely  bound  nuclei.  Quantum  rotational  and  vibrational  energy  levels  are 
calculated  as  part  of  this  analysis. 

Chapter  8  of  this  document  examines  the  kinetics  of  reaction  paths  that  could  lead 
to  the  formation  of  di-muonic  hydrogen  molecules.  Upper  bounds  on  the  fusion  yield 
that  could  result  from  different  reaction  paths  are  determined.  Additionally,  lower 
bounds  of  the  muon  flux  that  could  result  in  these  reaction  mechanisms  contributing 
significantly  to  the  overall  fusion  yield  are  determined.  From  this  information  the  scope 
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of  reaction  paths  which  have  the  possibility  of  enhancing  muon-catalyzed  fusion  is 
narrowed  and  reaction  paths  which  can  be  neglected  are  identified. 
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II.  Background 


This  chapter  gives  an  overview  of  some  of  the  properties  of  muons  and  the  muon- 
catalyzed  fusion  process  as  it  is  currently  understood.  Some  of  the  unique  properties  of 
muons  and  the  manner  in  which  they  have  been  used  are  discussed.  Emphasis  is  given  to 
the  history  of  muon-catalyzed  fusion  and  the  factors  which  affect  obtainable  yields.  The 
chapter  ends  with  a  short  discussion  of  possible  uses  of  muon-catalyzed  fusion  with 
currently  obtainable  fusion  rates. 

The  theoretical  systems  studied  in  this  dissertation  are  limited  to  muonic 
molecules;  however,  the  models  developed  can  be  applied  directly  to  other  types  of 
exotic  particle  systems.  Additionally,  the  methods  presented  in  this  document  allow 
properties  of  conventional  molecules  that  can’t  be  calculated  using  most  ab-initio  code  to 
be  studied.  Examples  of  this  are  K  x-ray  isotope  shifts  and  differences  in  the  binding 
energy  of  different  isotopes  of  the  same  type  of  atom.  [1:51] 


2.1  Muons 


There  are  two  types  of  muons:  one  which  has  a  -1  charge,  the  same  as  an  electron, 
and  its  antiparticle  which  has  a  +1  charge,  the  same  as  a  proton.  The  rest  mass  of  a  muon 


is  206.7682823  times  that  of  an  electron. 


207/n  or  ~  —  m„ 

e  9  p 


.[2;  3]  As  a  result,  a 


negative  muon  is  often  referred  to  as  a  heavy  electron. 

Muons  are  unstable  particles  with  relatively  short  half-life’s  of  about 
2.20  x  lO^A1  .[2;  3]  The  actual  half-life’s  of  the  particles  varies,  depending  on  if  the 
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particle  is  negative  or  positive,  bound,  or  free,  etc.  Muons  can  decay  by  several  paths. 
By  far  the  most  common  is  for  positive  muons  ( ju+ )  to  decay  into  a  positron  (e  ),  electron 
neutrino  ( ve )  and  muon  anti-neutrino  (v  ),  and  for  negative  muons  (ju~)  to  decay  into  an 

electron  ( e  ),  an  electron  anti-neutrino  ( ve ),  and  a  muon  neutrino  (v  ).[4:2] 

ju+ - >e++v^+ve  (2.1) 

jLi  - > e~+vfi  +  ve  (2.2) 

There  are  two  primary  sources  of  muons,  decay  products  of  accelerator  produced 
particles  (>140  MeV  particles)  [4:17;  5:367-385]  and  decay  products  from  particles 
produced  when  primary  cosmic  rays  interact  with  the  atmosphere. [6:323]  Accelerator 
produced  muons  can  be  produced  and/or  moderated  to  yield  muons  with  energies 
between  about  an  eV  and  100  GeV.  Cosmic  ray  muons  have  typical  energies  in  the  GeV 
to  TeV  range.  The  moderation  of  cosmic  ray  muons  is  problematic.  Due  to  their  energy, 
cosmic  muons  are  extremely  penetrating,  and  when  moderated,  the  resulting  flux  is  so 
low  as  to  make  them  unusable  for  most  applications  requiring  muons  in  the  thennal  and 
near  thermal  range.  As  a  result,  research  which  requires  low  energy  muons,  such  as 
studies  involving  the  absorption  of  muons  by  matter,  are  almost  always  perfonned  using 
accelerator  produced  muons.  Research  which  can  use  higher  energy  muons,  such  as 
studies  involving  the  scattering  of  muons  are  often  performed  using  cosmic  ray 
muons.  [4] 
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The  ability  of  muons  (both  positive  and  negative)  to  probe  the  composition  of 
material  and  the  ability  of  negative  muons  to  catalyze  nuclear  fusion  is  of  particular 
interest.[7;  8;  9;  10;  11;  12]  In  order  to  use  muons  to  study  the  composition  of  matter, 
muons  scattered  by  matter  can  be  analyzed,  or  molecular  compounds  containing  muons 
can  be  studied  by  many  of  the  same  methods  used  to  study  conventional  molecules.[l] 
The  ability  of  muons  to  catalyze  nuclear  fusion  will  be  discussed  in  great  detail  later  in 
this  document. 

While  muon  scattering  experiments  using  both  negative  and  positive  muons  can 
be  perfonned,  those  using  positive  muons  are  more  common.  Material  density  as  a 
function  of  position  can  be  studied  by  comparing  the  energy  and  angle  of  incoming  and 
outgoing  muons.  Cosmic  muons  penetrate  material  so  well  that  it  is  possible  to  use  them 
to  generate  3-D  density  diagrams  of  objects  as  large  as  the  Mayan  pyramids. [13;  14;  15; 
16;  17;  18;  19]  For  obvious  reasons,  interest  has  been  generated  in  using  muons  to  search 
cargo  ships  and  buildings  for  contraband. [20;  21;  22;  23] 

In  addition  to  analyzing  material  by  scattering  muons  through  it,  slow  ( i.e low 
energy)  muons  can  be  absorbed  by  material,  forming  muonic  molecules.  These  muonic 
molecules  can  be  studied  by  many  of  the  same  methods  used  to  study  conventional 
molecules.  By  studying  the  properties  of  these  muonic  molecules,  it  is  often  possible  to 
determine  infonnation  about  the  parent  molecule,  which  absorbed  the  muon.  As  an 
example,  when  forming  muonic  molecules,  negative  muons  will  most  often  initially  form 
a  compound  in  which  the  principle  quantum  number  (n)  is  very  large  (i.e.,  >14).  As  the 
muon  loses  energy,  transitioning  from  one  muonic  state  to  another,  characteristic  x-rays 
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are  emitted  which  can  be  used  to  identify  the  atoms  and/or  molecules  to  which  they  are 
bound. 

Perhaps  the  most  studied  method  of  using  muons  to  characterize  matter  is  muon 
spin  rotation/relaxation/resonance  (juSR ).  A  good  review  of  this  method  can  be  found  in 
Introductory  Muon  Science,  by  Kametada  Nagamine. [4: 1 00- 1 66;  24] 

2.2  Muon-Catalyzed  Fusion 

2.2.1  History 

Muon-catalyzed  fusion  was  first  proposed  by  F.  Charles  Frank  in  1947. [25]  The 
concept  was  simple;  if  a  heavy,  negatively  charged  particle  replaces  an  electron  in  a 
hydrogen  molecule,  the  heavy  negative  particle  will  spend  most  of  its  time  between  the 
nuclei  pulling  them  together  and  shielding  them  from  each  other.  If  one  or  both  of  the 
closely  bound  hydrogen  nuclei  are  heavy  isotopes  of  hydrogen  ( i.e .,  deuterium  and/or 
tritium),  this  process  would  allow  the  nuclei  to  get  so  close  together  that  nuclear  fusion 
could  occur. 

Unaware  of  Charles  Frank’s  prediction  that  muon-catalyzed  fusion  could  occur;  in 
1948  Andrei  Sakharov  made  the  same  prediction  and  estimated  the  fusion  rate  for  dpit 
fusion,  where  d  and  t  represent  deuterium  and  tritium  nuclei  respectively.  [26]  The 
symbol  pi  is  used  to  represent  negative  muons.  Due  to  the  short  half-life  of  muons,  (i.e., 
~2.2xl0’6  s)  he  estimated  that  the  average  number  of  fusions  catalyzed  per  muon  would 
be  slightly  greater  than  one. 

About  10  years  later  Luis  Alvarez  from  the  University  of  California,  Berkeley  ran 
some  experiments  in  which  he  observed  muon-catalyzed  fusion.  At  first  Luis  Alvarez 
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and  his  colleagues  were  very  excited.  They  thought  they  had  solved  the  world’s  energy 
problems.  A  few  quick  calculations,  however,  convinced  them  that  the  energy  released 
by  the  process  is  much  less  than  the  energy  necessary  to  produce  muons. [27]  For  the 
next  several  years,  muon-catalyzed  fusion  was  looked  at  as  an  interesting  phenomenon 
that  was  unlikely  to  result  in  practical  applications. 

In  1977  Semen  Gershtein  from  the  Institute  of  High  Energy  Physics  in 
Serpukhov,  USSR  and  Leonid  Ponomarev  from  the  Joint  Institute  for  Nuclear  Research, 
Dubna,  USSR  predicted  that  d/it  fusion  could  fonn  via  a  resonant  process  in  which  the 
energy  released  on  binding  could  be  divided  between  molecular  vibrations  and  rotations. 
This  process  was  predicted  to  enhance  the  d[it  formation  rate. [28] 

Steve  Jones,  who  worked  at  the  Idaho  National  Engineering  Laboratory  (INEL) 
read  the  predictions  of  Gerstein  and  Ponomarev  and  was  eager  to  test  their  hypothesis.  In 
1982  this  desire  was  realized  at  the  Los  Alamos  Meson  Physics  Facility  (LAMPF).  The 
results  of  Jones’s  and  coworkers’  experiments  showed  that  the  dfit  formation  rate  and  the 
d-t  fusion  rate  were  even  higher  than  Gerstein  and  Ponomarev  predicted. [29]  This 
discovery  generated  great  excitement  amongst  the  researchers. [30] 

Over  the  next  several  years  the  single-muon,  catalyzed  fusion  process  was  studied 
in  great  detail.  It  was  shown  that  muon-catalyzed  fusion  could  occur  at  temperatures  as 
low  as  4  K.[31]  Yields  as  high  as  150  fusions  per  muon  were  observed  in  the 
laboratory.  [32] 

Despite  all  of  the  positive  results,  there  was  one  overwhelming  factor  that 
continued  to  limit  obtainable  yields.  Some  of  the  muons  would  stick  to  helium  nuclei 
fonned  during  fusion  and  remain  attached  until  they  decayed.  [3  3]  Attempts  have  been 
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made  to  find  ways  to  strip  these  muons  from  the  helium  nuclei,  but  so  far  no  satisfactory 
solution  to  this  problem  has  been  found.[31;  34;  35;  36;  37] 

As  of  the  writing  of  this  document,  no  attempt  to  optimize  the  reaction 
temperature,  which  is  predicted  to  be  somewhere  between  1200  K  and  4000  K,  nor  to 
study  the  process  at  reasonably  achievable  high  gas  pressure  has  been  made.  [7]  The 
fusion  yield  under  optimized  conditions  is  calculated  to  be  around  300  fusions  per 
muon.[38;  39]  This  is  short  of  the  approximately  1000  fusions  per  muon  that  is  necessary 
for  this  process  to  be  a  practical  source  of  energy.  [3 8;  40;  41;  42;  43;  44] 

While  there  are  ways  to  increase  the  muon-catalyzed  fusion  reaction  rate  ( e.g ., 
increasing  the  pressure  or  optimizing  the  reaction  temperature),  and  thereby  the  fusion 
yield;  there  are  two  additional  factors  that  need  to  be  addressed  if  single-muon  catalyzed 
fusion  is  going  to  be  useful  as  an  energy  source:  (1)  the  energy  required  to  produce 
muons,  and  (2)  the  number  of  muons  stuck  to  helium  nuclei.  It  is  generally  held  that  one 
or  both  of  these  factors  must  be  addressed  in  order  for  single-muon  catalyzed  fusion  to 
become  a  practical  source  of  energy. 

2.2.2  Single-Muon,  Catalyzed  Fusion  Process 

When  a  muon  is  inserted  into  a  reaction  chamber  containing  protium  ('ll  or  P), 
deuterium  (~H  or  D),  and/or  tritium  (  //  or  T),  a  muonic  atom  (i.e.,  pp,  dp  or  Ip)  is  fonned 
where  p,  d  and  t  represent  the  nuclei  of  protium,  deuterium  and  tritium  respectively. 
When  the  muonic  atom  is  first  fonned  it  is  in  an  excited  state  with  a  principle  quantum 
number  (n)  approximately  equal  to  14  (see  Chapter  7  of  this  document  or  reference  [45]). 
The  excited  muonic  hydrogen  atom  loses  energy  via  x-ray  emission.  While  in  an  excited 
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state  the  muon  may  transfer  to  other  hydrogen  nuclei  via  exchange  reactions.  When  the 
muonic-atom  reaches  or  approaches  the  ground  state  it  can  participate  in  additional 
reactions  which  fonn  tri-nuclear  molecules  (see  Appendix  A). 

When  a  muonic  hydrogen  atom  in  or  near  its  ground  muonic  state  ( i.e .,  n  =  1), 
collides  with  a  diatomic  hydrogen  molecule  a  tri-atomic  muonic  molecule  sometimes 
forms.  For  example, 


p/j  +  DT - >(pd  p^tee  (2.3) 

Two  of  the  nuclei  will  be  bound  very  close  to  each  other  (i.e.,  ~  0.005  A)  by  the  muon. 
The  third  nucleus  is  bound  to  the  first  two  with  a  bond  length  approximately  equal  to  the 
nuclear  bond  length  of  a  diatomic  hydrogen  ion  H*  . 

For  some  of  the  tri-atomic  muonic  molecules  [i.e.,  (ddp)xee,  (dip)xee ,  and 

( ttp)xee ,  where  x  represents  the  nuclei  of  any  hydrogen  isotope]  bound  excited 

vibrational  and  rotational  states  between  the  two  tightly  bound  atoms  exists  (e.g.  v=l, 
J=  1).  The  small  amount  of  binding  energy  present  in  these  excited  exotic  molecules  is 
distributed  between  vibratonal  and  rotational  energy.  This  process,  known  as  resonant 
formation,  is  the  reason  some  muon-catalyzed  fusion  reactions  occur  much  more  rapidly 
than  originally  predicted.  Most  of  the  time  this  weakly  bound  tri-atomic  muonic 
molecule  will  quickly  dissociate,  the  muon  staying  with  the  heavier  of  the  two  nuclei  to 
which  it  is  closely  associated.  For  example, 
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Pju  +  D2 


ju^dee - >dp  +  HD 


(2.4) 


Some  of  the  tri-atomic  muonic  molecules  will  de-excite  by  Auger  de-excitation  rather 
than  dissociate,  resulting  in  a  more  tightly  bound,  more  stable  compound.  For  example, 

i^pd  p^dee - >[pd  p}de  +  e~  (2.5) 

Once  the  tri-atomic  muonic  molecule  is  formed,  the  nuclei  bound  together  by  the 
muon  can  fuse;  if  the  isotope  pair  is  one  for  which  fusion  is  possible.  Although 
rotationally  and  vibrationally  excited  muonic  hydrogen  molecules  sometimes  fuse,  it  is 
more  likely  that  this  will  happen  when  the  molecule  is  in  its  ground  state. 

Once  fusion  has  taken  place,  the  catalytic  muon  can  be  freed,  allowing  it  to 
catalyze  another  fusion  event,  or  it  can  stick  ( i.e .,  bind)  to  one  of  the  particles  fonned 
during  fusion.  When  negative  muons  stick  to  a  fusion  product  they  most  commonly  stick 
to  the  product  of  highest  charge  ( e.g 4 He).  When  fusion  results  in  multiple  products  of 
the  same  charge  being  fonned,  such  as 

dpd - >p  +  t  +  p  (2.6) 

if  the  muon  sticks  to  one  of  the  product  nuclei,  it  will  most  often  stick  to  the  fusion 
product  of  greatest  mass  (e.g.,  t). 
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The  probability  of  a  muon  sticking  to  a  fusion  product  during  the  fusion  process 
depends  not  only  on  the  charge  and  mass  of  the  fusion  products,  but  also  on  the  energy 
and  isotopic  makeup  of  the  tri-atomic  species  that  fuse.  Sticking  may  be  thought  of  as 
being  due  to  the  matching  of  the  muon  kinetic  energy  in  the  initial  state  with  that  of 
bound  final  states  on  the  recoiling  daughter  nucleus.  The  closer  these  energies  match,  the 
higher  the  probability  of  sticking.  [46]  The  effective  probability  of  a  muon  sticking  to  a 
helium  nucleus  during  dpt  fusion  (i.e.,  0)s  -  0.0043  )  is  much  less  than  the  probability  of  it 

sticking  during  ppd  (i.e.,  copd  =  0.99  ),  ppt  (i.e.,  copt  =  0.94  ),  dpd  (i.e.,  cod  -  0.122  ),  or  tpt 
(i.e.,  cot  —  0.14  )  fusion. [7;  10;  12;  31:47;  48;  49;  50;  51;  52;  53;  54;  55;  56;  57;  58]  dpt 
forms  more  rapidly  and  fuses  quicker  than  dpd  and  tpt,  therefore  the  reaction  paths 
leading  to  dpt  fusion  have  predominated  in  the  study  of  muon-catalyzed  fusion.  [7;  8;  9; 
10;  47] 

ppd  and  ppt  do  not  have  bound  excited  vibrational  states; [4:76-77]  therefore 
fusion  paths  which  lead  to  the  fusion  of  these  isotope  pairs  are  much  slower  than  those  in 
which  exclusively  heavy  isotopes  of  hydrogen  fuse. 

Figures  2-1  to  2-5  depict  the  most  common  muon-catalyzed  fusion  reaction  paths. 
The  diagrams  do  not,  however,  show  all  of  the  reactions  which  occur.  Appendix  A  gives 
a  more  complete  list  of  single  muon  catalyzed  fusion  reaction  steps  and  rate  equations 
which  facilitate  the  study  of  muon-catalyzed  fusion  kinetics. 
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2.3  Muon-Catalyzed  Fusion  Yield 


The  initial  impetus  for  most  scientists  pursuing  muon-catalyzed  fusion  was  a  hope 
that  it  would  provide  a  cheap,  clean,  and  abundant  source  of  energy.  Due  to  a-sticking 
(i.e.,  muons  sticking  to  4 He  nuclei),  the  energy  required  to  produce  muons,  and  to  a  lesser 
degree  the  tri-atomic  fonnation  rates  (e.g.,  ddud),  the  goal  of  using  muon-catalyzed 
fusion  as  a  clean  and  efficient  energy  source  has  not  been  realized. 

In  order  for  single-muon  catalyzed  fusion  to  become  an  efficient  source  of  energy, 
the  energy  required  to  produce  muons  must  decrease,  the  number  of  fusions  catalyzed  per 
muon  must  increase,  or  a  hybrid  reactor  must  be  used.  Research  in  all  three  of  these 
directions  has  proceeded,  with  some  success.  There  has  not,  however,  been  sufficient 
success  for  a  pure  fusion  reactor  based  on  muon-catalyzed  fusion  to  be  feasible. 

In  1978  Yu.  Petrov  presented  the  idea  of  a  hybrid  fusion-fission  reactor  based  on 
muon-catalyzed  fusion.[59;  60]  It  is  generally  believed  that  using  existing  technology, 
such  a  reactor  could  be  built.  In  addition  to  energy  production,  muon-catalyzed  fusion 
could  be  used  to  produce  reactor  or  weapons  grade  plutonium-239  (see  Section  2.5). 

2. 3. 1  Muon  Production 

If  the  energy  required  to  produce  muons  is  decreased  significantly,  the  feasibility 
of  muon-catalyzed  fusion  as  a  pure  fusion  source  of  energy  could  change.  A 
considerable  amount  of  research  aimed  at  reducing  the  energy  required  to  produce  muons 
continues,  and  some  improvements  have  been  attained;  yet  not  enough  for  a  pure  fusion 
reactor  to  be  designed. [  1:17-39]  While  muon  production  is  an  important  and  interesting 
area  of  research,  it  is  not  dealt  with  in  any  detail  in  this  document.  This  document 
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emphasizes  the  study  of  mechanisms  that  have  the  potential  of  increasing  the  number  of 
fusions  per  muon  obtainable. 

2. 3.2  Increasing  the  Number  of  Fusions  per  Muon 

There  are  three  factors  that  limit  the  average  number  of  fusions  muons  can 
catalyze:  1)  the  muon  half-life,  2)  the  muonic  molecular  reaction  rates,  and  3)  the 
probability  that  muons  are  pulled  out  of  a  reaction  cycle  through  a-sticking  or  muon 
scavenging  by  non-reactive  molecules  (The  terms  a-sticking  and  scavenging  will  be 
defined  in  sections  2. 3. 2.2  and  2. 3.2. 3  respectively).  There  is  no  proven  way  of 
significantly  affecting  the  muon  half-life,  but  there  are  things  that  can  be  done  to  affect 
the  other  two  factors. 

2. 3. 2.1  Muonic  Molecular  Reaction  Rates.  The  slowest  steps  in  the 
muon-catalyzed  reaction  cycles  are  the  fonnation  rates  of  ppd,  p/it,  dpd,  dpt,  and  tpt  (see 
figures  5-1  through  5-5).  As  was  mentioned  in  section  2.3.2,  due  to  the  existence  of 
bound  excited  vibratonal  and  rotational  states,  resonance  stabilization  occurs.  This 
results  in  dpd,  dpt,  and  tpt  forming  more  rapidly  than  ppd,  and  ppt  which  do  not  have 
bound  excited  vibrational  states. [61]  In  order  to  significantly  increase  the  muonic 
molecular  reaction  rates  the  rate  these  muonic-molecules  fonn  must  increase,  or  a  novel 
reaction  path  must  be  followed. 

In  actuality,  the  muonic  hydrogen  molecules  ppd,  ppt,  dpd,  dpt,  and  tpt  shown  in 
Figures  2-1  through  2-5  are  not  normally  fonned  in  low-temperature  muon-catalyzed 
fusion  (see  Appendix  A).  What  is  fonned  are  tri-nuclear  species  in  which  two  hydrogen 
nuclei  are  closely  muonically  bound  (e.g.,  0.005  A)  and  a  third  hydrogen  atom  is 
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electronically  bound  to  the  first  two  with  a  bond  distance  of  approximately  1.06  A.  For 
example, 


tju  +  D2 


■> 


e  +  e 


(2.7) 


Although  the  electronically  bound,  nucleus  does  not  have  a  significant  impact  on 
the  rate  at  which  the  closely  bound  hydrogen  nuclei  fuse,  it  does  have  a  significant  impact 
on  the  formation  rate  of  the  tri-nuclear  molecular  species  ddup,  ddjud,  ddpt,  dtjup,  dtpd, 
dtpt,  ttpp,  ttpd,  and  ttpt  and  can  impact  the  rate  which  the  ground  rotational  and 
vibrational  states  are  reached.  This  impact  is  due  to  the  existence  of  weakly  bound 
rotational  and  vibrational  states.  The  energy  released  on  binding  can  be  distributed 
between  vibrations  and  rotations.  [62] 

23.2.2  Muon  Sticking  to  Fusion  Product.  During  muon-catalyzed  fusion, 
muons  will  sometimes  bind  (i.e.,  stick)  to  nuclei  fonned  during  fusion.  When  fusion 
results  in  products  of  different  charge  being  formed,  if  the  muon  sticks  to  one  of  the 
products,  it  will  almost  always  be  to  the  product  of  greatest  charge.  When  nuclei  of 
different  mass,  but  the  same  charge  are  fonned,  the  muon  will  most  often  stick  to  the 
heavier  product;  if  it  sticks  at  all.  It  is  currently  believed  that  the  probability  a  muon  will 
stick  to  the  lighter  fusion  product  is  so  small  that  it  can  be  neglected  in  describing  muon- 
catalyzed  fusion  mechanisms. [47] 

When  fusion  results  in  the  formation  of  a  muonic  atom,  or  a  muonic  atomic  ion, 
the  muonic  atom  or  ion  will  be  in  an  excited  muonic  state.  After  fonnation  of  the  muonic 
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atom  or  ion,  collisions  can  cause  the  release,  or  transfer  of  the  muon.  This  process  is 
referred  to  as  stripping.  When  a  stuck  muon  is  transferred  to  a  molecule  that  allows  it  to 
continue  catalyzing  nuclear  fusion,  the  process  is  known  as  regeneration.  Stripping 
occurs  readily  when  the  muonic  atom  or  ion  is  in  a  highly  excited  muonic  state.  It  is  less 
likely  to  occur  when  the  atom  or  ion  is  in  a  low  excited  state,  or  in  its  ground  state.  The 
probability  that  a  muon  will  stick  to  a  fusion  product  (oj)  minus  the  probability  that  it  will 
be  regenerated  (R)  is  referred  to  as  an  effective  sticking  constant  (co°). 

(0°  =  co- R  (2.8) 

Being  that  the  probability  of  regeneration  is  dependent  upon  the  probability  of 
collisions  occurring  and  on  the  energy  of  the  collisions;  there  must  be  some  density  and 
temperature  dependence  on  the  effective  sticking  constant.  As  temperature  and  pressure 
increase  the  effective  sticking  probability  decreases. [43;  50;  63;  64;  49] 

Many  methods  of  enhancing  regeneration  have  been  considered,  but  as  of  the 
writing  of  this  document,  no  effective  method  of  preventing  sticking  or  of  regenerating 
muons  stuck  during  fusion  has  been  demonstrated. [63;  65;  64;  43;  66;  67;  68;  69] 
Regeneration  in  electric  fields  has  been  considered  and  may  hold  some  promise.  [70] 

2. 3. 2. 3  Muon  Scavensins.  In  addition  to  muons  sticking  to  fusion 
products,  they  can  be  removed  from  the  catalytic  fusion  cycle  by  binding  to  non-reactive 
atoms  within  the  reaction  chamber  or  which  make  up  the  chamber  walls.  This  process  is 
referred  to  as  scavenging.  The  probability  of  muon  scavenging  occurring  is  directly 
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proportional  to  the  probability  of  muons,  muonic  atoms,  and  muonic  ions  colliding  with 
atoms  incapable  of  undergoing  low  temperature  fusion.  Therefore,  the  pressure, 
temperature,  volume  and  inside  surface  area  of  the  reaction  chamber  affects  the 
scavenging  rate  and  probability.  Additionally,  the  concentration,  nuclear  charge,  and 
mass  of  non-reactive  atoms  within  the  chamber  affect  the  scavenging  process  and  rate. 

2.4  Applications  of  Muon-Catalyzed  Fusion  with  Currently  Obtainable  Yields 

While  a  pure  fusion  energy  source  is  the  goal  most  investigators  studying  muon- 
catalyzed  fusion  have  pursued,  it  is  not  the  only  application  of  the  process.  Each  d-t 
fusion  produces  a  14.1  MeV  neutron  and  about  17.6  MeV  of  total  energy.  If  a  Uranium- 
238  (U-238)  blanket  is  placed  around  the  reaction  chamber,  each  of  these  neutrons  could 
initiate  U-238  fission,  releasing  an  additional  200  MeV  of  energy.  This  may  well  be  high 
enough  for  muon-catalyzed  fusion  to  be  considered  a  practical  energy  source  using 
current  technology.  During  the  fission  process  an  average  of  about  4.1  neutrons  per 
fission  is  produced.  These  neutrons  are  lower  in  energy  than  the  dt  fusion  neutrons  and 
rarely  initiate  further  U-238  fission.  These  neutrons  can,  however,  be  absorbed  by  U- 
238,  producing  U-239  which  decays  to  Pu-239 

n  +  238U  — >  239U  — >  239 Np  +  e~+ve->  239 Pu  +  e~  +  ve 

This  method  could  be  used  to  produce  very  high  purity  Pu-239  which  could  then 
be  used  in  fission  reactors  or  weapons.  Assuming  a  maximum  muon-catalyzed  fusion 
yield  of  300  fusions  per  muon,  each  muon  could  catalyze  the  production  of  1230  atoms  of 


(2.9) 
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Pu-239.  Assuming  the  maximum  negative  muon  flux  to  be  equal  to  that  produced  at  the 
Gatchina  synchrocyclotron  (1.5  x  10 16  p'/s)  [1:34]  where  the  scientist  who  first  proposed 
this  method  worked, [59;  60;  47]  an  upper  limit  for  plutonium  production  is  231  kg/year. 
The  total  number  of  fissions  resulting  from  one  negative  muon  is  about  1500  fisions  (i.e., 
considering  both  U-238  fusion  and  Pu-239).  The  total  energy  produced  would  be  greater 
than  300  GeV  per  muon.  These  numbers  represent  an  upper  limit  on  the  energy  and 
number  of  fissions  that  could  be  produced,  but  even  considering  much  lower  efficiency, 
the  usefulness  of  this  process  can  be  seen. 

In  addition  to  energy  and  Pu-239  production,  there  is  another  practical 
application  of  muon-catalyzed  fusion.  It  can  be  used  to  produce  slow,  mono-energetic 
muons  which  are  useful  in  many  types  of  experiments. [4:32]  Most  sources  of  muons 
result  in  a  wide  spectrum  of  muons  being  produced.  If  moderated  muons  are  shot  into 
hydrogen,  the  muons  released  after  fusion  will  be  nearly  mono -energetic.  The  energy  of 
these  mono-energetic  muons  will  depend  on  which  muon-catalyzed  fusion  reaction 
occurs,  which  depends  on  the  isotopic  mix  of  the  reacting  hydrogen. 

2.5  Conclusions 

In  order  for  single-muon  catalyzed  fusion  to  be  a  useful  source  of  energy,  the 
energy  required  to  produce  muons  must  decrease,  the  number  of  fusions  per  muon  must 
increase,  or  a  hybrid  reactor  ( e.g fusion-fission)  must  be  used. 

Advances  in  methods  of  muon  production  have  been  made,  but  the  energy  cost 
per  muon  is  still  too  high  for  a  pure  fusion  reactor  based  on  currently  tested  muon- 
catalyzed  fusion  reactions. [4: 17-39]  While  the  design  of  a  feasible  pure  fusion,  muon- 
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catalyzed  fusion  reactor  is  likely  a  long  way  off,  viable  hybrid  fusion-fission  reactors 
could  be  built  using  existing  technology. 
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III.  Novel  Muon-Catalyzed  Fusion  Reaction  Paths 


In  this  chapter  three  novel  muon-catalyzed  fusion  reaction  paths  which  have  the 
potential  to  increase  the  number  of  fusions  per  muon  that  can  be  catalyzed  will  be 
presented:  1)  replacing  the  weakly  bound  hydrogen  nuclei  in  tri-nuclear  muonic 
molecules  with  a  positive  muon  or  a  positron,  2)  using  di-muonic  reactions  to  free  stuck 
muons,  and  3)  using  di  muonic-reactions  to  increase  the  formation  rate  of  tri-nuclear 
muonic  hydrogen.  The  second  and  third  of  these  novel  reaction  paths  will  be  addressed 
in  detail  in  Chapters  6  through  8  of  this  document. 

3.1  Addition  of  Positive  Muons  and  Positrons 

During  muon-catalysis  of  hydrogen  fusion,  molecules  of  the  form  (xjuy)zee  are 

fonned  where  jc,  y  and  z  represent  hydrogen  isotopes,  p  and  e  represent  negative  muons 
and  electrons  respectively,  x,  y  and  z  may  or  may  not  be  equivalent.  The  equilibrium 
bond  length  between  x  and  y  is  about  0.005  A.  The  jc-z  and  y-z  bond  length  is  around 
0.74  A  (see  Chapters  4  and  7). 

The  rate  at  which  the  isotopes  x  and  y  fuse  is  affected  very  little  by  the  mass  of  z; 
however,  the  rate  of  (xyp)zee  formation  is  strongly  influenced  by  the  mass  of  z  (see 

Table  3-1). [1;  2;  3;  4;  5].  When  z  represent  the  three  hydrogen  isotopes  p,  d,  and  t  (i.e., 
H,  ~H,  and  H  respectively)  the  lighter  the  mass  of  z,  the  faster  the  fonnation  rate.  From 
this,  the  question  arises,  what  happens  if  the  mass  of  z  is  decreased  even  further  by 
substitution  with  a  positive  muon  or  a  positron?  The  binding  energy  will  decrease  as  a 
result  of  the  smaller  mass.  In  the  case  of  a  positive  muon  being  added  to  the  system,  the 
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binding  energy  will  increase  as  a  result  of  correlation  energy  between  the  oppositely 
charged  particles  (the  closer  the  mass  of  the  particles  is  to  each  other,  the  greater  the 
correlation  energy).  In  the  case  of  a  positron,  the  correlation  energy  would  further 
decrease  due  to  the  relative  difference  in  mass  between  a  negative  muon  and  a  positron. 
Which  of  these  factors  predominates  and  how  they  affect  the  excited  state  binding  energy 
and  molecular  fonnation  rate  is  an  open  question.  The  author  is  unaware  of  any 
calculations  or  experiments  having  been  perfonned  to  address  these  questions. 


Table  3.1.  Formation  rate  (X)  of  ground 
vibrational  state  (dtp)d  and  (dtp)t  at  300  K.  Values 


are  experimental,  taken  i 

from  reference. [6;  1] 

Molecule 

(dtfi)p 

2xl010 

(dtju)d 

4xl08 

(dlfi)t 

2xl07 

A  Novel  Approach  to  Decrease  Muon  Loss  Due  to  Sticking 

One  of  the  products  of  p-d  and  d-d  fusion  is  He;  some 

catalyze  fusion  reactions  will  bind  to  a  3 He  nucleus: 


pd/up - >y  +  p^He+H' 


(3.1) 
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pd/ud - >y  +  /j^He  +  D' 


(3.2) 


ddjup - >  p3He  +  H*  +  nx 


(3-3) 


dd  pd - >ju3He  +  D *  +  nx 


(3-4) 


This  process  is  referred  to  as  muon  sticking.  In  the  case  of  d — d  fusion  the  probability  of 

•5 

the  muon  sticking  to  He  is  about  12%.  An  effective  method  of  stripping  the  muons  from 
the  'He  nucleus  so  that  the  muons  can  continue  to  catalyze  fusion  events  has  so  far  not 
been  demonstrated. 

An  alternative  to  stripping  the  muon  from  the  'He  nucleus  is  to  fuse  the  He  via 
the  reactions: 


D+3He - >4He  +  p 


(3.5) 


or 


3  He +  3  He - >4He  +  2p 


(3-6) 


The  problem  with  this  approach  is  that  the  Coulombic  repulsion  between  these  nuclei  is 
greater  than  the  repulsion  between  two  hydrogen  nuclei.  As  a  result,  a  single  muon 
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doesn’t  bring  the  nuclei  close  enough  together  for  them  to  fuse.  What  had  not  been 
studied,  prior  to  the  results  presented  in  this  document,  is  what  happens  when  there  are 
two  negative  muons  participating  in  the  reactions.  For  example: 


dp  +  p^He - >dp4Hep - >4He  +  p  +  2p  (3.7) 


and 


p3He  +  p^He - >p^He p  }He - >4He+2p  +  2p  (3.8) 


(see  Figures  3-1  and  3-2).  These  reactions  will  be  studied  in  detail  and  presented  in 
Chapter  6. 
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Figure  3-2  Reaction  cycle  for  using  3 Hep-3 Hep  fusion  to  regenerate  //'  during  dpd  muon- 
catalyzed  fusion. 


3.3  Two-Muon  Catalyzed-Fusion 

The  slowest  step  in  a  homogeneous  single-muon  catalyzed  fusion  reaction  chain  is 
often  considered  to  be  the  formation  of  tri-nuclear  molecules  that  result  in  fusion  (see 
Figures  2-1  to  2-5  and  Appendix  A).  In  actuality,  this  muonic  molecular  system  forms 
readily;  however,  it  most  often  forms  in  a  rotationally  and  vibrationally  excited  state  that 
is  only  slightly  bound.  In  excited  states,  the  molecules  are  much  more  likely  to  dissociate 
than  to  fuse,  or  fall  to  a  more  stable  vibrational  and  rotational  state  which  later  fuses.  The 
formation  rate  most  often  reported  in  the  literature  is  not  the  average  rate  of  fonnation; 
rather  it  is  an  effective  rate.  It  represents  the  average  rate  of  formation  of  a  tri-nuclear 
molecule  which  later  fuses  ( i.e those  which  dissociate  are  neglected). [6;  8] 

Prior  to  this  work,  it  was  not  certain  what  would  happen  if  exotic  hydrogen 
molecules  containing  two  negative  muons  were  formed  (see  Chapters  4,  5  and  7).  It  was 
predicted  that  the  bond  distance  between  the  hydrogen  molecules  would  be  shorter, 
resulting  in  a  system  which  fuses  more  rapidly,  possibly  from  an  excited  state.  [9]  Some 
di-muonic  hydrogen  molecules  are  bound  more  strongly  than  single  muon  molecules  and 
have  richer  rotational  manifolds.  Figures  3-3  and  3-4  gives  examples  of  how  the 
presence  of  di-muonic  hydrogen  molecules  could  change  a  muon-catalyzed  fusion 
reaction  cycle. 
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Figure  3-4  Reaction  cycle  for  t[it+  and  t ft  fit  fusion. 

3.4  Enhance  Fusion  Yields  with  Electromagnetic  Radiation 

Spatially  coherent  electromagnetic  radiation  could  be  used  to  selectively  break 
apart  undesirable  molecules  within  a  fusion  reaction  chamber.  [10;  11]  An  obvious  place 
where  this  could  be  applied  is  in  freeing  muons  bound  to  nonreactive  nuclei.  Not  having 
a  method  of  efficiently  regenerating  muons  bound  to  helium  nuclei  is  one  of  the  major 
factors  preventing  muon-catalyzed  fusion  from  being  an  efficient  stand  alone  source  of 
energy.  In  addition  to  freeing  bound  muons,  electromagnetic  radiation  could  be  used  to 
dissociate  undesirable  hydrogen  molecules.  An  example  of  when  this  would  be  desirable 
is  when  deuterium  and  tritium  are  both  present  in  a  reaction  chamber.  At  optimum 
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muon-catalyzed  fusion  reaction  temperatures  the  reaction  rate  of  DT  with  muonic 
hydrogen  atoms  is  slower  than  the  reaction  rates  of  D2  and  T2.  Due  to  the  presence  of 
tritium  betas,  at  optimum  reaction  temperatures,  mixtures  of  D2  and  To  quickly  reach 
equilibrium  with  DT.  The  use  of  photons  to  selectively  dissociate  DT  molecules  could 
have  benefits.  The  benefits  of  using  electromagnetic  radiation  in  this  manner  need  to  be 
weighed  against  the  energy  cost  of  its  production. 

While  the  use  of  lasers  to  enhance  muon-catalyzed  fusion  yields  is  worthy  of 
study;  it  is  not  the  only  method  of  decreasing  the  concentration  of  undesirable  molecules 
in  reaction  chambers,  or  of  freeing  bound  muons.  For  example,  synthetic  zeolites  could 
be  used  to  separate  molecular  isotopes  of  hydrogen  [12;  10]  as  could  thermal  diffusion 
and  gas  chromatography.  [10] 

3.5  Conclusions 

There  are  several  novel  reaction  paths  that  have  the  possibility  of  increasing  the 
number  of  fusions  catalyzed  per  muon.  This  includes  adding  positive  muons  or  positrons 
to  the  fusing  muonic  molecules  and  using  multiple  negative  muons  to  increase  the 
reaction  rate.  Two  methods  of  using  multiple  negative  muons  to  increase  the  number  of 
fusions  per  muon  that  can  be  obtained  were  discussed  in  this  chapter.  The  majority  of 
this  dissertation  is  devoted  to  the  study  of  these  methods. 
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IV.  General  Particle  Orbital  Method  of  Modeling 
Molecules  Made  of  Multiple  Types  of  Quantum  Particles 

4.1  Introduction 

In  order  to  model  physical  properties  of  muonic  hydrogen  and  helium  molecules  a 
General  Particle  Orbital  (GPO)  method  of  non-adiabatic  quantum  mechanics  was 
developed.  The  GPO  method  is  based  on  non-adiabatic  Hartree-Fock — configuration 
interaction  (HF  /Cl)  methods.  It  facilitates  the  modeling  of  particles  with  multiple  types 
(i.e.,  mass  and  charge)  of  quantum  particles.  Although  the  method  was  developed  to 
study  muonic-molecules,  its  usefulness  is  not  limited  to  molecules  containing  muons. 
Other  types  of  exotic  molecules  and  conventional  molecules  with  some  or  all  of  the 
nuclei  being  modeled  as  quantum  particles  can  be  studied. 

Non-adiabatic  methods  of  quantum  mechanics  can  be  used  to  study  the  properties 
of  molecular  systems  in  which  nuclei  are  not  accurately  modeled  as  fixed  points  under 
the  Born-Oppenheimer  approximation  (e.g.,  see  References  [1]  through  [11]).  Methods 
capable  of  modeling  non-adiabatic  systems  are  particularly  useful  when  studying 
molecular  systems  containing  exotic  particles  in  which  more  than  one  type  of  particle 
must  be  modeled  quantum  mechanically. [12;  13;  14;  15;  16]  Due  to  the  computational 
complexity  of  non-adiabatic  methods,  their  use  has  traditionally  been  limited  to 
comparatively  small  molecular  systems. 

Nuclear  Electronic  Orbital  (NEO)  methods  were  developed  in  order  to  facilitate 
the  study  of  quantum  chemical  effects  between  protons  and  electrons  in  larger  molecular 
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systems. [2;  17;  18;  7;  8;  9;  10;  11]  These  NEO  methods  treat  electrons  and  selected 
protons  as  quantum  particles,  while  treating  other  nuclei  classically.  In  2007  these 
methods  were  expanded  to  describe  systems  containing  electrons  and  positrons  in  the 
presence  of  classical  nuclei. [19;  20]  While  it  is  valuable  for  describing  positron 
chemistry,  this  modified  NEO  method  was  limited  to  the  study  of  systems  containing  two 
types  of  quantum  particles,  one  positive,  and  one  negative.  In  this  work,  the  NEO 
method  is  extended  to  model  systems  made  up  of  several  different  types  ( i.e mass  and 
charge)  of  quantum  particles.  Since  the  method  described  in  this  paper  is  not  limited  to 
the  chemistry  of  ordinary  molecules  having  electron  and  nuclear  orbitals,  but  can  be 
applied  to  Coulombic  orbitals  for  particles  of  any  mass  and  charge,  this  extension  is 
referred  to  as  a  “General  Particle  Orbital”  (GPO)  method. 

While  the  GPO  method  does  not,  in  principle,  limit  the  number  of  different  types 
of  quantum  particles  that  can  be  described  in  a  molecule,  computational  resources 
available  can  limit  its  applications.  Code  for  its  current  manifestation  is  not  an  efficient 
code  for  multi-processor  computers  nor  does  it  take  full  advantage  of  the  molecular 
symmetry  that  exists  in  many  molecular  systems.  Calculations  involving  more  than  three 
types  of  quantum  particles  are  currently  very  demanding  due  to  the  large  computational 
resources  needed  to  accurately  model  these  systems.  This  is  particularly  true  when 
correlation  energy  is  calculated  by  configuration  interaction  (Cl).  Correlation 
interactions  make  significant  contribution  to  the  molecular  stability  of  the  systems 
described  in  this  paper.  Future  modifications  of  the  code  will  likely  address  its 
deficiencies  and  increase  the  size  of  molecular  systems  which  can  be  accurately  modeled. 
A  computational  strength  of  NEO  and  GPO  methods,  that  can  enhance  efficiency,  is  that 
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some  particles  in  modeled  molecular  systems  can  be  treated  quantum-mechanically  while 
other  particles  are  treated  classically. 

The  ability  to  accurately  calculate  physical  properties  of  a  molecular  system  in 
which  some  of  the  particles  are  considered  quantum-mechanically,  and  others  are 
considered  classically  may  not  be  limited  to  systems  having  negative  and  positive 
particles  with  large  mass  differences.  Often  particles  do  not  interact  strongly  due  to  their 
locations  within  a  molecule.  As  a  result,  their  quantum  interactions  with  each  other  can 
be  ignored.  For  example,  in  some  muonic  molecules,  due  to  the  highly  localized  muon 
density,  it  may  be  possible  to  accurately  consider  weakly  interacting  nuclei  classically, 
even  when  their  mass  is  comparatively  close  to  that  of  a  muon. 

Hartree-Fock  and  Cl  molecular  orbital  theories  are  the  basis  for  the  GPO 
extension  of  NEO  presented  in  this  chapter.  For  HF  and  post-HF  methods,  the  choice  of 
basis  sets  and  basis  set  locations  is  of  great  importance  to  the  accuracy  of  calculations; 
this  GPO  method  is  no  exception.  Approaches  and  challenges  of  optimizing  basis  sets 
and  basis  set  center  locations  while  using  a  minimum  of  computational  resources  are 
discussed  in  this  chapter. 

There  are  several  factors  which  are  often  addressed  in  post-HF  methods,  which 
significantly  influence  Coulombic  interactions  in  many  exotic  molecular  systems. 
Among  these  are  correlation  energy  ( i.e .,  especially  correlation  between  different  types  of 
particles),  relativistic  effects,  nuclear  volume  effects,  vacuum  polarization,  and  charge 
density  effects  (i.e.,  effects  of  nuclear  electromagnetic  structure  and  nuclear 
polarization).[21;  22;  23;  24;  25;  26]  The  relative  importance  of  each  of  these  factors 
depends  on  the  particular  system  being  studied.  In  this  work  only  correlation  energy, 


48 


which  is  arguably  the  most  important  of  these  factors  in  low-Z  muonic  molecules,  is 
considered. 

While  there  are  many  methods  of  calculating  the  correlation  energy  between 
quantum  particles,  [8;  27:193-228,226,383,430,448-449,455;  28:61,231-269,238]  only 
configuration  interactions  (Cl)  will  be  addressed  in  this  paper.  A  Cl  application  of  the 
NEO  method  to  model  non-adiabatic  quantum  mechanics  for  quantum  protons  and 
electrons  was  published  previously.  [2]  In  this  paper  the  method  is  expanded  to  model 
any  number  of  different  types  of  quantum  particles. 

Accurately  calculating  the  correlation  energy  between  quantum  particles  using  Cl 
methods  can  be  computationally  demanding  and  require  significant  resources.  This  is 
particularly  true  when  calculating  the  interactions  between  different  types  of  particles. 
Calculations  in  this  chapter  illustrate  the  GPO  method  using  Cl  to  account  for  correlation 
interactions  between  particles.  For  molecular  systems  containing  negligible  static 
correlation,  it  is  not  necessary  to  account  for  dynamic  correlation,  between  all  particles, 
using  a  single  method.  For  example,  in  a  system  containing  quantum  protons  and 
electrons,  the  electron-electron  correlation  energy  and  the  proton-proton  correlation 
energy  could  be  calculated  by  one  method  ( e.g .,  Cl)  and  the  electron-proton  correlation 
energy  could  be  calculated  by  an  entirely  different  method  (e.g.,  explicit  correlation). [10; 
29;  30] 
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4.2  Theory 


In  this  section,  Hartree-Fock  (HF)  and  configuration  interaction  (Cl)  methods  are 
developed  in  a  formalism  useful  to  study  systems  that  involve  any  number  and  type  of 
quantum  Fennions  in  the  presence  of  classical  particles.  These  methods  are  particularly 
useful  in  the  study  of  exotic  particle  molecular  systems  that  contain  multiple  types  ( i.e 
various  mass  and  charge  combinations)  of  quantum  particles.  The  methods  presented 
were  developed  to  study  Fennions;  however,  they  can  be  used,  without  modification,  to 
study  systems  containing  Bosons,  as  long  as  all  quantum  Bosons  in  the  system  are 
distinguishable  (i.e.,  have  different  mass  and/or  charge). 

4. 2. 1  GPO  Hartree-Fock  method 

In  the  derivations  which  follow  an  open-shell  unrestricted  Hartree-Fock  method  is 
combined  with  a  closed-shell  restricted  method  to  allow  multiple  types  of  quantum 
particles  to  be  modeled  simultaneously  in  a  system  using  whichever  method  is  most 
appropriate  for  each  type  of  particle  ( e.g electrons,  muons,  positrons,  etc.). 

The  Hamiltonian  operator  can  be  written  in  the  form 


H  = 


N  1 

I- 

„  m. 


-V2  + 


N  N 

E2X 

a  P>a 


(4.1) 


where  N  is  the  total  number  of  particles  in  the  system.  Vuj,  can  be  any  function 
describing  a  Coulombic  potential  with  the  distance  between  particles  rufj  and  particles 
charge  Za  and  Z/( .  For  conventional  molecules  it  is  common  practice  to  use  the  Born- 
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Oppenheimer  approximation  and  to  break  the  Hamiltonian  operator  into  one  and  two 
quantum  particle  terms  ( i.e .,  one-electron  and  two-electron  terms).  When  one  of  the 
interacting  particles  is  at  a  fixed  location,  a  one-quantum  particle  term  hap  results.  When 
neither  interacting  particle  can  be  assumed  to  be  fixed,  the  expectation  value  of  the 
particle’s  contribution  to  the  Hamiltonian  is  a  two-quantum  particle  multi-center  integral 

(XtXfiXk x)  =  LLX  Mx]{xfi)Vapzk{xa)x,(,xp)dxadxp  (4.2) 


where  Xt  represents  spin  orbitals. 

If  h  represents  core-Hamiltonian  operators  for  quantum  particles  in  the  field  of 
classical  (i.e.,  fixed-location)  particles,  the  total  energy  of  the  system  [2;  19:27;  28:126] 
is 
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(4.3) 


where  /V  is  the  total  number  of  quantum  particles;  /V,  is  the  number  of  particles  of  type  i; 
and  n  is  the  number  of  types  of  particles  being  considered  quantum  mechanically. 
Therefore, 
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(4.4) 


N=tNi 

i= 1 

It  should  be  noted  that  (xiXf^XkXt)  's  no1  equivalent  to  the  common  notation 
{XiXj  I XkXi)  ■  {XiXj^XkXi)  allows  for  any  charge  {i.e.,  attraction  or  repulsion)  and  Vafi 

does  not  have  to  be  proportional  to  — .  Vap  can  be  any  appropriate  functions  of  ra/} ,  Za 

YaP 

and  Z p .  Potentials  which  are  not  proportional  to  —  are  needed  when  nuclear  volume  is 

V 

included  in  quantum  mechanical  calculations. 

[note:  All  of  the  ab  initio  calculations  presented  in  this  dissertation  use  a  potential 
proportional  to  —  .  A  correction  to  this  potential  is  included  in  some  of  the  calculations 

rafi 

using  a  post-Hartree-Fock  perturbation  method  which  is  presented  in  Chapter  6.  The 
methods  of  including  nuclear  volume  presented  in  Chapter  6  can  only  be  considered 
accurate  for  muonic  molecules  that  have  light  nuclei  {i.e.,  nuclei  from  the  first  few  rows 
of  the  periodic  table).  In  order  to  accurately  model  interactions  of  muons  with  larger 
nuclei  it  is  necessary  to  include  an  appropriate  potential  in  the  ab  initio  calculations.] 

It  is  possible  to  solve  for  the  Hartree-Fock  energy  in  terms  of  closed  or  open  shell 
configurations.  A  closed  shell  configuration  is  one  in  which  all  of  the  occupied  valence 
shells  are  full  (i.e.,  contain  indistinguishable  paired  particles).  An  open  shell 
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configuration  is  one  in  which  some  or  all  of  the  occupied  valence  shells  contain  only  one 
particle.  In  the  discussion  which  follows,  closed  shell  particles  are  modeled  using 
restricted  Hartree-Fock  methods  and  open  shell  particles  are  modeled  using  unrestricted 
Hartree-Fock  methods. [28;  27]  If  the  number  of  restricted  closed  shell  particles  of  type 
cl  is  defined  as  Nci  and  the  number  of  unrestricted  open  shell  particles  of  type  op  is 
defined  as  Nop,  where  cl  =  {\,2,---,nclosed)  and  op  =  {nclosed  +  l,nclosed  +  2, •••,«},  then  the 

Hartree-Fock  energy  of  a  system  containing  any  number  of  any  type  of  restricted  low- 
spin  and/or  unrestricted  high-spin  quantum  particles  is 


N  1  n  NoP  NoP 

£  =  2>, \hW)+-  £  V, 

i= 1  ^  op=(nrln-„j+l)  i= 1  /=1 


opy/op 


°P=(«chsed+l)  i=1  j=1 


r;>7 


Kl  Pa_ 

^ closed  2  2 


ciosea  z.  z.  i—  .  .  . 

+  £  T,T,[2(v'?i'cP  ry 


cl= 1  7=1  7=1 

Nod  N, 


y: 


(4.5) 


n  n  "oP  "d 

jU=l  V>jU  7=1  7  =  1 


where  y/  represents  spatial  orbitals,  n  is  the  number  of  particle  types,  and  nciose(j  is  the 
number  of  particle  types  in  which  all  of  the  particles  of  that  type  reside  in  closed 
shells. [2;  19:27;  28:134] 

Using  the  variational  method  to  minimize  energy,  Fock  operators  can  be 
expressed  as: 
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(4.6) 
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(4.7) 


where  Coulomb  and  exchange  operators  are  defined  as 

J,(l)  =  j^*(2)  Via  y/ j (2) dr2  (4.8) 


and 


^(1)^(1)  =  [J^(2)  Vl2 


'/C(2)  dr2  y/  (1) 


(4.9) 


The  Hartree-Fock-Roothan  equations,  coefficient  matrix  elements  and  charge- 
density  bond-order  matrix  elements  ( PXrJ  )  are  the  same  as  ordinarily  defined. [2;  28:137; 

27]  The  one-quantum  particle  terms  in  the  Fock  matrix  ihap)  are  defined  the  same  for 
restricted  closed  and  unrestricted  open  shell  particle  types. [2;  28:71]  The  two-quantum 
particle  terms  in  the  Fock  matrix  are  not  the  same  for  the  unrestricted  open  shell  case  and 
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the  restricted  closed  shell  case.  If  and  N"''  are  the  number  of  basis  functions 
combined  to  approximate  the  wave  functions  for  particles  of  type  cl  and  op  respectively, 
then  the  two  quantum  particle  terms  Gc‘  and  Gopv  are  defined  as: 


(4.10) 


and 


a Z  =  Z  Z  PZ 

A=\  (7=1 
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<p°>°:x  cpiwt 


(4.11) 


From  these  equations  the  Fock  matrix  elements  ( Ffkn, )  can  be  derived 

n  N‘bf 

Fk  =hk  +Gk  +YYY  Pi 

jliv  [iv  [iv  /  j  /  j  /  j  Act 
l^k  A= 1  <7=1 

where  k  represents  all  types  of  particles,  open  or  closed  shell. 

These  Hartree-Fock-Roothaan  equations  were  solved  iteratively  using 
convergence  accelerators  developed  previously  for  electronic  structure  theory.  [31]  The 
Fock  equations  for  each  particle  were  fully  converged  sequentially  after  each  step  in  the 
iterative  procedure. 


<PWi 


(4.12) 


55 


4.2.2  Configuration  Interaction 


A  configuration  interaction  (Cl)  wave  function  y/total  can  be  expressed  as  a  linear 
combination  of  configuration  functions,  d> .  For  a  system  of  n  types  of  quantum  particles 
this  is 


N'a 

. total  zz- 

*,=u,=i 


k =1 


(4.13) 


where  N‘CI  is  the  number  of  quantum  detenninants  of  particles  of  type  i.  If  N'so  is  the 
total  number  of  spin  orbitals  of  a  given  type  i, 


N'  = 
JVc/ 


N  ' 

so  ■ 


MNso -»,)'■ 


(4.14) 


for  Ni  particles  of  type  i .  The  total  number  of  quantum  detenninants  is 


N'a'1  =  FI  Na  (4.15) 

i=\ 


The  total  Cl  energy  Etotal  =  (wtotal  \  Htotal  \  Wtotal)  is 
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N;  is  the  number  of  particles  of  type  i  and  n  is  the  number  of  quantum  particles  the 
system  contains. 

If  we  define  spatial  shift  operators  Ea  p  between  particles  a  and  [>  of  type  m,  in 
terms  of  spin-creation  ( T  )  and  annihilation  (4<  )  operators  (a), 


then  HCI  can  be  expressed  in  tenns  of  spatial  orbitals  ( i//  ) 
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where  N'MO  is  the  number  of  molecular  orbitals  of  particle  i  in  the  configuration  space. 
From  this  it  follows  that  the  total  energy  of  the  system  is 
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where  the  density  matrix  elements  are  defined  as 
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and 
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If  the  variation  method  is  used  to  minimize  the  total  energy  with  respect  to  the  Cl 
coefficients  C(/t|  A  k  ) ,  then  the  Cl  coefficients  and  corresponding  vibrational-Coulombic 

ground  and  excited  state  energies  may  be  detennined  by  diagonalizing  the  Cl 
Hamiltonian  matrix.  This  method  allows  any  choice  of  multiplicity  for  all  of  the 
Fennionic  nuclei,  thereby  facilitating  the  study  of  systems  of  particles  having  different 
spin  states. [2] 

4.2.3  Basis  Set  Development 

As  originally  developed,  NEO  implemented  post-Hartree-Fock  methods  for 
molecules  having  two  types  of  quantum  particles,  one  positive  and  one  negative. [2;  11; 
18;  29;  30;  32;  33;  34;  35;  36;  37;  38]  Using  these  procedures,  the  total  number  of  self- 
consistent-field  (SCF)  iterations  that  must  be  performed  in  order  to  obtain  density 
convergence  is  the  product  of  the  number  of  iterations  (U)  needed  for  each  type  of 
quantum  particle  (/).  Expanding  the  procedure  to  n  types  of  quantum  particles;  the  total 
number  of  SCF  iterations  needed  to  obtain  density  convergence  is 
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(4.24) 


Total  #  of  SCF  iterations  =  ]^[  Ui 

i=i 

As  a  result,  minimizing  the  time  required  to  perform  individual  SCF  iterations  becomes 
more  important  as  the  number  of  types  of  quantum  particles  in  modeled  systems 
increases.  The  importance  of  having  the  most  efficient  basis  functions  possible  increases 
exponentially  with  the  number  of  quantum  particles,  and  the  relative  importance  of  the 
time  required  to  perform  the  integration  decreases.  Choosing  small,  efficient  sets  of  basis 
functions  to  represent  a  molecular  orbital  necessitates  fewer  HF-SCF  iterations  to  obtain 
converged  density  for  each  type  of  particle.  As  the  number  of  basis  set  functions 
decreases  not  only  does  the  time  required  for  each  SCF  iteration  decrease,  but  most  often, 
the  number  of  iterations  required  usually  decreases  as  well.  The  advantage  of  using 
Gaussian  basis  sets  diminishes  as  the  number  of  types  of  quantum  particles  in  a  system 
increases. 

When  representing  molecular  orbitals  as  linear  combinations  of  Gaussian  type 
atomic  basis  functions,  optimization  of  the  orbital  exponents  becomes  ever  more  complex 
as  the  number  of  atomic  basis  functions  increases.  The  computational  resources  needed 
to  perform  these  optimizations  increases  and  the  difficulty  of  obtaining  SCF  convergence 
increases  with  increasing  basis  set  size.  Many  methods  of  optimization  result  in 
exponents  (a,)  converging  to  the  same  value  when  atomic  basis  functions  have  the  same 
azimuthal  quantum  number  (/).  This  is  a  particularly  significant  problem  with  the  muonic 
molecules  presented  in  this  paper,  due  to  the  small  separation  distance  between 
interacting  particles.  Several  methods  of  circumventing  this  problem  have  been 


60 


developed.  The  best  know  of  these  methods  are  even-tempered  [39;  40;  41;  42;  43;  44] 
and  well-tempered  [45;  46;  47;  48;  49;  50]  basis  set  methods. 

Considerable  improvement  over  even-tempered  and  well-tempered  basis  sets  can 
be  obtained  through  the  use  of  the  Legendre  polynomial  optimization  method. [51]  With 
this  method,  the  Gaussian  exponents  are  represented  as 


a,  =  exp 
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where  Ak  are  variational  parameters,  C  is  the  number  of  Gaussian  functions  combined 
and  Pk  are  orthonormal  Legendre  polynomials.  [51]  This  is  the  method  employed  to 

optimize  electron,  muon,  and  nuclear  basis  sets  for  the  results  presented  in  this  document. 

The  problem  of  coefficients  converging  to  the  same  values,  resulting  in  linear 
dependant  basis  sets,  is  more  common  when  uncontracted  Gaussian  basis  sets  are  used 
than  when  Slater  and  contracted  Gaussian  basis  sets  are  used.  This  is  primarily  a  result  of 
smaller  values  of  C  commonly  being  used  with  these  types  of  basis  sets. 

In  order  to  improve  the  basis  sets  used  to  represent  quantum  orbitals,  additional 
basis  functions  can  be  added  to  basis  sets,  or  auxilary  basis  sets  can  be  fonned  and 
centered  at  different  geometric  locations.  When  applying  limited  computational 
resources  there  is  no  hard-and-fast  improvement  rule  whether  it  is  better  to  increase  the 
number  of  auxiliary  basis  set  centers  or  to  increase  the  size  of  the  basis  sets.  The  answer 
to  this  dilemma  depends  on  the  system  being  studied  and  the  basis  sets  being  compared. 
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4.3  Applications 

The  effects  of  correlation  energy,  basis  set  center  coordinates,  and  basis  set  sizes, 
for  pppp  and  lupt  are  addressed  in  this  section.  Here,  p  and  t  represent  H  and  H  nuclei, 
respectively,  and  p  represents  a  negative-muon.  In  this  chapter,  dynamic  correlation 
energy  is  calculated  using  configuration  interactions  (Cl)  methods. 

One  of  the  difficulties  encountered  treating  molecular  systems  having  different 
types  of  quantum  particles  that  have  similar  masses,  using  GPO  methods,  is  that  basis 
sets  in  general  are  not  transferable  between  molecular  systems  or  even  different 
molecular  geometries.  In  general  basis  sets  and  system  geometries  are  strongly  coupled 
except  when  all  oppositely  charged  quantum  particle  types  have  significantly  different 
masses.  As  a  result,  the  basis  set  coefficients  and  molecular  geometry  must  be  optimized 
simultaneously  to  attain  accurate  results.  In  order  to  perform  a  full  optimization  of  the 
molecules  presented  in  this  document,  a  basis  set  for  each  type  of  particle  was  optimized 
in  turn,  then  the  geometry  was  optimized.  This  process  was  repeated  iteratively,  until  a 
minimum  energy  configuration  was  found. 

4.3.1  Contributions  to  the  Cl  Energy 

In  the  discussion  which  follows,  the  results  of  applying  the  GPO/CI  method  to 
muonic  hydrogen  molecules  is  discussed.  Calculations  of  separate  particle-particle 
correlation  energy  for  different  types  of  particle  interactions  is  presented.  The  effects  of 
basis  set  size  on  the  Cl  energy  of  these  systems  are  addressed  and  the  limitations  of 
optimizing  basis  sets  and  geometries  at  the  HF  level,  then  using  these  optimized 
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parameters  to  perform  Cl  calculations,  rather  than  optimizing  the  parameters  at  the  Cl 
level  is  evaluated. 

Full-CI  (FCI)  can  be  an  accurate  method  of  calculating  non-relativistic  correlation 
energy,  but  its  usefulness  is  limited  by  computational  requirement  that  can  be 
overwhelming  when  calculating  configuration  interactions  between  different  types  of 
particles.  The  calculation  of  Cl  energy  requires  the  solution  of  a  Slater  detenninant  ( i.e ., 
Cl  determinant).  The  majority  of  the  computational  resources  needed  to  calculate  Cl 
energy  results  from  the  storage  of  the  elements  of  Cl  matrices  and  the  solution  of  these 
detenninants.  If  the  calculation  of  the  Cl  energy  for  one  type  of  quantum  particle  ( e.g ., 
muons)  requires  finding  the  roots  of  an  (m  x  m )  Hamiltonian,  and  the  calculation  of  the 
Cl  energy  of  a  second  type  of  quantum  particle  (e.g.,  electrons)  requires  finding  the  roots 
of  a  (p  x  p)  Hamiltonian,  then  the  calculation  of  the  Cl  energy  between  the  different  types 
of  particles  (e.g.,  electrons  and  muons)  requires  finding  the  roots  of  a  Hamiltonian  that  is 
(m  •  p)x  ( m  •  p )  in  size.  As  a  result,  the  calculation  of  FCI  energy  for  systems 

containing  multiple  types  of  quantum  particles  requires  significantly  greater 
computational  resources  than  are  required  for  molecular  systems  containing  only  one 
type  of  quantum  particle. 
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Table  4-1.  Effects  of  basis  set  size  on  HF  and  FCI  energy  of  pppp.  The  FCI  energy  is 
calculated  separately  for  different  types  of  interactions.  Basis  sets  and  p-p  bond 
distances  were  optimized  at  the  HF  level.  For  columns  2  thru  4,  the  Cl  energy  was 
calculated  using  FCI  methods  and  the  basis  sets  indicated.  In  column  5  the  HF 
energy  was  calculated  using  4s3p  muon  basis  sets  and  a  2s6p  proton  basis  sets,  pp  Cl 
energy  and  pp  Cl  energy  was  calculated  at  the  FCI  level  also  using  4s3p  and  2s6p 
basis  sets  for  muons  and  protons  respectively.  The  pp  Cl  energy  was  calculated 
using  FCI  with  4s3p  muon  and  2s2p  proton  basis  sets.  All  of  the  particles  were 
treated  quantum  mechanically. 


Muon  basis  sets 

2s 

4s3p 

4s3p 

Combined 

Proton  Basis  Sets 

2s2p 

2s2p 

2s6p 

Results 

HF  energy3 

-2.955674 

-2.962424 

-2.962480 

-2.962480 

p — p  Cl  energy3 

-0.112286 

-0.193636 

-0.193667 

-0.193667 

p — p  Cl  energy3 

-0.076096 

-0.082950 

-0.083391 

-0.083391 

p — p  Cl  energy3 

-0.227296 

-0.845830 

-0.845830 

Total  Cl  energy3 

-0.415678 

-1.122416 

-1.122888 

HF/CI  energy3 

-3.371352 

-4.084839 

-4.085368 

“Energy  is  reported  in  keV 


Fortunately,  there  is  no  requirement  that  the  same  method  of  calculating  dynamic 
correlation  energy  be  used  to  calculate  the  interactions  between  all  types  of  particles  in  a 
molecular  system.  Different  methods  of  calculating  correlation  energy  ( e.g .,  perturbation 
theory,  explicit  correlation,  etc.)  can  be  used  for  different  interactions  in  the  same 
molecular  system.  In  Table  4-1,  the  correlation  energy  as  a  function  of  basis  set  size  and 
particle  type  for  pppp  is  presented.  Several  significant  factors  can  be  noted  from  an 
analysis  of  these  calculations.  1)  The  pp ,  pp,  and  pp  correlation  energy  all  contribute 
significantly  to  the  overall  energy  of  the  systems.  2)  When  comparing  the  correlation 
energy  between  particles  of  the  same  type  (e.g.,  pp  or  pp)  the  smaller  the  mass  of  the 
particles  ( i.e .,  greater  the  quantum  character)  the  greater  the  correlation  energy.  3)  The 
correlation  energy  between  particles  of  different  types  (e.g.,  pp)  is  greater  than  the 
correlation  energy  between  particles  of  the  same  type.  4)  For  the  systems  studied,  the 
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calculation  of  correlation  energy  using  GPO/CI  methods  is  affected  significantly  by  basis 
set  size,  but  not  as  significantly  as  is  the  Hartree-Fock  energy.  As  a  result,  it  is  possible 
to  improve  computational  efficiency  by  calculating  the  HF  energy  with  larger  basis  sets 
than  are  used  for  the  more  expensive  Cl  calculations  and  then  add  the  HF  and  Cl  energies 
together.  An  example  of  where  this  has  been  done  is  given  in  Column  5  of  Table  4-1. 

Studies  were  performed  which  optimized  geometry  and  basis  sets  at  the  HF  level, 
then  used  these  optimized  parameters  to  calculate  the  Cl  energy.  [2]  The  Cl  results 
presented  thus  far  in  this  paper  are  the  result  of  this  type  of  analysis.  While  it  can  be 
argued  that  this  type  of  analysis  can  be  used  to  efficiently  optimize  geometries  and  basis 
sets,  particularly  when  only  one  type  of  quantum  particle  is  present  or  when  there  is  a 
large  mass  difference  between  the  types  of  quantum  particles  studied  ( e.g nuclei  and 
electrons),  for  the  muonic  molecular  systems  presented  in  this  paper  sizable  errors 
resulted  from  this  approach  (see  Table  4-2).  While  more  accurate  results  are  produced, 
optimizing  the  geometry  and  basis  sets  at  the  FCI  level  can  present  challenges  for 
systems  of  the  type  presented  in  this  chapter.  In  order  to  obtain  accurate  results,  the 
basis  sets  and  geometries  must  be  optimized  simultaneously.  Often  hundreds,  or  even 
thousands  of  individual  GPO  calculations  must  be  carried  out  in  order  to  perform  these 
optimizations.  As  a  result,  there  are  significant  computational  limitations  on  how  large 
and  how  many  basis  sets  can  be  used  when  FCI  optimization  is  perfonned.  Table  4-2 
shows  the  results  of  two  FCI  level  calculations  performed  using  the  same  size  basis  sets, 
one  with  HF  optimized  basis  sets  and  geometry,  and  one  with  Cl  optimized  parameters. 
As  can  be  seen  from  the  table,  the  difference  in  the  calculated  HF/CI  energy  is 
noteworthy,  but  the  difference  in  the  calculated  bond  length  is  of  even  greater 
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significance.  For  muonic  molecular  systems  of  the  type  presented  in  this  paper,  using 
limited  computational  resources,  it  is  better  to  use  a  relatively  small  basis  set  and  perfonn 
a  FCI  optimization  than  to  use  a  larger  basis  set  and  optimize  at  the  HF  level. 


Table  4-2.  HF  optimized  verses  HF/CI  optimized  Cl  calculations  of 
ppp+.  The  results  in  column  2  were  obtained  by  optimizing  the  basis  sets 
and  bond  length  at  the  HF  level,  then  calculating  the  FCI  energy  of  the 
muonic  molecular  system.  The  results  in  column  3  were  obtained  by 
optimizing  the  basis  sets  and  bond  length  at  the  FCI  level.  The  muon 
and  proton  basis  sets  were  centered  at  the  same  locations,  the  positions 
of  greatest  proton  density.  4s3p  muon  basis  sets  and  a  2s2p  proton  basis 
sets  were  used  for  the  calculations.  All  of  the  particles  were  treated 
quantum  mechanically. 


HF  optimized 

FCI  optimized 

Bond  length3 

0.006088 

0.003655 

HF  energy13 

-1.4061 

-1.3959 

p — p  correlation  energy13 

-0.06006 

-0.06712 

p — p  correlation  energy13 

-0.47942 

-0.50431 

HF/CI  energy13 

-1.9404 

-1.9674 

aThe  equilibrium  proton  separation  distance  is  reported  in  angstroms  (A) 
bEnergy  is  reported  in  keV 


4.3.2  Di -muonic  Hydrogen  Molecules 

In  this  section  two  studies  which  used  the  GPO  method  are  presented:  1)  the 
geometry  and  binding  energy  of  pppp  is  optimized  using  various  basis  sets,  and  2)  the 
bond  length  of  tppt  was  calculated  and  compared  using  fixed  and  quantum  tritons. 
Molecular  geometries  ( i.e .,  equilibrium  nuclei  positions)  were  determined  from  the 
expectation  values  of  particle  density.  The  term  bond  length  refers  to  the  distance 
between  the  expectation  values  of  nuclei  density.  For  the  symmetric  muonic  molecules 
presented  in  this  paper,  the  coordinates  of  highest  nuclear  density  are  the  same  as  the 
coordinates  of  the  optimized  basis  set  positions  to  at  least  4  significant  figures  (A)p 
This  is  not,  however,  expected  to  be  the  case  for  non-symmetric  molecular  systems. 
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In  order  to  accurately  calculate  the  geometry  and  binding  energy  of  the  muonic 
molecules  presented  in  this  section  the  basis  sets  and  the  basis  set  center  locations  were 
simultaneously  optimized.  Use  of  basis  sets  optimized  at  other  center  locations,  or  the 
use  of  basis  set  center  locations  optimized  with  other  basis  sets  often  resulted  in 
differences  of  several  electron  volts  (eV)  in  the  calculated  binding  energy.  Likewise, 
changing  the  basis  set  size  of  one  particle  necessitated  the  optimization  of  the  basis  sets 
of  the  other  particles,  if  accurate  results  were  to  be  obtained.  Due  to  the  symmetry  of  the 
molecular  systems  presented  in  this  paper,  basis  sets  for  equivalent  particles  were 
constrained  to  have  equivalent  basis  sets. 

Due  to  the  computational  resources  needed  to  accurately  include  correlation 
energy  and  other  post  Hartree-Fock  corrections  to  the  calculations;  it  is  necessary  to  limit 
the  size  of  basis  sets.  Results  obtained  using  different  size  basis  sets  are  compared  in 
Tables  4-3  through  4-5.  Table  4-3  compares  the  binding  energy  for  various  size  basis 
sets.  Basis  sets  and  bond  distances  were  optimized  at  the  HF  level,  then  HF/FCI 
calculations  were  performed  under  the  optimized  conditions.  For  comparison  purposes, 
the  bond  distance  and  binding  energy  of  ppp+ ,  calculated  at  a  comparable  level  of  theory, 
is  also  shown.  Table  4-4  compares  three  FCI  calculations  which  have  similar  time  and 
memory  requirements,  but  for  which  the  number  of  Gaussian  equations  used  are  divided 
up  differently  between  muons  and  protons.  Table  4-5  compares  results  in  which  different 
numbers  of  molecular  orbitals  are  included  in  the  Cl  calculations  [/.£.,  complete  active 
space  method  (CAS)].  The  molecular  orbitals  which  result  in  the  largest  contributions  to 
the  HF  energy  are  those  which  were  included  in  the  Cl  calculations. 


67 


The  calculated  pppp  binding  energy  of  410.5  eV  can  be  compared  to  previously 
published  values  of  374.5  eV  [15]  (see  Tables  4-4  and  4-5).  It  should  be  noted  that  the 
previously  published  results  neglected  pp  correlation  interactions.  The  significance  of  pp 
correlation  interactions  will  be  discussed  in  detail  in  Chapter  5. 

The  calculated  binding  energy  of ppp+  presented  in  Table  4-3  is  260.6  eV. 
This  can  be  compared  to  published  values  of  253.2  eV.[52]  The  previously  published 
ppp+  calculations  used  larger  basis  sets  than  those  used  in  the  calculations  presented  in 
this  chapter.  The  pppp  and  ppp+  results  presented  in  Table  4-3  demonstrate  weaknesses 
of  using  small  basis  sets  and  optimizing  at  the  HF  level. 


Table  4-3.  Equilibrium  bond  length  and  binding  energy  of  pppp  and  ppp+. 
The  muon  basis  sets,  proton  basis  sets  and  equilibrium  geometry  were 
optimized  simultaneously  for  each  calculation  at  the  HF  level.  The  binding 
energy  was  calculated  at  the  FCI  level  using  HF  optimized  basis  sets  and 
bond  distance.  All  of  the  particles  were  treated  quantum  mechanically. 


pppp 

pppp 

pppp 

PPP ' 

Muon  Basis  Sets 

4s3p 

4s3p 

2s 

4s3p 

Proton  Basis  Sets 

2s6p 

2s2p 

2s2p 

2s2p 

HF  Bond  Length3 

0.003639 

0.003609 

0.005377 

0.006088 

Binding  Energyb 

402.4567 

441.8192 

260.6276 

aThe  bond  length  is  reported  in  angstroms  (A) 
bThe  binding  energy  is  reported  in  eV 
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Table  4-4.  Equilibrium  bond  length  and  binding  energy  of pppp. 
The  muon  basis  sets,  proton  basis  sets  and  equilibrium  bond 
distance  were  optimized  simultaneously  for  each  calculation  at 
the  FCI  level  using  the  variational  principle.  The  binding 
energy  was  calculated  at  the  FCI  level  using  FCI  optimized 
basis  sets  and  bond  distance.  All  of  the  particles  were  treated 
quantum  mechanically. 


Muon  Basis  Sets 

2slp 

2s 

2s2p 

Proton  Basis  Sets 

2slp 

2s2p 

2s 

Bond  Length3 

0.005012 

0.004376 

0.005988 

Binding  Energyb 

410.525 

406.968 

362.173 

aThe  bond  length  is  reported  in  angstroms  (A) 
bThe  binding  energy  is  reported  in  eV 


Table  4-5.  Equilibrium  bond  length  and  binding  energy  of  pppp.  2s2p 
muon  basis  sets  and  2s2p  proton  basis  sets  were  used.  The  muon  basis 
sets,  proton  basis  sets  and  equilibrium  bond  distance  were  optimized 
simultaneously  at  the  Cl  level  using  CAS  methods  and  varying  numbers  of 
active  molecular  orbitals  (MOs).  All  of  the  particles  were  treated 
quantum  mechanically. 


Active  p'  MOs 

10 

8 

6 

4 

Active  Nuc.  MOs 

10 

8 

6 

4 

Bond  Length3 

.005013 

.004420 

.004298 

.004046 

Total  Energyb 

-3,927.0 

-3,751.2 

-3,541.2 

-3,235.6 

Binding  Energy3 

410.5103 

aThe  bond  length  is  reported  in  angstroms  (A) 
bThe  energy  is  reported  in  eV 


From  Tables  4-3,  4-4,  and  4-5  it  can  be  seen  that  correlation  interactions  result  in 
an  increase  in  the  calculated  p-p  bond  distance  in  pppp  molecules. 

The  Bom-Oppenheimer  approximation  is  not,  of  course,  an  accurate  assumption 
when  used  to  model  interactions  between  oppositely  charged  quantum  particles  of  similar 
size.  While  understanding  this,  the  author  was  uncertain  about  how  much  the  calculated 
bond-length  would  vary  if  one  or  both  of  the  nuclei  were  fixed.  Notwithstanding  the 
error  resulting  from  application  of  the  Born-Oppenheimer  approximation,  fixing  the 
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nuclei  could  eliminate  other  errors  in  the  calculations.  It  has  been  shown  that  fixing  one 
nuclei  eliminates  translational  “contamination,”  and  fixing  two  nuclei  eliminates 
translational  and  rotational  “contamination”  in  the  calculated  results. [29;  30;  53;  54] 
Table  4-6  shows  results  for  //i/i/  where  all  four  particles  are  considered  quantum 
mechanically,  where  one  of  the  nuclei  is  fixed  and  where  both  nuclei  are  fixed. 


Table  4-6.  Equilibrium  Geometry  of  tfifit.  4s3p  muon  and  2s2p 
triton  basis  sets  were  used.  Correlation  energy  corrections  were 
included  by  CAS/CI  methods  with  five  active  molecular  orbitals 
(MOs)  per  quantum  particle  (10  muon  MO’s  and  5  MO’s  for 
each  quantum  triton).  The  muon  basis  sets,  triton  basis  sets  and 
separation  distance  were  optimized  simultaneously  for  each 
calculation  at  the  Cl  level.  Both  muons  were  treated  quantum 
mechanically.  Tritons  were  treated  classically  and  quantum 
mechanically,  as  shown. _ 


Number  of  Quantum 
Protons 

0 

1 

2 

Bond  Length3 

0.003002 

0.003754 

0.004032 

aThe  bond  length  is  reported  in  angstroms  (A) 


4.4  Conclusions 

The  non-adiabatic  ab-initio  methods  outlined  in  this  paper  allow  the  study  of 
molecular  systems  containing  any  number  of  any  type  (/.<?.,  mass  and  charge)  of  quantum 
particles  in  systems  that  may  also  contain  classical  (i.e.,  fixed)  particles.  The  size  of  the 
systems  studied,  the  number  of  quantum  particles,  and  the  types  of  quantum  particles  are 
limited  only  by  the  computational  facilities  available. 

The  ab-initio  HF/CI  methods  described  in  this  paper  have  been  applied  to  some 
muonic  molecular  systems.  It  was  found  that  basis  sets  optimized  for  these  systems 
could  not  in  general  be  accurately  transferred  to  similar  systems.  In  order  to  obtain 
accurate  results,  the  basis  sets  for  each  type  of  particle  and  the  molecular  geometry 
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needed  to  be  optimized  simultaneously.  Failure  to  do  this  resulted  in  errors  that  were 
sometimes  several  tens  of  electron  volts.  It  was  additionally  shown  that  for  muonic 
hydrogen  molecules,  of  the  type  studied  in  this  paper,  the  optimization  must  include 
correlation  energy  if  accurate  results  are  to  be  obtained.  Optimizing  at  the  HF  level,  then 
adding  Cl  to  the  final  results  is  not  sufficient.  The  basis  sets  used  in  the  studies  presented 
in  this  paper  were  optimized  using  Legendre  polynomial  optimization,  because  this 
method  of  optimization  yields  better  results  (i.e.,  lower  variational  energy)  than  does  the 
more  common  even-tempered  and  well-tempered  basis  set  methods. [51] 

For  di-muonic  hydrogen  molecules  the  pp,  pp,  and  pp  correlation  interactions  all 
contribute  significantly  to  the  overall  energy  of  the  systems.  When  comparing  the 
correlation  energy  between  particles  of  the  same  type  ( e.g .,  pp  or  pp)  the  greater  the 
quantum  character  {i.e.,  smaller  the  mass)  of  the  particles,  the  greater  the  correlation 
energy.  The  correlation  energy  between  particles  of  different  types  {e.g.,  pp)  is  greater 
than  the  correlation  energy  between  particles  of  the  same  type.  The  calculation  of 
correlation  energy  using  GPO/CI  methods  is  affected  significantly  by  basis  set  size,  but 
not  affected  as  much  as  is  the  Hartree-Fock  energy  (at  least  for  the  systems  studied)  (see 
Tables  4-1  and  4-2). 

The  feasibility  of  using  different  methods  {e.g.,  different  Cl  levels)  to  calculate 
different  contributions  to  the  dynamic  correlation  energy  was  demonstrated.  When  static 
correlation  is  negligible,  as  it  is  for  all  of  the  calculations  presented  in  this  chapter,  there 
is  no  requirement  that  the  same  basis  sets,  or  even  the  same  methods  of  calculating 
dynamic  correlation  energy  be  used  for  interactions  between  all  particles,  or  all  types  of 
particles.  The  use  of  different  basis  sets  to  calculate  the  HF  energy  and  different 
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contributions  to  the  correlation  energy  was  demonstrated.  Situations  where  this  can 
improve  HF/CI  energy  calculations  were  discussed. 
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V.  Di-Muonic  Hydrogen  Molecules 


5.1  Introduction 

This  chapter  uses  the  General  Particle  Orbital  (GPO)  HF-SCF/CI  method 
presented  in  Chapter  4  to  calculate  equilibrium  bond  lengths  and  binding  energies  of  di- 
muonic  hydrogen  molecules  and  to  study  the  effects  of  correlation  interactions  on  these 
systems.  Both  parallel  and  anti-parallel  particle  spin  states  are  addressed.  Only  ground 
muonic  state  molecules  will  be  presented  in  this  chapter.  The  di-muonic  hydrogen 
molecules  studied  in  this  chapter  do  not  have  bound  excited  muonic  states.  Excited 
rotational  and  vibrational  states  of  di-muonic  hydrogen  molecules  are  addressed  in 
Chapter  7. 

Dynamic  correlation  is  shown  to  have  large  affects  on  the  physical  properties  of 
di-muonic  hydrogen  molecules.  As  a  result,  correlation  effects  of  these  molecules  were 
studied  in  detail  and  are  presented  in  this  chapter.  Dynamic  correlation  between  quantum 
particles  results  in  the  total  energy  of  molecular  systems  being  lower  than  that  which  is 
calculated  using  Hartree-Fock  (HF)  methods.  In  most  electronic  molecules  the  electron 
correlation  results  in  bond  distances  being  shorter  than  those  calculated  using  HF 
methods,  without  dynamic  correlation.  In  this  chapter  the  exotic  molecule  pimp,  where  p 
and  p  represent  protons  and  negative  muons  respectively,  is  compared  with  the  ppp+ 
molecular  ion.  Qualitative  differences  in  correlation  interactions  between  muonic 
molecules  and  corresponding  electronic  molecules  is  illustrated.  In  di-muonic  molecules 
such  as  pppp,  muon-proton  and  muon-muon  correlation  interactions  result  in  an  increase 
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in  the  equilibrium  bond  distance  compared  with  that  calculated  using  HF  methods  alone. 
This  comparison  highlights  the  dominant  influence  of  pp  correlation  energy  in  pppp 
molecules  and  analyzes  the  major  effects  of  particle  spin  multiplicity  on  the  binding 
energy  and  equilibrium  nuclear  separation. 

5.2  Methods 

In  order  to  model  equilibrium  bond  distances  and  binding  energies  using  GPO 
methods  it  is  first  necessary  to  optimize  basis  functions  for  particle  wavefunctions  and 
the  coordinates  of  basis  set  centers.  For  most  electronic  molecular  systems  these 
optimizations  are  not  strongly  dependent  on  each  other.  It  is  possible  to  employ 
optimized  basis  sets  between  different  molecular  geometries  and  similar  molecular 
systems.[l]  Such  is  not  possible  with  pppp.  The  optimum  coefficients  for  both  proton 
and  muon  wavefunctions  depend  strongly  upon  each  other  and  upon  the  coordinates  of 
the  basis  set  centers.  For  the  GPO  calculations  presented  in  this  chapter,  basis  sets  and 
basis  set  coordinates  were  optimized  simultaneously  by  an  iterative  scheme  employing 
the  variational  principal  (see  Chapter  4).  The  proton  basis  sets  were  optimized  first;  the 
optimized  proton  basis  sets  were  then  used  to  optimize  muon  basis  sets.  These  optimized 
basis  sets  were  then  used  to  optimize  the  locations  of  the  basis  set  centers.  This  process 
was  repeated  iteratively  dozens  to  hundreds  of  times  until  the  lowest  energy  configuration 
was  found.  For  all  of  the  results  presented  in  this  chapter  the  muon  and  proton  basis  sets 
were  collocated,  therefore,  only  the  separation  distance  between  the  basis  set  centers  was 
optimized.  The  basis  set  coefficients  were  optimized  using  the  Legendre  polynomial 
optimization  method  [2]  using  a  code  written  by  Gary  Kedzoria.  This  method,  which  has 
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been  shown  to  have  advantages  over  the  more  common  even-tempered  [3;  4;  5;  6;  7;  8] 
and  well-tempered  [9;  10;  11;  12;  13]  methods  is  described  elsewhere  (see  Chapter  4). [2] 

After  the  basis  sets  and  basis  set  centers  were  optimized,  it  was  possible  to 
calculate  equilibrium  bond  distances  and  molecular  binding  energies.  The  average  proton 
separation  distances  for  pppp  and  ppp  were  detennined  from  the  expectation  values  of 
proton  density.  For  these  symmetrical  molecules,  the  distance  between  the  optimized 
basis  set  centers  and  the  expectation  values  of  the  protons  densities  were  shown  to  be 
equal  to  at  least  4  significant  figures  ( i.e 1  pA)  (see  Chapter  4). 

In  order  to  study  different  contributions  to  predictions  of  the  dynamic  correlation 
energy,  the  bond  distance,  the  basis  sets  and  basis  set  centers  were  optimized  at  different 
levels  of  theory:  1)  at  the  HF  level  with  no  correlation  correction,  2)  at  the  HF  /Cl  level, 
including  only  pp  correlation,  3)  at  the  HF  /Cl  level  including  only  pp  correlation,  and  4) 
at  the  HF/CI  level,  including  all  three  types  of  particle  correlation  interactions,  pp,  pp, 
and  pp.  The  equilibrium  proton  bond  distance  was  then  calculated  for  each  of  these 
cases.  The  calculations  were  performed  using  2s  lp  proton  and  muon  basis  sets  with  full- 
CI  (FCI)  and  with  2s2p  proton  and  muon  basis  sets  using  complete  active  space  (CAS) 
methods,  where  the  Cl  calculations  used  10  active  molecular  orbitals  for  each  of  the  types 
of  particles.  The  molecular  orbitals  chosen  were  those  which  contributed  most  strongly 
to  the  HF  energy.  The  calculations  presented  in  this  chapter  treat  all  of  the  muons  and 
nuclei  quantum  mechanically. 

In  addition  to  studying  how  different  types  of  particle  correlation  contribute  to 
nuclear  separation,  their  individual  contributions  to  the  total  correlation  energy  were  also 
studied.  For  these  studies  the  variational  principal  was  used  to  minimize  the  Cl  energy 
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by  optimizing  the  basis  set  coefficients  and  center  coordinates  with  pp,  pp,  and  up 
correlation  all  three  being  considered.  These  optimized  basis  sets  and  center  locations 
were  then  used  to  calculate  the  individual  contributions  of  pp,  pp,  and  pp  to  the  total 
correlation  energy. 

The  differences  between  the  binding  energy  and  bond  distances  for  singlet  and 
triplet  muonic  and  protonic  states  of  the  pppp  molecular  system  was  also  studied. 

5.3  Results 

Table  5-1  summarizes  a  comparison  of  pppp  and  pup  bond  distances  calculated 
by  optimizing  the  coordinates  of  basis  set  centers  while  including  different  contribution 
of  proton  and  muon  correlation  in  the  calculations.  At  the  HF  level,  the  pppp,  bond 
distance  is  shorter  than  that  of  ppp+',  however,  when  particle  correlation  is  considered,  the 
pppp  bond  length  is  shown  to  be  significantly  longer  (i.e.  37%)  than  that  of  pup  .  For  the 
pup  molecular  ion,  both  pp  and  pp  correlation  contribute  to  reduce  the  pp  bond  distance. 
For  the  pppp  molecule  pp  correlation  affects  the  bond  distance  more  strongly  than  does 
pp  correlation,  shortening  the  equilibrium  bond  distance,  whereas  pp  and  pp  correlation 
both  contribute  to  an  increased  pp  separation.  For  the  analogous  Hi  molecule,  ee  and  pp 
correlation  interactions  shorten  the  equilibrium  pp  bond  distance. 
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Table  5-1.  Nuclear  bond  lengths  of  /?/////?  and  ppp+  with  respect  to  correlation 
energy.  The  bond  distances  were  optimized  assuming  1)  no  correlation 
(HF),  2)  pp  correlation  only,  3)  pp  correlation  only,  and  4)  all  three  included: 
pp,pp,  and  pp  correlation.  Row  1  was  calculated  using  2slp  muon  and  2slp 
proton  basis  sets  and  FCI.  Row  2  was  calculated  using  2s2p  muon  and  2s2p 
proton  basis  sets.  Ten  active  muon  molecular  orbitals  and  10  active  proton 
molecular  orbitals  were  used  for  the  CAS/CI  calculations.  Row  3  was 


calculated  usin' 

g  4s3p  muon  basis  sets  and  2s2p  proton  basis  sets  and  FCI. 

Bond  Distance  (A) 

HF 

pp 

correlation 

PP 

correlation 

pp,  pp,  pp 

correlation 

1 

pppp 

(2s lp) 

0.004573 

0.004700 

0.003176 

0.005012 

2 

pppp 

(2s2p) 

0.004112 

0.004700 

0.003167 

0.005013 

3 

PUP 

0.006975 

NA 

0.004349 

0.003655 

For  the  pppp  system,  individual  contributions  to  the  total  correlation  energy  were 
calculated  separately  (see  Table  5-2).  pp  correlation  energy  is  about  2!4  times  larger  than 
pp  correlation  energy.  This  is  interesting  given  the  fact  that  pp  correlation  was  shown  to 
affect  the  equilibrium  bond  distance  more  than  does  pp  correlation.  The  pp  correlation 
energy  is  about  four  times  larger  than  the  pp  correlation  energy.  These  ratios  are 
compared  to  those  for  H2  in  which  the  author’s  calculations  show  that  ee  correlation 
energy  and  ep  correlation  energy  are  about  equal  (i.e.,  ~0.9  eV)  while  pp  correlation 
energy  in  H2  is  about  20  times  less  than  ee  or  ep  correlation  energy. 
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Table  5-2.  /?/////?  correlation  energy.  Row  1  was  calculated  using  2slp  muon  basis 
sets  and  2slp  proton  basis  sets  with  FCI.  Row  2  was  calculated  using  2s2p  muon 
and  2s2p  proton  basis  sets  with  10  active  muon  molecular  orbitals  and  10  active 
proton  molecular  orbitals  used  for  the  CAS/CI  calculations. 


Correlation  Energy  (eV) 

pp 

correlation 

PP 

correlation 

pp 

correlation 

pp,  pp,  pp 

correlation 

1 

2slp 

FCI 

-183.2721 

-77.9020 

-724.1032 

-985.2773 

2 

2s2p 

CAS/CI 

-183.3006 

-77.9064 

-724.0254 

-985.2324 

The  difference  in  calculated  binding  energy  and  proton  equilibrium  bond 
distances  for  puup  when  particles  of  the  same  type  have  paired  spin  (i.e.,  singlet)  and 
when  they  have  parallel  spin  (i.e.,  triplet)  can  be  seen  in  Table  5-3.  Note  that  for  the 
system  studied,  there  are  comparatively  small  proton-proton  spin-orbit  influences.  There 
is  no  significant  difference  in  the  binding  energy  or  equilibrium  bond  distance  for  singlet 
or  triplet  proton  states.  This  can  be  attributed  to  the  larger  mass  and  separation  distance 
of  the  protons.  The  multiplicity  of  the  muons  does,  however,  have  a  significant  impact 
on  both  the  binding  energy  and  the  equilibrium  bond  distance.  When  the  pppp  molecule 
is  in  the  single  muonic  state,  pppp  is  strongly  bound,  with  a  binding  energy  of 
approximately  410  eV.  When  the  molecule  is  in  a  triplet  muonic  state  pppp  is  very 
weakly  bound  (i.e.,  -0.074  eV).  The  pp  bond  distance  increases  by  a  factor  of  almost 
five  when  going  from  a  singlet  to  a  triplet  muonic  state. 

For  conventional  molecules  (i.e.,  molecules  containing  only  nuclei  and  electrons), 
using  the  Born-Oppenheimer  approximation,  the  dynamic  correlation  energy  is  limited  to 
about  1%  of  the  corresponding  HF  energy.[14:265]  For  pppp  dynamic  correlation  energy 
is  calculated  to  be  as  much  as  60%  of  the  total  energy  (see  Table  5-3). 
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Table  5-3.  The  effect  of  particle  spin  states  on  the  binding  energy  and  equilibrium 
pp  bond  length  of  /?/////?.  FCI  calculations  were  performed  using  2slp  basis  sets. 
The  calculations  utilizing  2s2p  basis  sets  used  10  active  muon  molecular  orbitals  and 
10  active  proton  molecular  orbitals  for  the  Cl  calculations. 


p 

P 

eV 

A 

2slp 

2slp 

Ehf 

Eci 

Binding 

Energy 

Bond 

Distance 

singlet 

triplet 

-2941.7077 

-3926.9850 

410.5254 

0.005012 

triplet 

singlet 

-1389.2584 

-3543.7386 

0.074040 

0.023955 

triplet 

triplet 

-2783.2935 

-3543.7386 

0.074040 

0.023955 

p 

P 

eV 

A 

2s2p 

2s2p 

Ehf 

Eci 

Binding 

Energy 

Bond 

Distance 

singlet 

triplet 

-2941.7289 

-3926.9613 

410.5103 

0.005013 

triplet 

singlet 

-1389.2719 

-3543.7387 

0.074094 

0.023953 

triplet 

triplet 

-2783.2778 

-3543.7387 

0.074091 

0.023953 

5.4  Conclusions 

The  important  roles  that  pp,  pp  and  pp  correlation  interactions  play  in  the  binding 
of  pppp  molecules  was  demonstrated.  These  interactions  account  for  as  much  as  60%  of 
the  total  energy  of  the  muonic  system.  Of  the  three  correlation  contributions,  pp 
correlation  is  by  far  the  most  important.  This  indicates  that  the  results  of  previously 
published  calculations,  which  neglect  pp  correlation  interactions,  cannot  be  expected  to 
yield  accurate  binding  energies.[15]  As  with  other  types  of  molecules,  correlation  effects 
can  strengthen  molecular  bonding.  Unlike  most  electronic  molecules,  however, 
correlation  effects  in  the  di-muonic  hydrogen  molecule  pppp  result  in  an  increase  in  the 
nuclear  bond  distance  compared  with  ppp+ . 

Whether  nuclear  spin  is  paired  or  unpaired  does  not  have  a  significant  impact  on 
the  binding  energy  or  nuclear  bond  distance  of  pppp.  Muon  spin  does,  however 
significantly  affect  both  the  binding  energy  and  pp  bond  distance,  pppp  molecules  in  a 
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singlet  muon  state  have  a  binding  energy  of  approximately  410  eV.  When  the  muons  are 
in  a  triplet  spin  state  the  binding  energy  is  much  weaker  (i.e.,  -.074  eV). 
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VI.  Di-Muonic  Helium  Molecules 


6.1  Introduction 

Several  factors  limit  the  maximum  obtainable  muon-catalyzed  fusion  yield. 
Preeminent  amongst  these  factors  in  the  mono-muonic  mechanisms  is  the  sticking  of 
muons  to  helium  nuclei  formed  during  the  fusion  processes  outlined  in  Chapter  2  and 
Appendix  A.  This  chapter  examines  the  possibility  of  fusing  di-muonic  He  molecules, 
fonned  during  p-d  and  d-d  fusion,  liberating  muons  which  are  stuck  to  these  nuclei. 

Regardless  of  what  modifications  to  the  muon-catalyzed  fusion  reaction  cycle  are 
made,  mono-muon  catalyzed  fusion  cannot  be  a  practical  source  of  energy  production 
unless  muons  can  be  prevented  from  sticking  to  helium,  or  can  be  stripped  from  helium 
after  sticking  occurs.  Many  attempts  have  been  made  to  find  a  method  of  manipulating 
the  sticking  probability  and  of  freeing  stuck  muons.  As  of  the  writing  of  this  document, 
no  effective  methods  have  been  developed.  It  has  been  shown  that  muon  sticking  and 
stripping  are  affected  by  temperature  and  pressure.  Research  indicates,  however,  that 
adjustment  of  these  factors  cannot,  by  themselves,  solve  the  muon-catalyzed  fusion 
sticking  problem.[l] 

In  addition  to  sticking  to  helium  during  the  fusion  process,  muons  can  be 
“scavenged”  by  helium  in  a  reaction  chamber.  This  can  occur  either  as  a  result  of 
collisions  between  helium  and  free  muons  or  as  the  result  of  muon  exchange  reactions 
between  helium  and  muonic  hydrogen  molecules.  While  it  is  possible  to  address  this 
problem  by  removing  helium  to  maintain  a  low  concentration  of  helium  in  the  reaction 
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chamber;  being  that  helium  is  one  of  the  products  formed  during  muon-catalyzed  fusion, 
its  concentration  and  scavenging  need  to  be  addressed. 

Two  possible  mechanisms  for  releasing  some  of  the  muons  stuck  to  He  during 
fusion  are  the  fusion  of  two  He  nuclei  and  the  fusion  of  a  He  nuclei  and  a  deuteron. 
These  reactions  are  represented  in  Equations  6.1  and  6.2.  Helium-3,  which  is  a  product 
of  p-d  and  d-d  fusion  can  go  on  to  fuse  via  the  following  nuclear  reactions: 

3  He + 3  He - >4He  +  2p  (6.1) 

and 

3 He + d - » 4 He  +  p  (6.2) 

where  p,  and  d  represent  protons  and  deuterons  respectively.  Some  of  the  muons  that 
participate  in  these  reactions  will  be  freed,  enabling  them  to  catalyze  further  fusion 
events.  It  is  possible  that  helium  nuclei,  which  are  stuck  to  muons,  can  fuse  and  free  the 
muons  to  catalyze  additional  fusion  reactions.  Due  to  the  magnitude  of  the  Coulombic 
repulsion  of  helium  nuclei  with  other  nuclei,  single  muons  cannot  be  expected  to  catalyze 
helium  fusion  reactions.  While  most  of  the  muons  which  participate  in  reactions  6.1  and 
6.2  are  expected  to  be  liberated  by  the  reaction,  some  of  the  muons  will  stick  to  the  4 He 
nuclei  formed  during  the  reaction. [2;  3:72;  4] 


90 


If  formed,  He^iHe^i  and/or  'He/j/ud  molecules  have  the  possibility  of  fusing 
and  liberating  previously  bound  muons  in  the  process.  This  chapter  examines  the  bond 
lengths  and  vibrational  properties  of  di-muonic  He  -  He  and  ' He  -  d  molecules  to 
determine  whether  bound  forms  of  these  molecules  exist,  and  if  so,  their  equilibrium 
bond  lengths  and  vibrational  properties.  The  fusion  rate  of  these  molecules  will  depend, 
in  part,  on  the  separation  distance  between  the  nuclei. 

6.2  Methods 

In  order  to  determine  the  equilibrium  bond  lengths  of  di-muonic  ^  He  -  3 He  and 
3 He  -  d  molecules  the  General  Particle  Orbital  (GPO)  method  of  non-adiabatic  quantum 
mechanics  described  in  Chapter  4  was  used.  Dynamic  correlation  was  included  in  the 
calculations  through  the  use  of  configuration  interaction  (Cl)  methods. [5;  6]  The 
equilibrium  bond  length  was  calculated  from  the  expectation  value  of  nuclei  density.  For 
the  He  -  He  calculations  the  Born-Oppenheimer  approximation  was  used  and  the 
nuclei  positions  were  fixed  for  individual  HF/CI  calculations.  The  bond  length  for 
He  -  d  was  determined  by  fixing  the  3 He  nuclei  and  modeling  the  deuteron  and  muons 
as  quantum  particles.  The  results  of  these  calculations  were  compared  to  results  in  which 
both  of  the  nuclei  were  fixed.  The  basis  set  parameters  and  basis  set  center  locations 
were  optimized  simultaneously,  at  the  Cl  level,  using  five  active  molecular  orbitals  per 
quantum  particle. 

When  solving  the  Schrodinger  equation  for  ordinary  molecules  it  is  most  often 
assumed  that  the  volume  of  nuclei  is  negligible.  In  the  case  of  an  electron  and  a  proton 
the  error  generated  by  this  approximation  is  about  6x10  9  eV.  Because  a  muonic 
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molecular  hydrogen  radius  is  about  200  times  smaller  than  that  of  an  electronic  molecular 
hydrogen  radius,  significant  errors  can  result  from  assuming  the  nuclear  volume  of 
muonic  molecules  to  be  zero.  A  first  order  approximation  to  the  energy  correction  can  be 
calculated  using  perturbation  theory.  This  method  of  correcting  for  nuclear  volume 
effects  is  considered  accurate  as  long  as  the  perturbation  is  limited  to  a  small  fraction  of 
the  total  energy,  as  it  is  for  nuclei  from  the  first  few  rows  of  the  periodic  table.  The 
potential  energy  V(r)  is  assumed  to  result  from  a  point  charge  [5;  7:49-50;  8: 1 141-1 147] 


V(r)  =  - 


Ze2 


4n  e0  r 


(6.3) 


when  the  negative  particle  radius  (r)  is  greater  than  the  mean  nuclear  radius  (R).  When 
r  <  R  the  nuclei  can  be  considered  to  be  a  uniformly  charged  sphere  which  results  in  a 
potential  energy 


V(r)  = 


Ze 2  3 

An  e0  R  2 


1  (  r  ^ 


v*/ 


(6.4) 


where  Z  is  the  charge  of  the  nucleus  of  interest.[7:49-51;  8:1141-1147]  If  (f)  is  the 

expectation  value  of  the  potential  energy,  assuming  a  point  mass,  and  (V')  the  expectation 

value  calculated  using  equations  (6.3)  and  (6.4)  as  described,  then  the  potential  energy 
correction  ( AF)  for  nuclear  volume  is 
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AV  =  (V)-(V) 


(6.5) 


The  radial  wave  function  used  to  calculate  the  expectation  values  can  be 
approximated  by  Coulomb  hydrogenic  radial  wave  functions: 


KA)  = 
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where 


a0  = 


4 n  e0  n 
/ue1 


(6.7) 


and  I}'*  represents  associated  Laguerre  Polynomials,  n  and  /  represent  principle  and 
azimuthal  quantum  numbers  respectively.  For  the  ground  state  (Is)  this  yields: 
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and 
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which  reduces  to 


AV^  = 


8;r  e0  a0R  Z 


3 al  -3 a0R  Z2  +2 RiZi  -3a0(a0  +  RZ)~  exp 


r  2RZ  4 


lo  ) 


(6.10) 


Table  6-1  uses  the  methods  described  in  this  chapter  to  compare  the  potential 
energy  of  muonic  hydrogen  and  helium  atoms  with  and  without  nuclear  volume  included 
in  the  calculations. 


Table  6-1.  Potential  energy,  as  compared  to  free 
particles,  of  hydrogen  and  helium  muonic-atoms. _ 


Atom 

V  with  no  nuclear 
volume3 

V  with  nuclear 
volume3 

PH 

-2528.58 

-2528.54 

dp 

-2663.29 

-2663.22 

tp 

-2711.34 

-2711.24 

s  Hep 

-10845.42 

-10843.84 

4  Hep 

-10943.18 

-10941.22 

aEnergy  is  reported  in  eV 


A  potential  energy  correction  can  be  calculated  for  molecules  containing  more 


than  one  nucleus  by  first  calculating  the  atomic  correction  terms  At)  for  each  nucleus  in 


the  molecule.  The  total  molecular  correction  (AVm)  is 


N 

AV  =  V  — At) 

m  i 


i=l 


(6.11) 
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where  C,  represents  the  average  muon  charge  density  in  each  nuclei  i  ?  and  D  is  the 

muon  charge  density  in  the  nuclei  i  when  the  nuclei  and  muon  are  isolated  from  all  other 
nuclei.  Throughout  this  document,  the  nuclear  radii  ( R )  have  been  approximated  as 

R  =  1.25  femtometers. [7:48, 122] 

The  vibrational  energy  levels,  vibrational  frequency  and  magnitude  of  the  nuclear 
vibrations  were  determined  using  quasi-classical  methods.  The  quantum  vibrational 
energy  levels  between  nuclei  were  calculated  using  a  quantum  mechanical  energy  grid. 
The  quantum  vibrational  energy  levels  between  two  particles  were  calculated,  starting 
with  [9]: 


H  =  f  +  V  (6.12) 

/V  /V  A 

where  H  is  the  Hamiltonian  operator  of  the  system,  T  is  a  kinetic  energy  operator  and  V 
is  a  potential  energy  operator.  Being  that  the  kinetic  energy  operator  for  a  hannonic 
approximation  to  the  system  is  the  same  as  the  true  kinetic  energy  operator,  the 
Hamiltonian  operator  for  a  harmonic  approximation  (Hho)  to  the  system  can  be  written 
as: 


H HO  -  T  +  Vharm 


(6.13) 
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where  Vharm  is  a  harmonic  approximation  to  the  potential  energy  of  the  system.  Adding 
both  sides  of  the  above  equation  to  H  yields 

H  =  HHO  +  V-Vharm  (6.14) 

Since  it  is  possible  to  determine  analytical  expressions  for  the  expectation  values  of  Hm 

and  Vharm ,  and  the  expectation  value  of  V  can  be  determined  numerically  by  integrating  a 

fit  of  a  quantum  mechanically  determined  energy  grid,  the  expectation  value  of  H  can  be 
detennined. 

Through  the  use  of  ladder  operators  it  is  possible  to  detennine  an  n  x  n  matrix 
which  approximates  H 

^1/fjl^ 

(2  T/jl^  (2^|2)  •  •  •  (2#|«) 

H=  •  •  (6.15) 

(nH^  ■  ■  •  (nH\nj 

The  matrix  elements  of  H  are  defined  as: 
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where 


a  = 


v  n 


(6.17) 


and 


co  = 


(6.18) 


V"(Req  )  is  the  second  derivative  of  the  potential  energy  with  respect  to  bond  length, 
evaluated  at  the  equilibrium  geometry  R  .  De  is  the  equilibrium  dissociation  energy;  fi 
is  the  reduced  Planck’s  constant;  and  jxred  is  the  reduced  mass  between  the  particles  that 
were  fixed  for  each  quantum  mechanical  calculation. 

If  y/;(z)  is  defined  as  the  harmonic-oscillator  wave  function  (z  |  i\ 
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(6.19) 


!//,(z)  =  (z | i)  =  NiHi(az)e  2 


where  H^az)  represents  Hennite  polynomials: 


H0(az)  =  l 
Hx(az )  =  2  az 

H2(az)  =  4a2z2 -2  (6.20) 

Hj  ( az )  =  2ctrz  -  2n  H j  2  (az) 


and  Nj  refers  to  nonnalization  constants: 


N..  = 


a 


7T22‘i\ 


(6.21) 


z  is  the  nuclei  separation  distance  (x)  minus  the  equilibrium  bond  length  (Rcq) 


z 


x  R 


eq 


(6.22) 


If  eigenvalues  of  H  are  less  than  the  molecules  binding  energy,  then  they  represent 
quantum  vibrational  energy  levels  of  the  molecule.  The  larger  n  in  Equation  6.15,  the 
greater  the  accuracy  of  a  given  eigenvalue  i. 
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By  taking  a  numerical  first  derivative  of  the  potential  energy  with  respect  to  bond 
length  and  setting  the  derivative  equal  to  zero,  the  potential  energy  minimum  ( i.e ., 
equilibrium  geometry)  can  be  found.  The  hannonic  vibrational  frequency  of  the  nuclei  is 


i  f" K) 

V  Pred 


(6.23) 


In  order  to  observe  the  anhannonic  effects  of  internal  kinetic  energy  on  the 
vibrational  frequency,  Hamilton’s  equation  of  motion 


GH 

dq 


=  ~P 


(6.24) 


was  solved  for  position  ( q )  and  momentum  (p)  (i.e.,  in  this  case  q  is  equal  to  the 

internuclear  separation  (jc)).[  10:992- 1023]  The  symbol  H  represents  a  classical 
Hamiltonion  and  is  equal  to  the  sum  of  the  potential  and  kinetic  energy  operators.  The 
vibrational  spectra  (i.e.,  power  spectral  density)  of  the  system  was  generated  by  taking 
the  Fourier  transform  of  the  position,  momentum,  kinetic  energy  and  potential  energy  as 
functions  of  time.[l  1:600-717;  12;  13] 
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6.3  Results  and  Discussion 


6.3.1  u  3 He  +ju  3//e 

Due  to  its  small  size,  (  He^j.)  interacts  with  charged  particles  in  a  manner  similar 
to  how  hydrogen  nuclei  interact,  except  when  molecular  separation  distances  are  small  ( 
<0.01  A).  Due  to  the  larger  muon  mass  relative  to  an  electron,  muonic  helium  orbitals 
are  approximately  l/200th  the  size  of  a  helium  electronic  orbital.  A  significant  portion  of 
the  muon  density  is  calculated  to  be  located  within  the  helium  nucleus  ( i.e .,  ~3%).  As  a 
result  of  its  small  size,  two  (  He/j.)  ions  interact  with  each  other  similarly  to  two  +1  point 
charges,  except  when  the  nuclei  are  very  close  to  each  other  (  <0.01  A).  When  the 
molecular  separation  is  small,  both  muons  interact  with  both  nuclei  resulting  in  a 
localized  potential  energy  minimum.  These  results  can  be  seen  in  Figures  6- la,  6- lb  and 
6-lc. 
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(a) 


(b) 


0.2  0.4  0.6 

Bond  Length  (A) 


(c) 


Figure  6-la,  6-lb  and  6-lc.  ( 3  He  jj,3  He  /uj  potential  energy  verses  bond  length  with 

respect  to  free  (3He^)+  ions.  4s3p  muon  basis  sets  were  centered  on  classical  (i.e.,  fixed) 3 He 
nuclei.  The  localized  energy  minimum  occurs  at  0.003393  A. 


The  eigenvalues  of  H  were  calculated  for  the  localized  energy  minimum.  There 
are  no  bound  vibrational  states  (see  Figure  6-2).  From  these  results  it  can  be  concluded 

that  the  di-muonic  helium  ion  (3 He /u3 He //)  +  is  unstable  and  cannot  be  a  pathway  for 


producing  significant  quantities  of  fusion  energy. 


It  is  likely  that  (  and 
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(  '  He  p'  Hep)  2e  form  bound  molecules,  however,  the  He-He  bond  lengths  for  these 

molecules  are  expected  to  be  approximately  equal  to  the  bond  lengths  of  '//2‘  and  H 2 

respectively.  Therefore  these  molecules  are  expected  to  have  negligible  fusion  rates. 
Including  nuclear  volume  effects  in  the  calculations  does  not  change  these  results.  Rather 

it  results  in  a  calculated  eigenvalues  of  H  v=0  and  v=l,  which  are  slightly  more 
negative  than  when  nuclear  volume  is  neglected  in  the  calculations. 


*  /  3  3  \2+ 

Figure  6-2.  Eigenvalues  v=0  and  v=l  of  H  for  y  Hep  Hep)  .  This 

shows  there  are  no  bound  vibrational  energy  levels  of  the  localized 
minimum. 

6.3.2  3He  +  d 

(3 Heppd )  was  modeled  two  ways:  1)  using  classical  nuclei  and  quantum  muons 

and  2)  using  fixed  He  coordinates  (i.e.,  classical  He),  a  quantum  deuteron  and  quantum 
muons.  Both  methods  yielded  similar  results.  All  of  the  calculations  were  performed  at 
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the  FCI  level  with  4s3p  muon  basis  sets.  A  2s2p  deuteron  basis  set  was  used.  The  muon 
basis  sets  were  centered  on  the  classical  nuclei  and  collocated  with  the  quantum  deuteron 
basis  set  center  coordinates.  The  basis  sets  and  basis  set  center  coordinates  were 
optimized  simultaneously. 

A  potential  energy  curve  was  calculated  using  the  Bom-Oppenheimer 
approximation.  The  energy  was  corrected  to  account  for  the  muon  density  which  is 
located  within  the  nuclei.  The  results  of  this  correction  can  be  seen  in  Figure  6-3.  The 
3He-d  equilibrium  bond  length  was  calculated  to  be  0.003732  A.  There  are  two  bound 
vibrational  states,  the  ground  state  (v=0)  and  the  first  excited  vibrational  state  (v=l)  (see 
Figure  6-4). 


—  —  Uncorrected 


Nuclear  Volume 
Corrected 


Figure  6-3.  Comparison  of  potential  energy  curves  with  and  without 

nuclear  volume  being  considered.  The  energy  represents  the  total  energy  of  the 
muonic  ion  with  respect  to  free  particles.  The  nuclei  were  considered  classically  in 
the  calculations  shown. 
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Figure  6-4.  ^Hejujud^ potential  energy  curve  with  nuclear  volume  and 

bound  vibrational  energy  levels  included  in  the  calculations.  The  ground 
(v=0)  and  first  (v=l)  quantum  vibrational  energy  states  are  shown.  The 
potential  energy  minimum  occurs  at  0.003732  A.  The  quantum 
mechanical  points  were  interpolated  between  using  cubic  spline 
interpolation  with  fixed  endpoints.  The  nuclei  were  considered  classically 
in  the  calculations  shown.  A  plot  of  the  potential  energy  surface  and 
vibrational  energy  levels  calculated  without  including  nuclear  volume 
effects  looks  almost  identical  on  this  scale. 


Inclusion  of  nuclear  volume  effects  in  the  calculations  results  in  a  decrease  in  the 
calculated  binding  energy  (see  Table  6-1  and  6-3).  This  is  particularly  the  case  for  the 
v=l  state,  in  which  inclusion  of  nuclear  volume  effects  results  in  a  calculated  binding 
energy  that  is  almost  an  order  of  magnitude  less  than  when  nuclear  volume  effects  are 
neglected. 
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Table  6-2.  Binding  energy  of  the  v=0  and  v=l 
vibrational  states  of  (}Hejujudy  . 


Binding  Energy  (eV) 

Nuclear  Volume 

No  Volume 

Corrected 

Correction 

v=0 

364.9872 

370.6573 

V=1 

0.1633 

1.4430 

Table  6-2  shows  the  minimum  and  maximum  bond  lengths  which  occur  during 
nuclear  vibrations.  As  can  be  seen,  inclusion  of  nuclear  volume  in  the  calculations  has 
very  little  effect  on  the  calculated  minimum  separation  distance  between  molecules.  The 
maximum  bond  length  of  the  v=l  state  is  not  reported  since  very  small  errors  in  the 
calculated  binding  energy  result  in  large  errors  in  the  maximum  bond  length  (see  Figure 
6-4). 


Table  6-3.  The  magnitude  of  nuclear  vibrations  of  calculated  with 

classical  nuclei,  with  and  without  nuclear  volume  corrections. _ 


Nuclear  Volume  Corrected 

No  Volume  Correction 

Minimum 
Bond  Length 

Maximum 
Bond  Length 

Minimum 
Bond  Length 

Maximum 
Bond  Length 

v=0 

0.00255026 

0.00715116 

0.00255702 

0.00714663 

v=l 

0.00228465 

0.00228737 

The  vibrational  spectra  and  fundamental  vibrational  frequency  of  (3 He /i fid ) 

nuclei  in  the  ground  vibrational  state  (v=0)  was  calculated  by  propagating  a  classical 
trajectory  on  the  quasi-classical  potential  energy  curve  described  previously  (see  Figure 
6-4).  The  fundamental  vibrational  frequency  was  calculated  to  be  1.515  x  10  Hz.  This 

1 7 

is  considerably  smaller  than  the  hannonic  vibrational  frequency  of  2.408  x  10  Hz 
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calculated  from  the  Hessian  at  the  potential  energy  minimum.  The  large  amount  of 
anhaminisity  which  occurs  in  the  bound  vibrational  states  of  this  ion  is  the  result  of  large 
changes  in  energy  which  occur  with  small  shifts  in  nuclei  position  of  closely  bound 
muonic  molecules.  The  excited  state  vibrational  frequency  contains  a  great  degree  of 
anharmonicity  and  is  highly  dependent  on  the  energy  of  the  v=l  state.  Small  errors  in  the 
calculation  of  the  potential  energy  result  in  large  errors  in  the  calculated  frequency.  As  a 
result,  the  methods  presented  in  this  paper  cannot  be  used  to  determine  the  vibrational 
frequency  of  the  v=l  state. 

10 
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Figure  6-5.  Calculated  ground  state  (v=0) 

vibrational  spectra  of  nuclei. 

j 

Using  a  fixed  He  nucleus  and  a  quantum  deuteron,  the  bond  length  (/)  was 
calculated  as  the  distance  between  the  3 He  coordinate  and  the  expectation  value  of  the 
deuteron  density.  The  deuteron  basis  set  was  transferred  from  an  optimized  d/u/ud 

j 

molecule.  The  muon  basis  sets  were  centered  on  the  He  and  on  the  coordinates  of  the 
expectation  value  of  the  deuteron  density.  The  muon  basis  sets  were  optimized 
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simultaneously  for  the  (3 He/i/id, )  molecular  ion.  The  results  yielded  an  equilibrium 

bond  length  r  =  .002375  A  which  is  shorter  than  that  calculated  with  two  fixed  nuclei  r  = 
0.003732  A.  These  bond  lengths  are  shorter  than  the  bond  lengths  of  muonic  hydrogen 
molecules  which  have  been  shown  to  have  large  fusion  rate  constants  (>  107  s'1). [14] 

6.4.  Conclusions 

^3Heju3He/u3j  does  not  fonn  a  bound  molecule.  A  quasi-classical  potential 

energy  surface  indicates  that  a  localized  potential  energy  minimum  exists  when  the  nuclei 
are  separated  by  about  3.4  mA;  however,  no  bound  vibrational  states  occur  in  this  region. 
These  results  indicate  that  3  He/u3  He/u  will  not  participate  in  a  feasible  muon-catalyzed 
fusion  reaction  path,  nor  will  it  result  in  a  viable  muon  stripping  mechanism. 

(  'Hepipid )  forms  a  bound  system  with  two  bound  vibrational  states  (v=0  and 

v=l).  The  binding  energy  of  these  states  is  fi(v=0)  =  365  eV  and  E(v=  1)  =  0.16  eV.  The 
calculation  of  the  excited  state  binding  energy  is  strongly  influenced  by  nuclei  volume. 
Neglecting  nuclei  volume  results  in  almost  an  order  of  magnitude  greater  binding  energy 
being  calculated.  The  bond  length  between  the  nuclei  was  calculated  using  the  Bom- 
Oppenheimer  approximation  to  be  3.7  mA  and  was  calculated  to  be  2.4  mA  when  the 
deuteron  was  considered  non-adiabatically.  The  ground  state  vibrational  frequency  is 
152  PHz. 
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In  conclusion,  the  fonnation  of  ( "  He pipid )  has  the  potential  of  enhancing  the 
muon-catalyzed  reaction  rate  (i.e.,  number  of  fusions  catalyzed  per  muon). 
(  '  He ju  ' He ju^  cannot  significantly  change  the  overall  fusion  rate. 
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VII.  Reactions  paths  of p/i/ip(T>  and  ppppp<T> 


7.1  Introduction 

If  formed  in  a  muon-catalyzed  fusion  chamber,  di-muonic  hydrogen  molecules 
are  expected  to  affect  the  muon-catalyzed  fusion  rate.  Whether  the  effect  will  be  positive 
or  negative  depends  on  the  rate  the  di-muonic  molecules  fuse  or  transfonn  into  molecules 
which  fuse.  If  this  process  is  fast,  compared  to  single  fusion  reaction  paths,  the  fusion 
yield  will  be  enhanced.  If  the  process  is  slow,  di-muonic  reactions  will  have  a  quenching 
effect  on  the  muon-catalyzed  fusion  process. 

Di-muonic  hydrogen  molecules  can  form  in  a  triplet  muon  spin  state.  Due  to  the 
large  bond  length  of  these  molecules,  they  cannot  be  expected  to  have  a  significant  fusion 
rate  (see  Chapter  5).  As  a  result,  if  these  molecules  are  to  enhance  the  muon-catalyzed 
fusion  rate,  reaction  paths  that  transform  these  molecules  into  species  which  do  fuse 
rapidly  must  exist.  The  probability  of  these  reactions  occurring  depends  on  the  energy  of 
the  excited  vibrational  and  rotational  states  and  upon  the  methods  available  to  transfer 
this  energy.  In  this  chapter  a  calculation  of  the  energy  of  these  modes  will  be  presented 
and  methods  of  transferring  this  energy  will  be  examined. 

While  there  are  many  possible  isotopic  combinations  of  muonic  hydrogen 
molecules  which  can  be  fonned,  molecules  containing  only  protons,  muons,  and  electrons 
will  be  used  as  proxies  in  this  chapter  for  molecules  containing  other  isotopic  hydrogen 
combinations.  Deuterons  and  tritons  are  expected  to  undergo  similar  muonic  reactions  to 
those  experienced  by  protons,  with  the  exception  that  some  of  the  muonic  molecules 
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formed  with  the  heavier  nuclei  may  have  appreciable  fusion  rates.  Due  to  the  lower  mass 
and  increased  quantum  character  of  protons,  as  compared  to  deuterons  and  tritons,  the 
binding  energy  will  be  lower,  and  the  bond  distances  greater  than  for  molecules 
composed  exclusively  of  the  larger  nuclei. 

7.2  Methods 

In  order  for  a  reaction  to  occur  it  must  be  energetically  favorable  and  a  viable 
reaction  path  must  exist  which  results  in  an  overall  positive  change  in  entropy.  The 
energy  of  a  reaction  (A E)  can  be  defined  as  the  difference  in  the  total  energy  between 
reactants  and  products.  Activation  energy  (Eact)  is  the  energy  of  any  barrier  between 
reactants  and  products  that  must  be  sunnounted  or  tunneled  through  in  order  to  get  a 
reaction  to  proceed.  If  A E  is  negative  the  reaction  is  exothennic  and  will  proceed,  with 
its  rate  being  determined,  in  part,  by  Eact.  When  A E  is  positive  the  reaction  is 
endothennic  and  will  not  proceed  unless  energy  is  added  to  the  system. 

Basis  set  size,  methods  of  including  correlation  energy,  and  many  other  factors 
contribute  to  the  accuracy  to  which  the  energy  of  a  system  can  be  detennined.  Since  the 
same  types  of  computational  errors  typically  occur  amongst  reactants  and  products,  if  the 
same  computational  methods  and  levels  of  theory  are  used  to  calculate  reactants, 
intennediates  ( e.g transition  states)  and  products,  many  of  the  same  errors  will  occur  on 
both  sides  of  the  equation  and  will  therefore  cancel.  It  is  therefore  important  when 
calculating  A E  and  Eact  to  use  the  same  basis  sets,  methods,  and  levels  of  including 
correlation  energy  for  all  of  the  species  participating  in  the  reactions. 
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The  choice  of  basis  sets  significantly  affects  the  calculated  values  of  A E.  As  was 
mentioned  in  previous  chapters,  transferring  basis  sets  optimized  for  one  muonic 
molecule  to  another  molecular  system,  instead  of  re-optimizing  the  basis  sets  can  result  in 
significantly  higher  energies  being  calculated.  The  choice  of  basis  sets  optimized  can  be 
used  to  determine  error  limits  on  the  calculated  values  of  A E.  If  basis  sets  optimized  for 
reactants  are  used,  the  calculated  values  of  AE  and  Eact  are  upper  bounds  to  the  most 
accurate  values  obtainable  at  a  given  level  of  theory.  If  basis  sets  are  optimized  using 
products,  the  values  of  AE  and/or  Eact  calculated  will  be  lower  bounds.  Therefore,  while 
the  optimum  basis  set  for  a  given  reaction  may  not  be  achieved,  basis  set  optimization 
error  can  be  bounded.  Basis  sets  for  reactants  and  products  were  optimized  using  the 
procedures  presented  in  Chapter  4.  In  order  to  establish  error  limits,  AE  was  calculated 
using  both  reactant  and  product  optimized  basis  sets.  Bond  lengths  were  optimized 
individually  for  each  molecule  being  studied. 

The  calculations  presented  in  this  chapter  were  performed  using  2s lp  muon  and 
proton  basis  sets.  Correlation  interactions  were  included  in  most  of  the  di-nuclear 
calculations  using  full  Cl  (FCI)  methods.  Using  available  software  and  computational 
facilities,  it  was  not  possible  to  fully  optimize  the  ppppp+  basis  sets  at  the  FCI  level 
For  this  reason,  ppppp+  calculations  were  perfonned  using  basis  sets  optimized  for 
ptu p 1 .  ppppp  *  bond  length  calculations  were  perfonned  using  8  active  muon 
molecular  orbitals  and  8  active  proton  molecular  orbitals.  It  is  likely  that  future  software 
modifications  will  parallelize  and  take  advantage  of  symmetry  in  the  Cl  code,  thereby 
allowing  larger  basis  sets  to  be  used  and  larger  molecules  to  be  studied. 
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The  quantum  vibrational  energy  levels  of  the  molecules  presented  in  this  chapter 
were  calculated  using  the  methods  presented  in  Chapter  6.  The  rotational  energy  (c)  was 
determined  by  scaling  results  obtained  for  similar  molecules.  The  scaling  factors  were 
determined  by  assuming  the  molecules  to  be  rigid  rotors.  The  calculated  rotational 
energy  of  ppp+  was  scaled  to  determine  the  rotational  energy  levels  of  the  ppp p 
molecules,  using  the  relationship, 


^2  2  ) 


SiW  K 
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(7.1) 


which  was  derived  from  the  rotational  energy  eigenvalue  equation  for  a  rigid  rotor. 
Molecule  (1)  in  equation  7.1  is  ppp+  and  molecule  (2)  is  pppp .  /represents  the 
rotational  quantum  numbers  ( i.e J  =  0,  1,  2,  .  .  .),  m,  is  the  mass  of  the  molecules,  and  rl 
is  the  distance  between  the  expectation  values  of  the  nuclei  (/  =  1,  2).  For  the  muonic 
molecules,  the  mass  was  accounted  for  in  the  scaling  factor  by  collapsing  the  proton  mass 
onto  the  expectation  coordinates  of  proton  density  and  superimposing  the  muon  mass 
distribution  to  get  the  total  mass  distribution.  The  rotational  energy  of  the  oblate 
symmetric  top  ppppp+  was  similarly  detennined  by  scaling  H [  .  Being  that  all  of  the 

rotational  energy  levels  of  interest  are  known  for  lH^ ,  a  different  scaling  factor  was 
determined  for  each  rotational  energy  level,  using  the  same  rotational  quantum  levels.[l; 
2;  3;  4;  5] 
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(7.2) 


f  \ 

s2{J,K)  =  sl{J,K)d2  r^~ 

\m2  ) 

r 

Where  d  is  the  ratio  of  the  bond  distances  — .  J  and  K  represent  the  two  rotational 

r2 

quantum  numbers  associated  with  symmetric  top  molecules.  J  corresponds  to  the 
principal  axis  of  angular  momentum  and  K  corresponds  to  the  angular  momentum  along 
the  top  axis  ( i.e .,  the  axes  which  has  a  unique  moment  of  inertia). 

7.3  Results  and  Discussion 

Pairs  of  indistinguishable  Fermions  can  exist  in  both  singlet  {i.e.,  anti -parallel 
spin)  and  triplet  {i.e.,  parallel  spin)  states.  Due  to  the  mass  of  protons  and  their  separation 
distance  in  the  molecules  studied  in  this  section,  the  binding  energy  and  bond  length  for 
singlet  and  triplet  protons  are  approximately  equal  (see  Chapter  5).  Due  to  the  lighter 
mass  of  muons,  compared  to  protons,  their  spin  has  a  significant  impact  on  the  physical 
properties  of  the  di-muonic  hydrogen  molecules  which  were  studied  (see  Table  7-1). 

pppp  is  bound  in  both  singlet  and  triplet  muonic  spin  states  hereafter  designated 
as  pppp(S>  and  ppppn\  respectively.  Bound  excited  vibrational  states  of  pppp  do  not 
exist.  As  can  be  seen  in  Table  7-1,  the  bond  length  of  pppp(T>  is  more  than  four  times 
greater  than  that  of  pppp(S>.  In  its  ground  state,  pppp(S>  has  a  relatively  high  binding 
energy  and  is  very  stable  relative  to  dissociation  into  muonic  atoms  {i.e.,  pp).  In 
comparison  pppp(T>  is  more  weakly  bound.  Setting  the  binding  energy  equal  to  kB  Te, 
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where  kg  is  the  Boltzmann  constant  and  Te  is  temperature,  the  limiting  Te  for  pppp8'  and 
pppp(T>  is  4.76><106K  and  859  K  respectively. 


Table  7-1.  PCI  results  for  pppp  using  2slp  basis  sets.  Binding 
energy  (relative  to  pp)  and  equilibrium  p-p  bond  length  for  different 
particle  spin  states  for  pppp  are  shown.  The  binding  energy  of  the 
ground  rotational  states  (/=0)  of  the  ground  (v=0)  and  first  excited 
(v=l)  vibrational  quantum  states  are  shown.  Only  the  ground 
vibrational  states  are  bound. 


p 

P 

Binding  Energy  (eV) 

A 

2slp 

2slp 

v=0 

V=1 

Bond  Length 

singlet 

triplet 

410.525 

-56.938 

0.005012 

triplet 

singlet 

0.07404 

-737.782 

0.023955 

triplet 

triplet 

0.07404 

-737.722 

0.023955 

Upon  collision  with  a  proton,  hydrogen  atom  (lH)  or  hydrogen  molecule  ( 1 H2), 
pppp8’  and  pppp  T>  can  react  to  form  the  oblate  symmetric  top  molecules  ppppp  or 
pppppe  if  energy  is  efficiently  transferred: 


PPPP  +  P+ - >  PPPPP+  (7.3) 

pppp  +  lH - >pppppe  (7.4) 


or 
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p/LlJLip+{H1 - >  p/Lip/Lipe+  lH 


(7.5) 


The  fusion  rate  constants  as  a  function  of  distance  are  approximately  proportional 
to  e~cx ,  where  c  is  a  constant  and  x  represents  nuclear  bond  length.  As  a  result,  relatively 
small  changes  in  nuclear  separation  can  have  a  large  impact  on  the  rate  of  fusion.  The 
pppp(T>  bond  length  is  more  than  six  times  the  length  of  ppp  (see  Tables  5-1  and  7-1). 
The  p-p  bond  length  in  pppppn>  is  about  three  and  a  half  times  the  length  of  the  p-p  bond 
length  in  ppp+  (see  Tables  5-1  and  7-2).  As  a  result  of  the  relatively  long  bond  lengths  in 
puup'T>  and  pppppn\  fusion  cannot  be  efficient,  even  if  deuterons  or  tritons  replace 
protons  in  the  molecule.  The  bond  length  of  the  triplet  molecules  would  need  to  decrease 
several  fold  for  there  to  be  sufficient  overlap  of  the  nuclear  wave  functions  that  fusion 
would  be  expected  to  occur  at  a  rate  comparable  to  the  single  muon  catalyzed  fusion  rate. 

In  order  for  di-muonic  molecules  with  triplet  muon  spin  to  enhance  the  muon- 
catalyzed  fusion  reaction  rate,  their  reactions  must  result  in  muons  being  freed,  or  in 
species  which  rapidly  lead  to  fusion  being  fonned.  Otherwise  the  muon-catalyzed  fusion 
cycle  will  be  made  inefficient  by  removal  of  muons  from  the  reaction  cycle.  If  pppplT> 
collides  with  a  proton,  or  hydrogen  molecule,  several  reactions  are  possible.  Most  of 
these  result  in  products  and/or  intermediates  that  are  more  closely  bound  being  formed. 
Some  examples  of  possible  reactions  are: 

p/ujup(T)+H2 - >  pjujup(s)  +  AE  =  [-383.2, 14.4]  eV  (7.6) 
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pjUjupiT)  +  p+ - >ppp+  +  pp  AE  =  [—317.5,  —46.8]  eK  (7.7) 


pppp(T)+H1 - >2ppp+  +2e  AE  =  [-575.5,-92.9]  eV  (7.8) 


pppp(T)+H2 - >ppp+  +  pp  +  H*  +e  AE  =  [-301.2, -30. 5]eF  (7.9) 


Mechanisms  by  which  these  reactions  may  occur  are  presented  in  the  discussion  which 
follows. 

The  first  value  of  AE  listed  was  calculated  using  basis  sets  optimized  with  the 
products  and  can  therefore  be  considered  a  lower  bound  of  the  level  of  theory  being  used 
(i i.e HF/FCI  with  2slp  muon  and  proton  basis  sets).  The  second  value  of  AE  listed  was 
calculated  using  basis  sets  optimized  with  the  reactants,  and  can  be  considered  an  upper 
bound.  These  reactions  are  exothennic.  The  rate  they  occur  depends  on  the  reaction  path 
and  on  the  rate  of  competing  reactions.  Some  possible  reaction  paths  are: 


pppp(T)  +  p+ 


-> 


[ppppp+f]  ->  ppp+  +  pp 


(7.10) 


pppp(T)  +H2  — >  pppppe{T)  +H *  — »  pppe[T)  +  pp  +  H * 
— »  ppp+  +  pp  +  H*  +e~ 


and 
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pppp(T)  +  H2  — >  ppppppee(T)  — >  pjupe{T)  +  pjupe(T)  — >2 ppp+  +2e  (7.12) 


Some  possible  competing  reactions  are: 


pppp(T)  +H2  — »  p  p  p  p  ppee{,)  — >  ppppppe+  +e  (7.13) 


and 


pppp(T)  +  H2  — »  pppppe  +H*  — >  ppppp+  +H*  +e  (7.14) 

In  order  for  a  reaction  to  proceed  rapidly,  an  efficient  mechanism  for  the  transfer 
of  energy  must  operate.  Energy  may  be  transferred  or  released  through  collisions,  Auger 
transitions,  or  high-energy  x-ray  emission.  There  are  limits  to  the  amount  of  energy  that 
can  be  transferred  by  collisions,  Auger  transitions  or  x-ray  emission.  Collisional 
relaxation  is  only  effective  if  the  reaction  coordinate  couples  energy  into  particle  motion; 
reactions  that  rely  on  efficient  collisional  energy  elimination  usually  proceed  more 
rapidly  when  the  difference  in  energy  between  reactants  and  products  is  comparable  to 
relevant  vibrational  modes  accessed  in  the  reaction. [6;  7]  For  larger  energy  differences 
Auger  electron  transitions  are  more  efficient.  For  example,  t/ut+  can  go  from  an  excited 
rotational-vibrational  state  (Jy)  =  (1,1)  to  the  rotationally  excited,  ground  vibrational 
state  (1,0)  via  the  release  of  an  Auger  electron  ( t  represents  a  triton).  This  transition 
yields  an  energy  release  of  approximately  244  eV.[8]  The  somewhat  similar  reaction, 
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(7.12)  where  the  triplet,  ground  state  vibrational  and  rotational  molecule  pppppe 
transitions  to  the  singlet  ground  rotational  and  vibrational  state  of  ppppp  via  an  Auger 
transition  is  unlikely  to  proceed  at  a  significant  rate.  A E  for  the  reaction  is  about  -1.00 
keV,  and  the  energy  of  an  Auger  electron  would  need  to  be  about  718.5  eV  (see  Figure  7- 
1  and  Table  7-2).  To  a  first  order  approximation,  an  electron  can  be  expected  to  interact 
with  (ppppp+)(S>  and  (ppppp+)<T>  in  a  manner  similar  to  how  it  would  interact  with  a  triton. 
ppppp+  ions  are  small  with  a  concentrated  +1  charge  and  a  mass  similar  to  that  of  tritium. 
Being  that  Auger  transitions  with  energies  greater  than  700  eV  normally  do  not  occur 
with  nuclei  that  have  an  atomic  number  smaller  than  about  9,  the  energy  of  this  transition 
is  larger  than  could  be  expected  to  occur.[9;  10]  The  reaction  path  shown  in  reaction 
Equation  1.11  is  unlikely  to  proceed  at  an  appreciable  rate  for  ground  rotational  states,  for 
the  same  reason,  (note:  A E  and  the  energy  of  the  Auger  transition  was  calculated  using 
basis  sets  optimized  for  ppp+).  The  only  way  these  reactions  could  be  expected  to 
proceed  at  appreciable  rates  is  if  rapid  transitions  to  excited  rotational  states  occur. 
Figure  7-1,  and  Tables  7-2  and  7-3  compare  the  total  energy  ( i.e energy  compared  to 
free  unbound  particles)  of  several  muonic  hydrogen  molecules  and  molecular 
combinations.  For  the  calculations  presented  in  Figure  7-1,  and  Tables  7-1  and  7-3,  basis 
sets  and  bound  lengths  were  optimized  simultaneously  for  pp,  p/ip  \  pppp<S),  and  pppp(T> 
using  2s  lp  basis  sets  and  FCI  methods.  The  (ppppp+)<s>  and  {ptipup  )<T>  bond  lengths 
were  optimized  using  the  basis  sets  indicated  and  Cl  methods  with  8  active  muon  and  8 
active  proton  molecular  orbitals.  As  can  be  seen  from  Figure  7-1  and  Table  7-3, 
(ppppp+)(S)  has  bound  rotational  states  that  are  relatively  close  in  energy  to  those  of 
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(piupiup+)(T) ,  therefore,  rapid  transitions  to  the  (2,1)  and  (3,3)  rotational  states  of 
(pjuppp+)(S>  are  possible. 


Table  7-2.  Cl  results  for  ppppp+  with  singlet  and  triplet  particle 
spin  in  ground  vibrational  and  rotational  states.  Binding  energy 
(relative  to  pp)  and  equilibrium  p-p  bond  length  for  different 
particle  spin  states  are  shown.  Only  the  ground  vibrational  states 
are  bound. 


P 

Binding  Energy 

(eV) 

Bond  Length 

(A) 

singlet 

triplet 

961.773 

.005792 

triplet 

singlet 

243.208 

.01212 

triplet 

triplet 

243.208 

.01212 

Some  possible  avenues  to  enhance  fusion  indicated  by  these  calculations  are 
processes  that  connect,  with  small  energy  differences,  to  cation  species  that  can  relax  by 
Auger  relaxation,  for  example: 


pppp(T){0,®)  +  p+ - >{ppppp+ )(n(3,3)  (7.15) 


(ppppp+){r)  (3, 3)  {ppppp+fT)  (2, 1)  -A  ppp+  (1)  +  pju  (7.16) 


and 


{ppppp+){T\ 3,3)  {ppppp+){T\ 0, 0)  -A  ppp+{ 0)  +  pp  (7.17) 
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Table  7-3.  Binding  energy  with  respect  to  pp  for  the  bound  rotational  states  (/)  of 
ppp  and  (J,K)  of  pppp  and  ppppp  (see  Figure  7-1).  The  rotational  states  shown  all 
correspond  to  ground  vibrational  states.  There  are  no  excited  vibrational  states  of 
the  muonic-molecules  shown.  The  binding  energy  of  ppp+  was  taken  from 
reference.[8]  (S)  represents  singlet  muonic  spin  states  (i.e.,  anti-parallel  muon  spin) 
and  (T)  represents  triplet  muonic  spin  states  (i.e.,  parallel  muon  spin). 


(J) 

Bine 

ing  Energy  (eV) 

UK) 

Binding  Energy  (eV) 

JS) 

PPPP 

T) 

PPPP 

PPP 

(PPPPP+)W 

(PPPPP')'0 

(0) 

410.525 

0.07404 

253.152 

(0,0) 

961.773 

243.208 

(1) 

336.868 

107.266 

(1,0) 

745.275 

193.765 

(2) 

189.554 

(U) 

802.128 

206.749 

(2,1) 

370.865 

108.258 

(2,2) 

540.266 

146.945 

(3,3) 

176.716 

63.919 

Reactions  which  produce  pp  as  a  product  can  transfer  a  significant  amount  of 
energy  to  the  bound  muon.  The  quantum  energy  levels  of  the  bound  muon  in  pp  can  be 
calculated  from 


=  (7.18) 

2/7  n 

Where  Z  is  the  charge  on  the  proton  (+1),  p  is  the  reduced  mass  between  the  proton  and 
the  muon,  h  is  the  reduced  Planck’s  constant  and  n  is  the  principle  quantum  number  (see 
Figure  7-2).  When  forming  muonic  hydrogen  atoms,  negative  muons  will  most  often 
initially  form  a  muonic  atom  in  which  the  principle  quantum  number  (n)  is  very  large 
(i.e.,  >  14).  As  the  muon  loses  energy,  transitioning  from  one  state  to  another, 
characteristic  x-rays  are  emitted. 
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Figure  7-1.  Binding  energy  with  respect  to  /?//  of  the  bound  rotational  states  (J)  of 
ppp+  and  the  bound  rotational  states  ( J,K)  of pppp  and ppppp+  (see  Table  7-3).  Both 
singlet  and  triplet  muonic  spin  states  are  shown.  The  rotational  states  shown  all 
correspond  to  ground  vibrational  states.  There  are  no  excited  vibrational  states  of 
the  muonic-molecules  shown.  The  binding  energy  of  ppp+  was  taken  from 
reference.  1 8]  The  basis  sets  and  bond  lengths  of  pppp  and  pppp(T>  were  optimized 
using  2slp  basis  sets  and  FCI  methods.  The  ppppp+  and  (ppppp+)(T)  bond  lengths 
were  optimized  using  optimized  ppp+  basis  sets  and  Cl  methods  with  8  active  muon 
and  8  active  proton  molecular  orbitals.  The  ground  state  rotational  and 
vibrational  energy  of  ppppp+  and  (ppppp+/T>  was  calculated  using  the  indicated 
bases  sets  and  bond  lengths  at  the  FCI  level.  The  ground  rotational  state  (0,0)  does 
not  exist  in  the  ground  vibrational  state  of  triatomic  oblate  symmetric  top 
molecules  [ i.e .,  (ppppp+1  and  (ppppp+/T> ].[5;  11;  12;  13]  (1,1)  is  the  lowest 
obtainable  energy  level  for  these  molecules. 
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Figure  7-2.  The  first  few  muonic  energy  levels  of /?//.  The  ground  state  energy  [£(1)]  is  - 
2.52849  keV. 

7.4  Conclusions 

Di-muonic,  di-hydrogen  molecules  can  exist  in  singlet  and  triplet  muon  spin 
states.  The  triplet  molecules  are  very  weakly  bound,  and  due  to  their  bond  length,  cannot 
be  expected  to  fuse  at  an  appreciable  rate. 

Reaction  paths  by  which  triplet  pppp  (i.e.,  pppp<T))  can  be  transformed  into  singlet 
pppp  {i.e.,  pppp(S>)  or  single  muon  hydrogen  molecules  were  studied.  All  of  the  viable 
reaction  paths  involve  collisions  with  other  molecules  and/or  the  fonnation  of  rotationally 
excited  molecules.  Singlet  pppp  has  three  bound  rotational  states  (J  =  0,  1,  2).  pppp<T>  is 
only  bound  in  the  ground  rotational  state.  The  equilibrium  bond  lengths  of  pppp(T>  and 
pppp(T)  are  5.012  and  23.955  mA  respectively  (see  Table  7.1). 
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Both  singlet  and  triplet  (muon  spin)  forms  of  ppppp+  exist.  Both  of  these  forms 
of  the  molecule  have  five  bound  rotationally  excited  states:  (1,1),  (1,0),  (2,2),  (2,1),  and 

(3,3).  The  equilibrium  positions  of  (ppppp+YS^  and  ( ppppp+){T)  nuclei  form 
equilateral  triangles  with  bond  lengths  of  5.792  and  12.12  mA  respectively  (see  Table 
7.2). 

No  vibrationally  excited  fonns  of  ptiup  or  ptipup  exist.  If  heavy  isotopes  of 
hydrogen  are  substituted  for  the  protons  in  these  molecules,  bound  excited  vibrational 
states  may  exist. 
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VIII.  Di-Muoic  Hydrogen  Reaction  Kinetics 


8.1  Introduction 

Determining  a  mechanism  by  which  muon-catalyzed  fusion  can  be  used  as  a  pure 
fusion  source  of  energy  has  so  far  been  elusive.  A  study  of  current  experimental  and 
theoretical  literature  leaves  little  hope  that  mono-muonic  reaction  mechanisms  can  result 
in  sufficiently  high  fusions  yields  for  this  process  to  be  used  as  a  pure  fusion  source  of 
energy  (see  Chapter  2  for  a  review  of  this  literature).  If  a  mechanism  is  to  be  found, 
which  results  in  significantly  higher  yields  than  has  thus  far  been  observed,  novel 
reaction  paths  must  be  investigated.  As  of  the  writing  of  this  document,  most  studies  of 
muon-catalyzed  fusion  have  dealt  with  single  muon  hydrogen  molecules. [1;  2;  3;  4;  5;  6; 
7;  8;  9]  In  this  dissertation  a  novel  approach  to  muon-catalyzed  fusion  is  investigated. 
The  possibility  that  di-muonic  molecular  reactions  can  enhance  the  overall  fusion  yield  is 
considered. 

Di-muonic  molecules  have  not  contributed  significantly  to  the  observed  reaction 
rate  of  muon-catalyzed  fusion  in  any  of  the  experimental  studies  performed  to  date.  Few 
theoretical  studies  of  these  molecules  exist  in  the  literature,  and  no  studies  have  been 
published  which  analyze  the  effects  of  di-muonic  molecules  on  the  muon-catalyzed 
fusion  reaction  rate. [10]  This  can  be  attributed  to  the  low  concentration  of  near-thermal 
nruons  (i.e.,  low  muon  flux)  that  generate  hydrogenic  muonic  molecules  used  in  most 
muon-catalyzed  fusion  experiments  performed  to  date. [10]  For  a  given  reaction 
chamber,  as  the  concentration  of  thennal  muons  increases  the  contribution  of  di-muonic 
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molecular  reactions  will  increase.  To  the  author’s  knowledge,  the  concentration  of 
thennal  muons  required  to  significantly  influence  muon-catalyzed  fusion  reaction  rates, 
via  di-muonic  reactions,  has  not  been  published.  A  better  understanding  of  potential 
mechanistic  pathways  that  may  contribute  to  muon-catalyzed  fusion  is  needed  to  estimate 
whether  di-muonic  molecular  reactions  can  significantly  contribute  to  enhance  muon- 
catalyzed  fusion  kinetics.  Increasing  the  fonnation  rate  of  di-muonic  hydrogen  molecules 
will  not  increase  the  fusion  rate  unless  these  species  efficiently  form  and  fuse,  or  provide 
a  rapid  pathway  to  enhance  the  formation  of  d/j,d+ ,  duf,  or  t/it' ,  promoting  ordinary 
muon  catalysis.  As  was  shown  in  Chapter’s  4  and  5,  the  average  distance  between  nuclei 
in  di-muonic,  di-nuclear  hydrogen  molecules  is  greater  than  in  the  corresponding  mono- 
muonic  molecules.  Di-muonic  hydrogen  molecules  can  exist  in  singlet  ( S)  and  triplet  (7) 
muonic  spin  states.  The  triplet  molecules,  which  are  more  likely  to  form,  cannot  have  a 
significant  fusion  rate,  due  to  the  long  bond  lengths  of  these  molecules  (see  Chapter  7). 
Mechanisms  by  which  these  triplet  molecules  may  be  transfonned  into  states  or 
molecules  which  may  have  appreciable  fusion  cross-sections  were  discussed  in  Chapter 
7. 

In  previous  chapters  physical  properties  of  di-muonic  hydrogen  molecules  and 
reactions  they  could  undergo  have  been  studied.  In  this  chapter  the  reaction  kinetics  of 
paths  that  could  lead  to  the  formation  of  these  molecules  and  lower  bounds  on  the  muon 
flux  needed  to  form  them  in  significant  quantities  will  be  presented.  Although  it  is 
currently  unknown  if  the  reaction  of  di-muonic  molecules  will  enhance  or  decrease  the 
overall  reaction  yield,  by  studying  potential  reaction  paths  a  lower  bound  on  the  muon 
flux  required  to  yield  significant  di-muonic  molecular  effects  can  be  established. 
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Additionally,  it  is  possible  to  detennine  an  upper  bound  on  possible  yield  enhancement 
resulting  from  di-muonic  muon-catalyzed  fusion  pathways. 

8.2  Methods 

In  order  for  di-muonic  hydrogen  molecules  to  have  a  significant  impact  (negative 
or  positive)  on  the  rate  of  muon-catalyzed  fusion,  a  substantial  quantity  of  di-muonic 
molecules  must  be  present  in  the  reaction  chamber.  Additionally,  the  rates  these 
molecules  fuse,  or  go  through  a  series  of  reactions  leading  to  fusion,  must  be  significantly 
different  than  the  rates  of  competing  single  muon  reactions  and/or  the  sticking  fractions 
of  these  reactions  must  be  significantly  different  than  competing  single  muon  reactions. 

The  di-muonic  hydrogen  formation  rate  (iy,)  can  be  estimated  using  the  hard 
sphere  collision  frequency  (zAB) 


z 


AB 


=  nAnBno 


2 

AB 


'  %kBT  V 

V  71  ^AB  A 


(8.1) 


no 


AB 


8  kBT 
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Concentration  of  molecules  A  and  B  respectively 
Collision  cross-section 
Reduced  mass  of  particles  A  and  B 
Temperature  (K) 

Boltzmann’s  constant 
Activation  energy 
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An  upper  bound  of  the  reaction  rate  occurs  when  all  collisions  result  in  a  reaction.  In  this 

case,  Eact  =  0  and  zAB  =  rhi. 

In  order  to  calculate  the  di-muonic  hydrogen  fonnation  rates  there  are  three 
species  that  can  collide  which  need  to  be  considered:  fx,  xju,  and  (xfix)xee,  where  x 

j  2  3  12  3 

represents  H,~H,  and  //nuclei.  A  represents  H,  H,  and  '//atoms.  Following  are  some 
reactions  which  can  result  in  the  formation  of  di-muonic  hydrogen  molecules  that  will  be 
considered  in  this  chapter. 


(xjLix)xee  +  fi  - >x[ijux  +  X * 


(8.3) 


(xjLix\xee  +  ju  - >xjux/jx 


(8.4) 


XjLl  +  XjLl - >XJUJUX 


(8.5) 


(. xjux)xee  +  xju - >xJujux  +  X1 


(8.6) 


[xjux)xee  +  xn - >(x jux /ux)xee 


(8.7) 


( x  jux)  xee  +  (x  jux)  xee - >xjujux  +  2X2 


(8.8) 


and 
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[x/ux)xee  +  {x/ux)xee - >xjuxjux+ ljX2 


(8.9) 


These  reactions  cannot  form  products  efficiently,  even  though  they  are 
exothermic,  unless  energy  is  rapidly  dissipated  from  the  collision  complex.  In  the  cases 
of  reactions  8.3,  8.6,  8.8  and  8.9  the  non-muonic  product  could  carry  off  this  excess 
energy.  In  the  case  of  reactions  8.4  and  8.7,  the  energy  may  be  dissipated  by  Auger 
processes.  Reaction  8.5  could  occur  if  an  excited  rotational  and/or  vibrational  state  of 
x/i/vc  exists  in  which  the  energy  of  reaction  (A E)  is  dissipated  between  rotational  and/or 
vibrational  modes.  Such  a  reaction  mechanism  can  only  exist  if  A E  is  small. [4;  5]  If 
excited  vibrational  and/or  rotational  states  do  exist,  which  make  this  reaction  possible,  the 
reverse  reaction  rate  is  expected  to  be  close  to  that  of  the  reaction  rate.  Therefore  the 
only  potential  mechanisms  for  reaction  8.5  to  occur  involve  collisions  with  a  third 
body.[l  1] 

As  can  be  seen  from  Equations  8.1  and  8.2,  the  reaction  rate  depends  on  the 
number  of  collisions  and  upon  the  effectiveness  of  those  collisions.  The  effectiveness  of 
the  collisions  depends  on  the  orientation  and  energy  of  the  colliding  species.  It  is 
common  for  muons  to  be  input  into  a  reaction  chamber  with  energy  greater  than  will 
react  to  fonn  x/x.  When  this  is  the  case  muons  must  be  moderated  through  collisions 
before  they  can  react.  If  the  energy  of  effective  collisions  between  /X  and  Xi  is 
significantly  lower  than  that  required  to  fonn  products  when  muons  collide  with  x/vcx  the 
formation  rate  of  the  di-muonic  molecules  could  be  enhanced.  If  the  energy  of  effective 
collisions  between  /£  and  Xi  is  significantly  higher  than  that  required  to  form  products 
when  muons  collide  with  x/vcx  the  formation  rate  of  the  di-muonic  molecules  via  these 
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reaction  paths  would  be  slowed.  The  size  of  the  difference  in  the  effective  collision 
energy  of  these  molecules  is  unknown. 

The  maximum  reaction  rate  will  occur  when  the  total  number  of  muons  are 
equally  divided  between  the  colliding  species.  Dividing  all  of  the  muons  in  a  reaction 
chamber  equally  between  reactants  can  be  difficult  or  impossible.  Being  that  the  results 
of  such  calculations  result  in  upper  bounds  to  the  reaction  rate,  there  is,  however,  value  in 
perfonning  the  calculations,  even  if  these  reaction  conditions  cannot  be  obtained.  Any 
reaction  which  has  an  upper  bound  of  its  reaction  rate  which  is  too  low  to  significantly 
affect  the  overall  fusion  rate  can  be  neglected,  and  further  analysis  of  this  reaction  path 
can  be  ignored.  Taking  this  into  account,  upper  bounds  of  reaction  rates  were  calculated 
in  this  way  (see  Table  8-1). 

Table  8-1  lists  the  maximum  reaction  rate  (i.e.,  collision  frequency)  for  the 
reactants  in  Equations  8.3  through  8.7.  The  collision  frequency  is  calculated  for  a  muon 
flux  of  1.5  x  10 16  negative  muons  per  second.  This  negative  muon  flux  was  chosen  for 
the  calculations  because  it  is  the  highest  continuous  flux  known  to  be  available  at  a 
currently  operating  accelerator^  12:34] 
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Table  8-1.  Theoretical  maximum  reaction  rate  obtainable  with  a 
negative  muon  flux  of  1.5  x  1016  muons/s  focused  into  a  volume  of 


lcm3.  A  +  B - > Products.  An  excess  of  H2  in  the  reaction 

chamber  is  assumed. 


A 

B 

Max  reaction 
rate  at  300  K 
(s'1) 

Max  reaction 
rate  at  1000  K 
(s'1) 

Max  reaction 
rate  at  4000  K 
(s'1) 

PPPP 

PP 

9.24xl021 

1.69xl022 

3.37x1 022 

PP 

PP 

7.72xl017 

1.41xl018 

2.82xl018 

PPPP 

P 

2.53xl022 

4.61x1 022 

9.23xl022 

Due  to  the  low  collision  frequency  of  pp  +  pp  and  the  energy  dissipation 
problems  discussed  previously,  pp  +  pp  reactions  will  not  be  considered  further. 

Upper  bounds  for  reaction  rate  constants  (7)  for  reactions  8.3  through  8.9  can  be 
calculated  ( i.e Eact  =  0)  using  a  hard  sphere  model, 


/l 


=  no , 


8  kBT 


i 


Was) 


Eac, 

"Ur 


(8.10) 


These  values  vary  slightly,  depending  on  the  hydrogen  isotopes  being  considered  (see 
Table  8-2).  The  reaction  cross  section  (cr)  was  taken  as  the  sum  of  the  hard-sphere  radii 
of  the  reactants.  The  hard-sphere  radii  for  xpi  was  taken  as  hard-sphere  radii  of  xe. 


divided  by  the  reduced  mass  ratio 


f  u  ^ 

r^xfi 


\  /U  ) 


.  The  hard-sphere  radii  for  (. x/ux)x  was  assumed 


to  be  equal  to  the  radii  of  similar  size  molecules  X2. 
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Table  8-2.  Upper  bounds  of  the  reaction  rate  constants  for  reactions  8.3  through 
8.6  at  300  K  and  liquid  hydrogen  density.  All  hydrogen  isotope  combinations  are 

shown.  A  +  B - »  Products 


A 

B 

xao'V) 

A 

B 

XOO'V) 

PsP 

P 

69.7 

pd2p 

P 

69.2 

p3p 

pp 

25.5 

pd2p 

PP 

24.1 

P3P 

dp 

20.6 

pd2p 

dp 

18.9 

P3P 

tp 

18.5 

pd2p 

tp 

16.6 

d3p 

P 

69.1 

pt2p 

P 

69.0 

d.3p 

pp 

23.8 

pt2p 

PP 

23.5 

d3p 

dp 

18.4 

pt2p 

dp 

18.1 

d3p 

tp 

16.1 

pt2p 

tp 

15.7 

hp 

P 

68.9 

d2tp 

P 

69.0 

t3p 

PP 

23.2 

d2tp 

PP 

23.5 

Up 

dp 

17.6 

d2tp 

dp 

18.1 

t3fi 

tp 

15.2 

d2tp 

tp 

15.7 

P2dp 

P 

69.4 

dt2p 

P 

68.9 

P2dp 

PP 

24.7 

dt2p 

PP 

23.3 

P2dp 

dp 

19.5 

dt2p 

dp 

17.8 

P2dp 

tp 

17.3 

dt2p 

tp 

15.4 

p2tp 

P 

69.2 

d2tp 

d2tp 

45.7 

p2tp 

PP 

24.1 

d2tp 

dt2p 

47.2 

p2tp 

dp 

18.9 

dt2p 

dt2p 

48.8 

p2tp 

tp 

16.6 

xjujux  can  form  in  singlet  or  triplet  muonic  spin  states  depending  on  the  amount  of 
excess  energy  carried  off  by  X’or  X2  (see  Equations  8.3,  8.6,  8.8  and  8.9,  and  Figure  7- 
1).  The  singlet  spin  state  of  this  molecule  may  have  a  rapid  fusion  rate  when  heavy 
hydrogen  isotopes  are  involved  ( e.g X >  1  f)7.s'  1 ).  Triplet  state  molecules  cannot  have 

rapid  fusion  rates  due  to  the  large  nuclear  bond  distances  involved  in  these  molecules, 
(see  Section  7.3  of  this  document). 

In  order  for  (xpxpx)xee  to  fonn,  via  Equation  8.7,  the  binding  energy  may  be 
divided  between  excited  vibrational  and/or  rotational  quantum  states.  This  can  only 
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occur  for  reactions  with  very  small  values  of  AZf. [  1 3]  As  a  result,  ppppp+  cannot  form 
with  anti-parallel  muon  spin  via  reaction  8.6  (see  Figure  7-1)  (note:  The  binding  energy  is 
calculated  to  be  too  large  for  an  Auger  transition  to  carry  off  all  of  the  binding  energy). 
When  heavier  isotopes  of  hydrogen  are  involved,  singlet  molecules  similar  to  ppppp+ 
may  form  due  to  the  existence  of  weakly  bound  excited  vibrational  states.  The  singlet 
states  of  this  molecule  may  have  rapid  fusion  rates,  but  triplet  muon  spin  states  of  this 
molecule  cannot  have  rapid  fusion  rates  due  to  the  magnitude  of  the  nuclear  bond  lengths 
(see  Chapters  4  and  5) 

Although  xppx'1'  and  xgxgx'n  cannot  have  rapid  fusion  rates  (e.g.,  A  >  1 07,v  1 )  due 

to  the  relatively  large  nuclear  bond  length  of  these  molecules;  these  molecules  may  react 
and  form  molecules  with  significantly  shorter  bond  lengths.  These  reactions  are 
discussed  in  detail  in  Chapter  7. 

In  Appendix  A  kinetic  rate  equations  for  single  muon,  muon-catalyzed  fusion  are 
developed.  These  equations  can  be  simplified  considerably,  when  x  is  allowed  to 
represent  all  hydrogen  isotopes  for  which  a  given  reaction  is  possible.  In  the  discussion 
which  follows,  it  is  assumed  that  the  reaction  conditions  are  such  that  the  free  radicals 
H” ,  D‘ ,  and  T‘  do  not  exist  in  significant  quantities  (e.g.,  temperature  <  2000  K).  Xo 
represents  the  decay  rate  constant  of  negative  muons  (A0  =4.55xlOA_1).[3]  In  the 

equations  which  follow,  electrons  and  products  whose  formulation  doesn’t  have  a 
significant  impact  on  the  overall  muon-catalyzed  fusion  rate  have  not  been  shown.  The 
simplified  reaction  equations  are: 
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A.res 

XjLl  +  X0  <  r"  >  XjLlXX 


(8.12) 


XJU  +  X2 XjLlXX 


(8.13) 


and 


XjLixx — »  fusion  products  +  ju  (8.14) 

Xjuxx  can  form  via  reversible  resonant  pathways  (i.e.,  Equation  8.12)  and  non-resonant, 
non-reversible  pathways  {i.e.,  Equation  8.13).  Which  of  these  reactions  predominates  or 
significantly  contribute  to  the  overall  fusion  yield  depends  on  which  hydrogen  isotopes 
are  involved  in  the  reactions. 

It  is  common  practice  to  combine  the  rate  constants  that  affect  the  concentration 
of  xjuxx  to  form  an  effective  rate  constant  (Ay)  that  more  closely  represents  what  is 
observed  experimentally.  For  the  resonance  pathways  this  is  done  by  taking  the  rate 
equation 


8  [x//xr] 


Xu  [Xp\  [X2  ]-(r  +  kr)  [xpxx\ 


8 1 


(8.15) 


and  setting  it  equal  to 


f  [x//xr] 


8t 


=  kb  [x//]  [X2  ]  -  kc  [x//xx] 


(8.16) 


which  corresponds  to  the  effective  reaction  equation 


XjU  +  X2 — » XjUXX 


(8.17) 


where  x™ ,  2"  ,  Y ,  ka ,  kh,  and  kc  are  rate  constants  defined  in  Equations  8.12  thru  8.14. 

0 1  Jt  LLXX 

Using  the  steady  state  approximation  — - -  =  0  and 

dt 


k  -  2 res 

Kb  ~  \fi 


\Y  +  kcJ 


(8.18) 


When  non-resonant  paths  and  different  spin  states  of  xp  (F)  and  xpxx  (5)  are  added  to  the 
equation, 
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K  = 


2;;+ 1 2 


res 

'x/j,FS 


S 


_ K 

Xr„ 

V  F 


+  k 


J 


(8.19) 
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Then  it  follows  that  the  reaction  rate  equations  for  these  effective  reaction  paths 


are: 


=  K[X 2  ]  [/T  ]  -  [*//]  [kb  [x2  ] + Aq  }  (8.20) 

d^XX^  =  kb  [X2 ][x/u\-(kc+\) ) [xx/ux]  (8.21) 

^  ^=kc{\-Bca) [xnxx\ - [/T ] { ka  [X2 ]  +  A0 }  (8.22) 


and 


3  [  product 


dt 


=  kc  [x//xx] 


(8.23) 


where 


//  =  ntotal  ~[xju]-[xjuxx]- Byproducts] 


(8.24) 


and 
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(8.25) 


[/U/]  =  Me_v 

The  He  path  branching  ratio  B  equals  1 ,  except  for  d-d  fusion.  The  He  sticking  fraction  is 

co  and  Bo  represents  the  initial  muon  concentration  in  a  pulsed  system. 

If  a  continuous  muon  flux  (O)  is  added  to  the  reaction  chamber,  then 


^  Btotal  _  .  r  _  n 

ot  °~AoL^J 


(8.26) 


3  fit  } 

A  steady  state  (SS)  of  /u~  ,  can  be  determined  by  considering  Lim  L  lo'“  J  =  0  and 

00  dt 


consequently, 


[^ta/]ss=  Y  (8-27) 

/L0 

In  the  case  of  d-t  fusion  the  above  simplified  equations  assume  near  equal  concentrations 
of  deuterium  and  tritium  and  neglect  many  of  the  less  significant  reaction  paths.  A  more 
complete  analysis  could  be  perfonned  using  the  equations  presented  in  Appendix  A  and 
the  isotope  concentration  calculation  methods  presented  in  Appendix  B.  Using  the 
simplified  equations  a  maximum  muon-catalyzed  fusion  yield  of  160  fusions/muon  is 
calculated  if  only  single-muon  catalysis  occurs  in  significant  amounts.  This  compares  to 
an  observed  experimental  value  of  150  fusions/muon.[2;  3] 
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The  following  di-muonic  reactions  can  be  added  to  Equations  8.1 1,  8.14  and  8.17: 


XjUXX  +  jU  - - — >  XjUjUX 


(8.28) 


XJUXX  +  jU  — XJUXJUX 


(8.29) 


X  jUXX  +  XjLl - L — >XfifiX  +  X1 


(8.30) 


XjUXX  +  XjU - - — >  X  jLLX  jLLX  +  X 


(8.31) 


X/H/UX+  X2 — — >  X  /.LX  /.LX 


(8.32) 


xppx — — >  fusion  products +  2  p 


(8.33) 


k; 


xjuxjux  +  Xj  — ■'— >  2xpx 


(8.34) 


xpx — - — >  fusion  products  +  p 


(8.35) 


xpxpx — k — >  fusion  products  +  2 p 


(8.36) 


141 


2x  /uxx  — x  jujux  +  2x2 


(8.37) 


and 


2  XjUXX - ^-^XjUXjUX  +  \\x2 


(8.38) 


The  rate  equations,  with  the  di-muonic  reaction  paths  included  become: 


<3[x//] 

dt 


K  [X2  ]  [V  ]  -  [xp]  {kb[X2]  +  ; L0 }  -  (kf  +  kg )  [x/jxx]  [x/u] 


(8.39) 


C  [x/AXx] 


dt 


[X2][x//]-(&c  +  20)[xx jux\-{k  f  +  £g)[x//xx][x//] 
-{km+kn)[xpixxf 


(8.40) 


<3[x//x] 

dt 


-(20  +  k,  j[xjux]  +  2k j  [x//x//x][X2] 


(8.41) 


d  [x////x] 
dt 


-  ( A0  +  kj )  [x////x]  -  kh  [x////x]  [x2  ]  +  kf  [x//]  [x//xx] 
+  km  [x//xx]2 


(8.42) 
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It  may  be  noted  that  the  terms  kd  [x//xx][//  J  and  ke[xjuxx][ju  ]  have  not  been 
included  in  the  above  equations.  This  is  because  these  terms  are  negligible  when 
compared  to  the  terms  kf  [x//xx]  [x//]  and  kg  [x//xx]  [x//] .  If  rate  constants  significantly 

different  from  those  used  in  this  study  are  used,  then  the  appropriateness  of  neglecting 
these  tenns  should  be  reassessed.[l;  5] 

The  probability  of  xjuxx  +  xjuxx  collisions  resulting  in  the  formation  of  di-muonic 
molecules  depends,  in  part,  on  the  binding  energy  of  the  xpx  portion  of  the  molecule. 
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The  possibility  of  these  collisions  resulting  in  the  fonnation  of  di-muonic  molecules  only 
exists  when  the  xjux  binding  energy  is  very  small  (e.g.,  ~0.5  eV)  as  it  is  in  the  case  of  the 
excited  vibrational  and  rotational  state  (1,1)  of  djut.[  14]  dgd  and  till  are  bound  too  tightly 
to  expect  these  collisions  to  result  in  the  formation  of  significant  quantities  of  di-muonic 
molecules  via  this  reaction  path.  In  the  case  of  dgt  collisions  it  is  unknown  if  Eact  for  this 
reaction  is  low  enough  for  these  collisions  to  contribute  significantly  to  the  overall  fusion 
yield.  For  this  reason,  these  reactions  have  been  considered  in  the  discussion  which 
follows. 

The  reaction  rate  constants  used  for  the  calculations  presented  in  this  chapter  have 
been  taken  from  references  [1]  and  [5],  or  are  presented  in  Table  8-2  of  this  document. 
The  rate  constants  calculated  in  this  section,  as  well  as  those  taken  from  other  references 
have  all  been  normalized  to  liquid  hydrogen  density  and  300  K  ( i.e .,  4.25  x  10" 
molecules/cm  ).  The  rate  equations  presented  in  this  chapter  (i.e.,  k  and  /.)  can  be 
adjusted  to  account  for  pressure  [i.e.,  k(P )  and  A(P)]  by  multiplying  the  rate  constants 
which  depend  on  pressure  (i.e.,  those  which  are  multiplied  by  the  concentration  of  A3  in 
the  rate  equations)  by  5.75374  xlO'4  times  the  pressure  in  atmospheres.  This  includes  the 
effective  rate  constant  kb,  since  k,  and  T FS  are  independent  of  pressure  and 


K<p)-K'/p)+Ij^Zfs<p) 


=  5.75374x10  4 Pressure 


(  \ 

K 

X  r  fs  +k< 

V  F 

n  nr  .  X  1  n  res 

4h  +  2j/v. 


FS 


f  \ 

K 

TjTFS+kc 

V  F 


(8.47) 


=  5.75374x10  4  Pressure  k 
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8.3  Results  and  Discussion 


This  chapter  examines  reaction  paths  that  may  lead  to  the  formation  of  di-muonic 
hydrogen  molecules  and  the  muon  flux  needed  to  form  these  molecules  in  significant 
quantities.  The  reaction  paths  studied  are  limited  to  systems  containing  only  deuterium 
and/or  tritium.  The  reactions  conditions  assume  protium  concentrations  to  be  negligible. 
The  muon  flux  required  to  significantly  affect  the  muon-catalyzed  fusion  yield  depends 
on  which  reaction  paths  predominate  and  which  isotopes  are  present. 

The  muonic  hydrogen  molecules  in  greatest  concentration  in  a  muon-catalyzed 
fusion  chamber  are  xjuxx.  As  a  result,  if  xjuxx  +  xjuxx  collisions  are  efficient  at  forming 
di-muonic  hydrogen  molecules  these  reactions  will  predominate  relative  to  other  di- 
muonic  hydrogen  formation  paths.  As  was  discussed  in  the  previous  section,  due  to  the 
binding  energy  and  Coulombic  repulsion  of  the  xjux  portion  of  these  molecules,  d/udd  and 
tpitt  molecules  are  not  likely  to  have  large  reaction  rates  for  this  reaction  path  \i.e.,  the 
activation  energy  (Eact)  is  too  large].  As  a  result,  these  reaction  paths  can  be  neglected 
for  all  molecules,  except  dfxtx.  By  assuming  Eact  =  0  ( i.e all  dpitx  +  dpitx  collisions 
result  in  di-muonic  molecules  being  formed),  a  lower  bound  on  the  muon  flux  required  to 
significantly  affect  the  overall  fusion  rate  can  be  detennined.  If  this  effect  will  be 
positive  or  negative  depends  on  the  fusion  rate  of  the  molecules  formed  (see  Chapter  7). 

If  Eact  =  0  for  the  reactions  discussed  above,  and  the  di-muonic  molecules  fonned 
quickly  go  to  fusion  products,  then  a  continuous  negative  muon  flux  as  low  as  105  could 
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result  in  di -muonic  reactions  increasing  the  overall  fusion  yield  by  as  much  as  10%  ,  and 
a  flux  as  low  as  3  x  105  could  double  the  total  yield  (see  Figure  8-1). 


Maximum  Increased  Yield 
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Figure  8-1.  Maximum  possible  fusion  yield  increase  resulting  from  di- 
muonic  hydrogen  reactions  formed  by  dt/.ix  +  dt/.ix  collisions. 


On  the  other  hand,  if  Eact =  0  and  the  di-muonic  molecules  formed  do  not  result  in 
significant  fusion  reactions,  but  rather,  bind  up  the  muons  so  they  are  no  longer  available 
to  participate  in  catalysis,  a  negative  muon  flux  of  1.5  x  105  could  decrease  the  fusion 

o 

yield  by  as  much  as  10%,  and  a  flux  of  5  x  10l  could  quench  the  fusion  rate  by  as  much 
as  95%  (see  Figure  8-2). 
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Figures  8-2a  and  8-2b.  Maximum  possible  decrease  in  fusion  yield,  as 
a  function  of  negative  muon  flux,  resulting  from  di-muonic  hydrogen 
reactions  formed  by  dtp x  +  dtp x  collisions. 


147 


If  the  reaction  rate  for  diutx  +  d/itx  is  much  lower  than  the  collision  rate  (i.e.,  Eact 
is  large)  and  doesn’t  result  in  the  formation  of  significant  quantities  of  di-muonic 
hydrogen  molecules,  then  the  di-muonic  reactions  which  could  have  the  greatest  effect  on 
the  fusion  rate  are  those  involving  x//  +  xjuxx  collisions.  Assuming  Eact  for  these 
reactions  to  be  zero  (i.e.,  all  collisions  result  in  the  formation  of  di-muonic  molecules), 
the  lowest  muon  flux  which  can  result  in  a  significant  increase  in  the  fusion  yield  can  be 
calculated.  The  lowest  flux  for  which  these  reactions  could  be  observed  will  occur  in  an 
all  deuterium  system.  For  these  reactions  to  increase  the  fusion  yield  by  10%,  a 
continuous  negative  muon  flux  of  at  least  1.5  x  10  is  needed  (see  Figure  8-3a). 
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Maximum  Increased  Yield 
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Figures  8-3a,  8-3b  and  8-3c.  Maximum  possible  fusion  yield  increase 
resulting  from  di-muonic  hydrogen  reactions  formed  by  dd/.id  +  du 
collisions. 
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At  higher  flux  an  upper-bound  on  the  fusion  increase  that  could  be  obtained  via  these  di- 
muonic  reactions  can  be  determined.  Up  to  a  30  fold  increase  in  the  d-d  fusion  rate  is 
possible  (see  Figure  8-3). 

If  the  di-muonic  molecules  d^^d  or  d^dfid  form  readily,  but  don’t  lead  to  fusion 
products,  a  quenching  effect  occurs  (see  Figure  8-4). 
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Figures  8-4a  and  8-4b.  Maximum  possible  decrease  in  fusion  yield,  as  a 
function  of  negative  muon  flux,  resulting  from  di-muonic  hydrogen 
reactions  formed  by  ddfid  +  d/.i  collisions. 


For  d-t  fusion  the  same  effects  are  predicted,  but  they  will  occur  at  a  much  higher 
muon  flux  (see  Figures  8-5  and  8-6). 
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Figures  8-5a,  8-5b  and  8-5c.  Maximum  possible  fusion  yield  increase 
resulting  from  di-muonic  hydrogen  reactions  formed  by  dt/ix  +  xju 
collisions. 


As  can  be  seen  from  Figure  8-5c,  a  three-fold  increase  in  the  fusion  yield  is  an  upper 


bound  on  what  can  be  obtained  via  this  reaction  pathway,  and  this,  only  with  an 


extremely  high  muon  flux,  which  is  currently  unobtainable. 


If  quenching  occurs,  the  effects  will  not  be  noticed  at  a  muon  flux  lower  than 

24  25 

about  10  .  A  50%  decrease  in  yield  could  not  occur  with  a  flux  lower  than  about  10“ 


(see  Figure  8-6). 
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Figures  8-6a  and  8-6b.  Maximum  possible  decrease  in  fusion  yield,  as 
a  function  of  negative  muon  flux,  resulting  from  di-muonic  hydrogen 
reactions  formed  by  dt/ix  +  x/i  collisions. 
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8.4  Conclusions 


Di-muonic  hydrogen  molecules  have  the  potential  of  affecting  the  overall  fusion 
rate  positively  or  negatively  (see  Chapter  7).  No  di-muonic  effects  will  be  observed, 
however,  unless  a  significant  quantity  of  these  molecules  are  fonned.  In  this  chapter 
lower  bounds  on  the  muon  flux  needed  to  observe  these  effects  were  presented. 

In  order  for  di-muonic  hydrogen  reactions  to  contribute  significantly  to  the  overall 
muon-catalyzed  fusion  yield  the  muon  flux  must  be  greater  than  has  traditionally  been 
used  for  muon-catalyzed  fusion  experiments.  The  lowest  continuous  muon  flux  for 
which  di-muonic  reactions  could  be  important  to  the  overall  fusion  yield  is  about  I  05 
muons/s. 

Di-muonic  reactions  can  be  important  at  this  relatively  low  flux,  only  if  xjuxx  + 
xjuxx  collisions  are  highly  likely  to  result  in  the  fonnation  of  di-muonic  molecules.  Due 
to  the  binding  energy  and  Coulombic  repulsion  between  the  xjux  portions  of  the 
molecules,  djidd  and  tjxtt  molecules  cannot  be  expected  to  fonn  di-muonic  molecules  at  a 
significant  rate  via  this  mechanism.  The  only  apparent  possibility  for  this  reaction 
mechanism  to  be  important  is  if  dpitx  +  d/utx  collisions  readily  fonn  di-muonic  molecules. 

If  d/utx  +  d/utx  reactions  do  not  have  a  significant  probability  of  occurring  (i.e.,  the 
activation  energy  is  large),  then  xjuxx  +  x/u  reactions  could  be  significant.  For  d-d  fusion 
these  reactions  cannot  be  observed  at  significant  rates  if  the  muon  flux  is  below  about 
1.5xl020  fusions/muon.  A  maximum  increased  fusion  yield  of  up  to  3,000%  is  possible. 
For  d-t  fusion,  this  reaction  path  could  result  in  10%  quenching  at  a  flux  of  1024 
muons/fusion.  If  di-muonic  hydrogen  results  in  an  increase  in  the  d-t  fusion  yield,  the 
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26 

lowest  flux  which  could  result  in  a  10%  increase  is  10“  muons/s,  with  a  maximum 
possible  yield  increase  of  300%. 
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IX.  Conclusions 


In  order  for  muon-catalyzed  fusion  to  become  an  efficient  source  of  pure  fusion 
energy  novel  reaction  mechanisms  must  be  considered.  Since  the  discovery  of  muon-catalyzed 
fusion  in  1947  considerable  research  has  been  conducted  on  single-muon  catalytic  processes  (see 
Chapter  2  and  Appendix  A).  Although  great  strides  have  been  made  towards  understanding 
muon  catalysis,  the  yields  obtained  have  been  about  an  order  of  magnitude  lower  than  what  is 
required  to  produce  a  pure  fusion  reactor.  Two  factors  have  significantly  limited  the  maximum 
obtainable  yield:  1)  the  fonnation  rate  of  molecules  which  efficiently  fuse  and  2)  the  sticking  of 
muons  to  helium  nuclei.  It  has  been  calculated  that  under  optimum  conditions  ( e.g temperature 
and  pressure)  single  muon  catalysis  could  yield  up  to  300  fusions/muon.[l:97;  2;  3]  This 
continues  to  fall  short  of  the  at  least  1000  fusions/muon  required  to  produce  a  pure  fusion 
reactor.[3;  4;  5;  6;  7;  8]  The  greatest  hope  of  using  muon  catalysis  as  a  source  of  pure  fusion 
energy  lies  in  the  discovery  of  novel  reaction  mechanisms. 

When  the  thennal  muon  flux  in  a  reaction  chamber  is  low,  the  concentration  of  di- 
muonic  hydrogen  molecules  is  insignificant;  however,  as  the  flux  increases,  di-muonic  species 
are  expected  to  play  an  increasing  role  in  affecting  the  fusion  yield.  This  document  examines  di- 
muonic  hydrogen  and  helium  molecules  and  the  possibility  that  these  exotic  compounds  can 
enhance  the  number  of  fusions  per  muon  obtainable.  In  order  for  di-muonic  molecules  to 
enhance  the  fusion  yield  the  formation  rate  of  molecules  which  efficiently  fuse  must  increase, 
and/or  these  molecules  must  facilitate  the  liberation  of  muons  stuck  to  helium  nuclei.  The  work 
presented  in  this  dissertation  examines  both  of  these  possibilities. 
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In  order  to  facilitate  the  study  of  muonic  molecules,  and  in  general,  any  exotic  molecule 
with  a  sufficiently  long  half-life  for  its  chemical  properties  to  be  of  interest,  a  General  Particle 
Orbital  (GPO)  method  of  non-adiabatic  quantum  mechanics  was  developed.  This  non-adiabatic 
multi-configurational  ab-initio  method,  which  is  outlined  in  Chapter  4,  facilitates  the  study  of 
molecular  systems  containing  any  number  of  any  type  ( i.e .,  mass  and  charge)  of  quantum 
particles  in  systems  that  may  also  contain  classical  {i.e.,  fixed)  particles.  The  size  of  the  systems 
studied,  the  number  of  quantum  particles,  and  the  types  of  quantum  particles  is  limited  only  by 
the  computational  facilities  available.  The  results  of  using  this  GPO  method  to  study  di-muonic 
molecular  properties  are  presented  in  Chapters  4  through  7. 

Dynamic  correlation  interactions  between  particles  have  a  large  impact  on  the  binding 
energy  and  bond  length  of  the  muonic  molecules  presented  in  this  document  (see  Chapter  5). 
Configuration  interaction  (Cl)  methods  of  accounting  for  these  interactions  have  been  developed 
and  presented  (see  Chapter  4).  For  di-muonic  protium  molecules  the  pp,  pp,  and  pp  correlation 
interactions  all  contribute  significantly  to  the  overall  energy  of  the  systems.  These  interactions 
account  for  as  much  as  60%  of  the  total  energy  of  muonic  protium  systems.  Of  the  three 
correlation  contributions,  pp  correlation  is  by  far  the  most  important.  This  indicates  that  the 
results  of  previously  published  calculations,  which  neglect  pp  correlation  interactions,  cannot  be 
expected  to  have  accurate  binding  energies. [9]  The  calculation  of  dynamic  correlation  energy 
using  GPO/CI  methods  is  affected  significantly  by  basis  set  size,  but  not  affected  as  significantly 
as  is  the  Hartree-Fock  Energy  (at  least  for  the  systems  studied). 

The  work  presented  in  this  document  shows  that  di-muonic  hydrogen  molecules  have 
some  unusual  properties.  Diatomic  di-muonic  hydrogen  molecules  {e.g.,  pppp)  have  a  larger 
equilibrium  p-p  bond  length  than  does  corresponding  single-muon  molecules  {e.g.,  ppp+).  This 
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effect  is  the  result  of  strong  dynamic  correlation  interactions  between  particles  (see  Chapters  4 
and  5).  Additionally,  in  di-nuclear  hydrogen  molecules,  the  vibrational  excitation  energy  of  the 
nuclei  is  higher  than  the  transition  energy  between  singlet  and  triplet  muons  ( i.e anti-parallel 
and  parallel  muon  spin).  Using  pppp  as  an  example,  it  was  shown  that  bound  singlet  and  triplet 
muon  spin  states  and  multiple  excited  rotational  states  of  this  molecule  exist  (see  Chapter  7). 
The  equilibrium  bond  length  of  the  triplet  molecules  was  calculated  to  be  almost  five  times 
larger  than  that  of  the  singlet  molecules  (see  Table  6-3).  It  was  shown  that  bound  excited 
vibrational  levels  of  pppp  do  not  exist.  Although  significant  differences  between  molecules 
with  parallel  and  anti-parallel  muon  spin  were  identified,  the  same  cannot  be  said  about  nuclear 
spin.  The  calculated  physical  properties  of  the  molecule  were  indistinguishable  for  parallel  and 
anti-parallel  nuclear  spin  states.  This  does  not  mean  that  nuclear  spin  is  not  important  to  the 
muon-catalyzed  fusion  process.  The  possibility  that  nuclear  spin  affects  physical  properties  that 
were  not  calculated  as  part  of  this  work  ( e.g fusion  rates)  exists. 

The  formation  rate  of  muonic  hydrogen  molecules  which  have  little  energy  transfer 
between  reactants  and  products  is  predicted  to  be  much  higher  than  for  those  reactions  which 
involve  the  transfer  of  large  amounts  of  energy. [10;  11]  Given  this  fact,  it  is  more  likely  that 
pppp  will  form  in  the  triplet  muonic  spin  state  than  in  the  singlet  state  (see  Figure  7-1).  The 
nuclear  bond  length  of  the  parallel  muon  spin  molecules  is  too  large  for  them  to  have  a 
significant  fusion  rate. [4]  As  a  result,  if  the  existence  of  di-muonic  hydrogen  molecules  is  to 
enhance  the  muon-catalyzed  fusion  yield,  an  efficient  mechanism  for  transitioning  triplet  di- 
muonic  molecules  to  molecules  which  do  have  a  high  fusion  rate  must  exist.  No  efficient  means 
by  which  the  triplet  molecules  being  discussed  could  transition  into  more  tightly  bound  anti¬ 
parallel  spin  states  were  found. 
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Methods  by  which  the  triplet  state  diatomic  molecules  can  transition  to  single  muon 
molecules  ( e.g pup  )  or  singlet  state  triatomic  oblate  symmetric  top  molecules  (e.g.,  ppppp  and 
pppppe)  which  have  much  shorter  nuclear  bond  lengths  (e.g.,  -0.003  A)  have  been  identified 
(see  Chapter  7).  All  of  the  possible  reaction  paths  for  these  reactions  involve  collisions  with 
other  molecules  and/or  the  formation  of  rotationally  excited  molecules.  Singlet  pppp  has  three 
bound  rotational  states  (J  =  0,  1,  2).  ppppm  is  only  bound  in  the  ground  rotational  state.  Both 
singlet  and  triplet  (muon  spin)  forms  of  ppppp  exist.  Both  of  these  forms  of  the  molecule  have 
five  bound  rotationally  excited  states:  (1,1),  (1,0),  (2,2),  (2,1),  (3,3)  (see  Chapter  7).  No  bound 
vibrationally  excited  forms  of  pppp  or  ppppp  exist.  Vibrationally  excited  forms  of  comparable 
molecules  containing  heavy  hydrogen  isotopes  may  exist,  but  this  has  so  far  not  been  calculated. 

There  are  two  di-muonic  reaction  paths  which  have  the  possibility  of  increasing  the 
muon-catalyzed  reaction  yield  sufficiently  that  muon-catalyzed  fusion  could  become  a 
standalone  source  of  energy.  The  first  possibility  is  for  dptx  +  dptx  collisions  to  result  in  the 
fonnation  of  di-muonic  hydrogen  molecules  which  quickly  react  to  fonn  fusion  products.  The 
second  possibility  is  for  dpdd  +  dp  collisions  to  result  in  di-muonic  hydrogen  molecules  which 

j 

quickly  go  to  fusion  products,  while  at  the  same  time  p  He  +  dp  collisions  result  in  fusions 
which  release  muons,  thereby  decreasing  the  effective  sticking  constant  ( i.e fraction  of  muons 
bound  to  He  nuclei). 

In  order  for  dptx  +  dptx  collisions  to  significantly  increase  the  fusion  yield,  the 
activation  energy  (Eact)  leading  to  di-muonic  hydrogen  formation  must  be  low.  If  Eact  for  this 
reaction  is  near  zero,  then  muon  fluxes  as  low  as  105  may  result  in  significant  increases  in  the 
overall  fusion  yield  (see  Figure  8-1). 
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If  Eact  for  this  reaction  is  large,  then  a  muon  flux  of  greater  than  1024  will  be 
needed  to  observe  changes  in  the  d-t  fusion  rate  due  to  di-muonic  molecular  effects  and  an 
increase  in  yield  of  300%  is  an  upper  bound  on  what  is  possible. 

For  di-muonic  reactions  to  significantly  impact  d-d  fusion  a  negative  muon  flux  greater 
than  10'  is  needed.  An  upper  bound  on  possible  yield  enhancement  due  to  di-muonic  deuterium 
molecules  was  calculated  to  be  3,000%.  Although  this  reaction  alone  is  not  likely  to  increase  the 
fusion  yield  sufficiently  for  it  to  be  used  as  a  standalone  source  of  energy,  if  coupled  with  a 
decrease  in  the  effective  He  sticking  constant,  then  there  is  a  chance  that  d-d  muon-catalyzed 
fusion  could  become  viable. 

A  method  of  using  di-muonic  molecules  to  liberate  muons  stuck  to  helium  nuclei  was 
considered.  One  of  the  products  of  d-d  and  d-p  fusion  is  He.  The  possibility  of  using  two 
muons  to  catalyze  He-  He  and  d-  He  fusion  was  studied.  Although  a  localized  energy  minimum 
for  (3 Hep  Hep)2+  exists  when  the  nuclei  separation  is  around  0.0034  A,  no  bound  vibrational 
state  of  this  ion  exists  (see  Figures  6-1  and  6-2).  Therefore,  using  two  muons  to  catalyze  He- 
He  fusion  is  not  a  viable  option.  The  results  of  using  two  muons  to  catalyze  d-  He  fusion  were  a 
bit  better.  Two  bound  vibrational  states  were  identified  with  an  equilibrium  nuclear  bond  length 
of  about  0.0037  A  (see  Figure  6-4  and  Table  6-1).  In  order  to  accurately  calculate  the  vibrational 
energy  levels  it  was  necessary  to  include  nuclear  volume  in  the  calculations.  Assuming  the 
nuclei  to  be  point  charges  results  in  a  calculated  binding  energy  of  the  J=  1  vibrational  state 
almost  an  order  of  magnitude  larger  than  when  nuclear  volume  is  considered  (see  Figure  6-3  and 
Table  6-1).  Due  to  the  existence  of  a  very  weekly  bound  vibrational  state  of  this  molecule  ( i.e 
Binding  Energy-  0.16  eV)  there  is  a  stro  ng  possibility  that  a  resonant  formation  mechanism 
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exists  which  could  enhance  the  formation  rate  of  this  molecule  (see  References  [1:82-83],  [4], 
[10],  [11],  [12],  [13],  and  [14]  for  an  example  of  a  similar  resonance  formation  mechanism). 
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Appendix  A.  Single  Muon,  Muon-Catalyzed  Fusion  Reactions  and  Rate 

Equations 

While  the  basic  idea  of  a  negative  muon  replacing  an  electron  in  a  hydrogen  atom 
bringing  the  nuclei  closer  together  and  thereby  catalyzing  fusion  may  sound  like  a  simple 
concept,  the  number  of  reaction  steps  involved  and  the  complexity  of  the  rate  equations  is 
significant.  This  appendix  outlines  the  reaction  steps  involved  in  muon-catalyzed  fusion 
and  presents  kinetic  rate  equations  which  correspond  to  these  equations.  It  may  be 
possible  to  consider  systems  that  do  not  contain  appreciable  concentrations  of  tritium, 
thereby  simplifying  these  equations  considerably.  It  is,  however,  unlikely  that  any 
system  could  contain  sufficiently  low  concentrations  of  protium  or  deuterium  so  as  to 
make  complete  neglect  of  these  isotopes  possible.  Due  to  the  difference  in  magnitude  of 
certain  reaction  rates  and  relative  concentrations  of  isotopes,  much  can  be  done  to 
simplify  the  kinetic  equations  of  specific  systems  as  the  limit  for  some  of  the  terms 
approach  zero.  No  attempt  to  simplify  the  equations  in  this  manner  has  been  made  in  this 
appendix;  all  of  the  terms  are  presented,  regardless  of  the  size  of  the  reaction  rate 
constants. 

In  a  hydrogen  only  system,  there  are  13  reactions  that  fonn  the  monatomic 
muonic  molecules  pp,  dp,  and  tp.  The  probability  of  a  specific  muonic  molecule  forming 
via  a  specific  path  is  dependent  upon  the  collision  rate,  molecular  orientation,  energy  of 
the  reacting  particles  (i.e.,  temperature  of  the  system)  and  the  height  of  the  energy  barrier 
that  must  be  overcome  in  order  for  a  reaction  to  occur.  Following  is  a  list  of  the  reactions 
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between  hydrogen  and  muons  that  fonn  monatomic  muonic  molecules:  [note:  in  the 
following  equations  H  refers  to  protium  ( i.e ! II),  and  p  refers  to  a  proton  (i.e.,  ! II 
nucleus)].  It  should  be  noted  that  many  of  the  equations  in  this  section  are  not  balanced. 
In  most  cases  the  only  products  shown  are  the  muonic-molecules.  Most  papers  on  muon- 
catalyzed  fusion  ignore  the  contribution  of  atomic  hydrogen  isotopes.  While  this 
approximation  is  appropriate  for  cryogenic  systems,  it  is  not  at  temperatures  which  are 
likely  to  produce  optimum  muon-catalyzed  fusion  yields.  In  the  equations  which  follow, 
atomic  hydrogen  has  been  included: 


H*  +  fT — : >pju  (A.l) 

H2  +  /lT  h  >p)i  (A.2) 

HD  +  /lC  — — — >  pju  (A. 3) 

HT  +  p~  — >  pp  (A.4) 

£>•+//"  >dju  (A. 5) 

HD  +  //“  — d  //  (A.  6) 
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D2  +  jli~  ^  >dju  (A.7) 

DT  +  fT  — >  d  p  (A. 8) 

T%  +  p~  ^  >tp  (A. 9) 

HT  +  jLf  A'°  >  tin  (A.  10) 

DT  +  n~ — (A.l  1) 


T2+JU~  >tjH  (A.  12) 


Once  fonned,  there  are  several  reaction  paths  the  monatomic  hydrogen  muonic 
molecules  can  follow:  1)  The  muon  can  decay,  most  often  into  an  electron,  an  electron- 
antineutrino,  and  a  muon-neutrino. [1:2]  2)  The  muon  can  be  transferred,  forming  another 
monatomic  muonic  molecule  via  the  reactions: 


pp  +  D * — 


(A.  13) 
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PH  +  T •  Al4  >tjU  (A.  14) 

pp  +  HD — » dp  (A.  15) 

pp  +  HT — » tp  (A.  16) 

pp  +  D2  hl  >dp  (A.  17) 

pp  +  DT — dp  (A.18) 

pp  +  DT — ^2— > tp  (A.  19) 

pp  +  T2  >tp  (A.20) 

dp  +  H * — Al1  >  pp  (A.21) 

dp  +  r  122  >tp  (A.22) 

dp  +  H2  ^  >pp  (A.23) 
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dp  +  HD  124  >pp  (A.24) 

dp  +  HT — — — >  p/u  (A.25) 

dp  +  HT  ^  >t/.i  (A.26) 

dju  +  DT  ^  >tjLi  (A.27) 

dju  +  T2  ^-8  >  tp  (A.28) 

tp  +  H'  — ~ ^  pp  (A.29) 

tp  +  D*  i3°  >dp  (A. 30) 

t/u  +  H7 — ^-^pp  (A. 31) 

tp  +  HD — ^+^>pp  (A. 32) 

tp  +  HD  i33  >  d p  (A. 33) 
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tju  +  D2 — ju 


(A. 34) 


tjLl  +  HT »  PIU 


(A. 35) 


and 


tju  +  DT  >dju  (A. 36) 

The  probability  of  an  exchange  reaction  in  which  the  muon  is  transferred  from  a 
heavier  isotope  to  a  lighter  isotope  of  hydrogen  is  so  small  that  it  can  be  neglected.  This 
is  primarily  due  to  the  difference  in  the  binding  energy.  3)  The  tri-atomic  muonic 
molecules  (ppp)pee,  (ppp)dee,  (pdp)pee,  (pdp)dee,  (ppp)tee,  (ptp)pee,  (pdp)tee,  (ptp)dee, 
(ptp)tee,  (dpp)pee,  (dpp)dee,  (ddp)pee,  (ddp)dee,  (dpp)tee,  (dtp)pee,  (ddp)tee,  (dtp)dee, 
(dtp)tee,  (tpp)pee,  (tpp)dee,  (tdp)pee,  (tdp)dee,  (tpp)tee,  (ttp)pee,  (tdp)tee,  (ttp)dee  and 
(ttp)tee  can  form  via  the  following  reactions: 

+  ^  >(ppju)p  (A. 37) 

pp  +  H2{ T>lj — ^  >(ppju)p  (A. 38) 

pp  +  HD — ^-^(ppp)d  (A.39) 
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pp  +  HD— 

^40  >(pdp)  p 

(A.40) 

pp  +  D2  (TT ) 

— ^— >( pdp)d 

(A.41) 

pp+D2[\l} 

— *42  >{pdp)d 

(A.42) 

pp  +  HT  - 

>{PPP)t 

(A.43) 

pp  +  HT  — 

(A.44) 

pp+DT — 

^  >{pdp)t 

(A.45) 

pp+DT — 

*46  >{ptp)d 

(A.46) 

pp+t2[  TT) 

\  -ie-+(ptM)t 

(A.47) 

pp+t2[  1M<) 

A. 

oo 

(A.48) 
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dp  +  H^  (TTj — pdp )  p 

(A.49) 

dp+H2[\ slj — >(pdp) p 

(A. 50) 

dp  +  HD — ^— >( pdp)d 

(A. 51) 

dp  +  HD — ^2— >( ddp )  p 

(A. 52) 

dp  +  D >2  (TTj — >{ddp)d 

(A. 53) 

dp  +  D^{^  — ^-^(ddp)d 

(A. 54) 

dp  +  HT  — ^-^(pd  p)t 

(A. 55) 

dp  +  HT  — ^— >( dtp )  p 

(A. 56) 

dp  +  DT  — ^J-^(ddp)t 

(A. 57) 

dp  +  DT  — >{dtp)d 

(A. 58) 
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(A. 59) 


dju+T^  TT) 

| — >(dtju)t 

d  p-\-  7^  (T4) 

— >(dtju)t 

tji  +  H2(  TT) 

— — — >(pt  n)  p 

tju  +  H2(  TT) 

— ^—>(ptju)  p 

t/j  +  HD— 

>(ptp)d 

tjLi  +  HD— 

>(dtju)  p 

tju  +  Dn  (TT) 

— ^5-->{dtp)d 

tjU  +  Dn  (TT) 

— ^—>(dtp)d 

tp  +  HT  — 

Xbl  >{ptp)t 

(A. 60) 


(A. 61) 


(A. 62) 


(A. 63) 


(A. 64) 


(A. 65) 


(A. 66) 


(A. 67) 
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tju  +  HT  — >( ttju )  p 

(A. 68) 

tp  +  DT  — 1(39  >(dtp)t 

(A. 69) 

tp  +  DT  — ^2— »( ttju)d 

(A. 70) 

tp  +  T0  |TTj — ^ — >(ttp)t 

(A. 71) 

tjU  +  T 2  — 1,2  >{ttp)t 

(A. 72) 

where  (TT)  represents  parallel  nuclear  spin  and  (Ti)  represents  anti-parallel  nuclear 

spin.  When  a  monatomic  muonic  hydrogen  atom  collides  with  a  heavier  hydrogen  atom 
the  probability  that  an  exchange  reaction  will  occur  is  greater  than  the  probability  of  a  tri- 
atomic  muonic  molecule  forming. 

There  are  several  reaction  paths  these  tri-atomic  muonic  molecules  can  follow: 
The  muon  can  decay  away,  causing  the  molecule  to  break  apart  and  the  muon-catalyzed 
cycling  chain  to  cease.  An  exchange  reaction  can  occur,  forming  a  different  tri-atomic 
muonic  molecule.  The  tri-atomic  exotic  molecules  can  dissociate,  most  often  with  the 
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muon  being  bound  to  the  more  massive  of  the  two  nuclei  it  was  closest  to.  It  can  lose 
energy  through  an  Auger  transition  such  as 

[dt/j^dee - >(dt/j)de  +  e~ 

3 ) 

or,  nuclear  fusion  can  occur,  giving  off  energy  and  a  variety  of  particles.  After  Auger 
stabilization  the  molecule  can  decay,  fuse  or  collide  with  another  molecule  and 
dissociate. 

Before  presenting  the  reaction  paths  the  tri-atomic  muonic  molecules  can  follow 
there  is  another  term  which  needs  to  be  defined,  the  sticking  constant  (co,).  This  is  the 
probability  that  when  muon-catalyzed-fusion  occurs  the  muon  will  stick  to  a  fusion 
product. 

In  most  cases  muon  sticking  is  a  non-desirable  effect  since  it  precludes  the  muon 
from  catalyzing  further  fusion  reactions.  A  notable  exception  to  this  is  when  ddpp,  ddjud 
or  ddpt  fuse  forming  tritium  and  a  proton.  In  these  cases  if  the  muon  sticks  to  the 
products  they  can  immediately  enter  back  into  the  reaction  sequences  listed  in  this 
section.  Although  He  can  fuse  via  muon-catalyzed  fusion,  the  probability  is  small  due  to 
a  relatively  large  separation  distance  between  the  fusing  particles.  As  a  result,  He  fusion 
has  been  neglected  in  the  reactions  presented  in  this  section. 

Following  are  exchange  reactions  tri-atomic  muonic  molecules  can  participate  in, 
in  a  hydrogen  only  system 
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{ppm)p+d*— 

273  >(ppp)d  +H * 

(A. 74) 

(ppju)p+r~ 

^  >(ppp)t  +  H* 

(A. 75) 

( ppju)p  +  HD - 

^ 5  >{ppp)d  +H2 

(A. 76) 

(ppp)p  +  HT- 

Alb  >(ppp)t  +  H2 

(A. 77) 

(ppp)p  +  D1{\'[) 

| — - in—>(ppju)d  +  HD 

(A. 78) 

(. PPP)P  +  D2 ('N') 

— ^  >(ppju)d  +  HD 

(A. 79) 

(. ppm)p+dt - 

-^z2— >(/?/?, p)d  +  HT 

(A. 80) 

(. ppm)p+dt - 

-^-^■(pp p)t  +  HD 

(A. 81) 

{ppp)p+t2^) 

— ^-^{pp p)t  +  HT 

(A. 82) 

(. PPP)P+t2[ 

— ^3 -^y(pp  p)t  +  HT 

(A. 83) 
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( ppju)d  +  H * — ppju )  p  +  D * 

(A. 84) 

[ppp)d  +  Tm — >Q opju)t  +  D' 

(A. 85) 

(ppp)d  +  H2  (TTj — ^-^>(ppp) p  +  HD 

(A. 86) 

(ppju)d  +H2  (t -ij — ^  >{ppp) p  +  HD 

(A. 87) 

( ppp)d  +  HD — >[ppp) p  +  D2 

(A. 88) 

( ppp)d  +  HT  — -^— >(/?/?//)  p  +  DT 

(A. 89) 

(ppju)d  +  HT  — ^-^{pp  p)t  +  HD 

(A. 90) 

(ppp)d  +  DT  — ^-^>(ppp)t  +  D1 

(A. 91) 

(ppp)d  +  ro  (TTj — >(pp  p)t  +  DT 

(A. 92) 
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(ppju)d  +  T2  — ^SL^>(ppp)t  +  DT 

(A. 93) 

(ppp)t  +  H* — ^— >(/?/?//)  p  +  T * 

(A. 94) 

( ppp)t  +  D * — ^-^(ppju)d  +  T' 

(A. 95) 

(ppju)t  +  H2  (TTj — ^-^(ppju)  p  +  HT 

(A. 96) 

( ppp)t  +  H1  — ^-^(ppp)  p  +  HT 

(A. 97) 

( ppp)t  +  HD — Xc>1  >(ppp) p  +  DT 

(A. 98) 

[ppp)t  +  HD — ^— >( ppp)d  +  HT 

(A. 99) 

( ppp)t  +  D >2  |TTj — >(ppju)d  +  DT 

(A.  100) 

(ppp)t  +  D1  |T  — ^m-^>(ppp)d  +  DT 

(A. 101) 

(ppp)t  +  HT  — ^01  >{ppp)  p  +  T 

(A.  102) 
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(pp/u)t  +  DT  - 

^°2  >{ppp)d  +  T1 

(A.  103) 

( pdp )  p  +  D'  — 

^°3  >{pdp)d  +H* 

(A.  104) 

( pdp )  p+r- 

pdp)t  +  H * 

(A. 105) 

( pdp )  p  +  HD- 

Al05  >{pdp)d  +  H2 

(A.  106) 

( pdp )  p  +  HT  - 

Ai06  >{pdp)t  +  Hn 

(A.  107) 

(pdju)p+D2[  TT) 

l — ^mL~^{pdp)d  +  HD 

(A.  108) 

(pdju)p+D2(  INI.) 

l — i|08  >{pdp)d  +  HD 

(A.  109) 

[pdp]  p  +  DT — 

Al09  >{pdp)d  +  HT 

(A.  110) 

{pdp)  p  +  DT  - 

;l1I()  >{pdp)t  +  HD 

(A.  Ill) 

181 


(pdp)  p  +  T 2  (TTj — ^ 11  >(pdju)t +HT 

(A.  112) 

(pdju) p  +  T0  — ^—^(pd p)t  +  HT 

(A.  113) 

(pdp)d  +  H' — ^113  >{pdp) p  +  D * 

(A.  11 4) 

[pdp)d  +  T* — ^— >( pdju)t  +  D * 

(A.  11 5) 

(pd ju)d  +  (TT) — ^15  >(pdju) p  +  HD 

(A.  11 6) 

(pdp)d  +  |T — ^-^(pd ju) p  +  HD 

(A.117) 

(pdp)d  +  HD — A‘17  >(pd p) p  +  D2 

(A.  118) 

(pd ju)d  +  HT  — il18  >(pdju) p  +  DT 

(A.  11 9) 

[pdp)d  +  HT  — >( pdp)t  +  HD 

(A.  120) 

(pdp)d  +  DT  — ^2S~^{pdp)t  +  D2 

(A. 121) 
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(pdp)d  +  T2^  TTj — ^—^(pd  p)t  +  DT  (A.122) 

(pdp)d  +  T^  |Tsl) — ^~^{pdp)t  +  DT  (A. 123) 

(pdp)t  +  H* — ^2^(pdp)p  +  T'  (A.  124) 

( pdp)t  +  D * — ^m-^>(pdp)d  +  T'  (A.  125) 

[pdp)t  +  H1  |TTj — ^-±{pd  p)  p  +  HT  (A.  126) 

[pdp)t  +  H1{\\-  j — ^m~^>{pdp)p  +  HT  (A.  127) 

(pdp)t  +  HD — ^-^(pd p)  p  +  DT  (A.  128) 

( pdp)t  +  HD — ■ >{pdp)d  +  HT  (A. 129) 

[pdp)t  +  A  (tTj — ^m^>(pdp)d  +  DT  (A.  130) 
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(pd/j)t  +  D2  — ^m-*(pdp)d  +  DT  (A. 131) 

( pdju)t  +  HT — ^m-^(pdp)p  +  T1  (A.132) 

( pdp)t  +  DT — ;-132  >( pdp)d  +  T2  (A. 133) 

( ptp)p  +  D * — A‘33  >{ptp)d  +  H *  (A.  134) 

( ptp)p  +  T * — il34  >{ptp)t  +  H *  (A. 135) 

( ptju)p  +  HD — ^m-^[ptp)d  +  H1  (A.136) 

( ptp)p  +  HT — ^m~^>(ptp)t  +  H1  (A.  137) 

(ptju)p  +  D0  (TTj — ^m—>(ptp)d  +  HD  (A. 138) 

(/>?//) p  +  D2 1' T^j — ^m-^(ptp)d  +  HD  (A.  139) 

( ptju)p  +  DT — ^m-^{ptp)d  +  HT  (A.  140) 
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( ptju)p+DT  - 

ptp)t  +  HD 

(A. 141) 

{ptfi)p+T2(  TT) 

^  >(ptju)t  +  HT 

(A.  142) 

(ptp)P  +  T2(  tl) 

A42_>( ptju)t+HT 

(A.  143) 

(ptju)d  +  H*  — 

il43  >{ptp) p  +  D * 

(A.  144) 

(ptfi)d  +  T*  — 

-^^{pt  ju)t  +  D* 

(A.  145) 

(ptju)d  +  H2  |TTj 

— - im—>(ptp)p  +  HD 

(A.  146) 

(ptju)d  +  H2  (T 4'j 

— ^46  >(ptju)p  +  HD 

(A.  147) 

( ptp)d  +  HD - 

Al47  >{ptp)  p  +  D1 

(A.  148) 

( pt/u)d  +  HT  — 

^14S  >(ptju) p  +  DT 

(A.  149) 
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(ptju)d  +  HT  - 

ptju)t  +  HD 

(A.  150) 

( ptju)d  +  DT  - 

^ 50  >(ptjLl)t+D2 

(A. 151) 

(ptp)d  +t2  |TT  | 

| — >(ptp)t  +  DT 

(A.  152) 

(ptp)d  +T2 

| — ^2— >(ptju)t  +  DT 

(A.  153) 

(pt p)t  +  H*  — 

il53  >(ptp)p  +  r 

(A.  154) 

(ptju)t  +  D'  — 

ptp)d+r 

(A. 155) 

(ptju)t  +  H2  |TT  j 

— ^~^(ptp)p  +  HT 

(A.  156) 

(ptju)t  +  H2 

— ^-^(pt p)  p  +  HT 

(A.  157) 

( ptp)t  +  HD — 

^157  >(ptju) p  +  DT 

(A.  158) 

(ptju)t  +  HD— 

A‘58  >{ptju)d  +  HT 

(A.  159) 
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{ptp)t  +  D0  (TTj — ?m—>{ptp)d  +  DT 

(A.  160) 

(ptp)t  +  D1  — ^m-^>{ptp)d  +  DT 

(A. 161) 

( ptp)t  +  HT  — 2)61  >(ptp)  p  +  T-i 

(A.  162) 

( ptp)t  +  DT  — ^62  >{ptp)d  +  T2 

(A.  163) 

{ddp) p  +  D * — /L|63  >{dd p)d  +  H * 

(A.  164) 

(i ddp)p  +  T * — il64  >(ddp)t  +  H * 

(A.  165) 

{ddp)  p  +  HD — ^-^{dd  p)d  +  H1 

(A.  166) 

{ddp) p  +  HT  — il66  >{ddp)t  +  H1 

(A.  167) 

{ddp)  p+D2i  TTj — 1167  —>{ddp)d  +  HD 

(A.  168) 
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( ddju)p  +  D2  |T>lj — ^-^{dd  /j)d  +  HD 


(A.  169) 


(ddju)p  +  DT — 

-^>( ddp)d  +  HT 

(A.  170) 

(ddju)  p  +  DT — 

il7°  >(ddju)t  +  HD 

(A. 171) 

(ddju)  p  +  T 2  (TTj 

— >(dd  ju)t  +  HT 

(A.  172) 

(. ddju)p  +  T2 

A72— >  (ddju)t  +  HT 

(A.  173) 

(ddju)d +H*  — 

>(ddju) p  +  D% 

(A.  174) 

(ddju)d  +  r- 

>(ddju)t  +  D * 

(A.  175) 

(ddju)d  +  H2  |TTj 

— ^  ->(ddju) p  +  HD 

(A.  176) 

{ddp)d  +  H2 

— A'76  >(dd/u)  p  +  HD 

(A.  177) 

(ddju)d  +  HD- 

Al77  »( ddju )  p  +  D >2 

(A.  178) 
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(A.  179) 


( ddju)d  +  HT  — il78  >(ddju) p  +  DT 
(< ddju)d  +  HT — ^-^(dd  p)t  +  HD  (A.  180) 

(< ddp)d  +  DT — ^  >{ddju)t  +  D1  (A.181) 

(ddju)d  +  T2^  TTj — ^m-^>(ddp)t  +  DT  (A.182) 

(dd]u)d +  T2[\\ l-j — ^-^{dd  ju)t  +  DT  (A.183) 

( ddp)t+H* — ^183  >(ddju)d  +  T *  (A.  184) 

(< ddju)t  +  D' — ^m-^(ddp)p  +  T*  (A.  185) 

(ddju)t  +  H2  |TTj — ^m^{ddp)p  +  HT  (A. 186) 

(ddp)t  +  H2^l-l ) — ^86  > (dd p)p  +  HT  (A. 187) 
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(i ddju)t  +  HD 


^187 


*{ddp)  p  +  DT 


(A.  188) 


{dd  p)t  +  HD — ^88  >{ddp)d  +  HT 

(A.  189) 

{dd  p)t  +  |TTj — ^S9  >{ddp)d+DT 

(A.  190) 

{ddp)t  +  D2  — ^m-^{ddp)d  +  DT 

(A. 191) 

(dd  p)t  +  HT  — Al91  >(dd  ju)p  +  T1 

(A.  192) 

(ddju)t  +  DT  — 1,92  >(dd p)d  +  T2 

(A.  193) 

{dtp) p  +  D * — /1‘93  »( dtp)d  +  H * 

(A.  194) 

{dtp) p  +  T ’ — Al94  >{dtp)t  +  H * 

(A.  195) 

{dtp)p  +  HD — ^9S  >{dtp)d  +  H2 

(A.  196) 

{dtp) p  +  HT  — 2,96  >{dtp)t  +  H1 

(A.  197) 
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(A.  198) 


(< dtju)p  +  D >2  |TTj — Al97  - >{dtju)d  +  HD 

(dtju)p  +  D2( T^j — ^98  >(dtju)d  +  HD  (A.199) 

( dtfi)p  +  DT — 1199  >(dtp)d  +  HT  (A.200) 

(dtju)p  +  DT — 1200  >(dtp)t  +  HD  (A.201) 

(dtju)p  +  T2^  TTj — — — >(dtp)t  +  HT  (A.202) 

[dtp)p  +  T1  — ^~^(dtp)t  +  HT  (A.203) 

(< dtp)d  +  H * — ^°3  >(dtju)p  +  D *  (A.204) 

(^//)c/  +  r* — ^  >(dtp)t  +  D'  (A.205) 

(<^//)t/  +  //2|TTj — ^-^{dt  ju)  p  +  HD  (A.206) 

(fifr//)t/  +  //2|TNlj — dtp)p  +  HD  (A.207) 


191 


(A.208) 


(dtju)d  +  HD — ^°7  >{dlju)  p  +  D1 
( dtju)d  +  HT — ^os  >(dtju)p  +  DT  (A.209) 

(dtp)d  +  HT — ^m-^(dtp)t  +  HD  (A.210) 

( dtp)d  +  DT — '*210  >{dtfj)t  +  D2  (A.211) 

(dtju)d  +  T2^  TTj — ^m-^>(dtp)t  +  DT  (A.212) 

(dtp)d +  T2^l-l  j — ^m^{dtp)t  +  DT  (A.213) 

( dtp)t  +  H * — ;-13  >{dtp)d  +  T *  (A.214) 

(dtp)t  +  D* — ^  >(dtju)p  +  r  (A.215) 

(dtju)t  +  H0  |TTj — ^-^(dt p) p  +  HT  (A.216) 
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(dtp)t  +  H2  |TnL| — Xll6->(dtp)  p  +  HT 

(A.217) 

(< dtju)t  +  HD — 2217  >(dtp)p  +  DT 

(A.218) 

(< dtp)t  +  HD — ^m^(dtp)d  +  HT 

(A.219) 

(dtp)t  +  D1  |TTj — ;"2'9  ->  (dtp)  d  +  DT 

(A.220) 

(i dtp)t  +  D >2  — 2/1 20->(dtp)d  +  DT 

(A.221) 

(dt  p)t  +  HT  — ^21  >(dtp)  p  + 

(A. 222) 

(dtp)t  +  DT  — — — >(dtp)d  +  T2 

(A.223) 

(l ftp ) p  +  D * — ^23  >(ttp)t  +  H * 

(A. 224) 

(ttp) p  +  T * — 2/224  >(ttp)d  +  //* 

(A.225) 

(ttp)p  +  HD — ^25  >(ttp)d  +  Hn 

(A.226) 
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(A.227) 


( ttju ) p  +  HT  — ^26  >{ttp)t  +  H2 
(ttju)p  +  D2^  TTj — - m—>(ttp)d  +  HD  (A.228) 

(ttju)p  +  D2  (TTj — ^->(ttp)d  +  HD  (A.229) 

(. ttjLl)p  +  DT — >( ttfj)d  +  HT  (A. 230) 

( ttju)p  +  DT — ^  >{ttp)t+  HD  (A.231) 

{ttp)p  +  T2i TTj — ^  >{ttp)t  +  HT  (A.232) 

{ttp)p  +  T2  (TTj — ^32  >(ttju)t  +  HT  (A.233) 

( ttp)d  +  H * — 1233  >{ttp)p  +  D *  (A.234) 

{ttp)d  +  r  1234  >{ttp)t  +  D'  (A.235) 
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( ttju)d  +  H 2  (TT) — ^  ->(ttju)  p  +  HD 

(A. 23 6) 

(ttju)d  +  H2  — 2236  >(ttju)d  +  HD 

(A.237) 

(ttju)d  +  HD — 2/17,1  >(ttjLi)  p  +  D1 

(A.238) 

(ttju)d  +  HT  — ^m-^(ttp)p  +  DT 

(A.239) 

(i ttju)d  +  HT — 2239  >(ttjU)t  +  HD 

(A. 240) 

(ttju)d  +  DT  — 2240  >{tt p)t  +  D2 

(A.241) 

(ttju)d  +  r2(ttj  *(ttju)t  +  DT 

(A.242) 

( ttfi)d  +T2  — 2/142  ->{tt  p)t  +  DT 

(A.243) 

(ttju)t+H* — 2243  >{ttp)p+r 

(A. 244) 

(ttju)t+D* — 2244  >{ttp)d  +  Tm 

(A.245) 
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(ttju)t  +  H2  (TT) — ^45  >(ttju) p  +  HT 

(A. 246) 

(ttp)t  +  H1  |T — ^46  >(ttp)  p  +  HT 

(A. 247) 

(ttp)t  +  HD — ^47  >{ttp)  p  +  DT 

(A.248) 

( ttju)t  +  HD — ;-48  >(ttju)d  +  HT 

(A. 249) 

(ttju)t  +  D2  (TTj — ^49  >(ttp)d  +  DT 

(A.250) 

(, ttp)t+D2  >(ttju)d+DT 

(A.251) 

( ttp)t  +  HT  — ^  >{ttp)  p  +  T2 

(A.252) 

( ttju)t  +  DT  — ^ 5 2  >{ttju)d  +  D 


(A.253) 


Regeneration  reactions  which  “free”  muons  bound  to  helium  nuclei  will  not  be 


considered  directly  in  the  rate  equations  presented  in  this  appendix 


H 3  He - > 3 He  +  /u 


(A.254) 


/uAHe - >  AHe  +  /u 


(A.255) 


In  actuality,  regeneration  reactions  are  not  as  simple  as  is  shown  in  equations  (A.254)  and 
(A.255),  rather  they  involve  collisions  with  other  particles  and  most  often  involve  the 
transfer  of  a  muon  rather  than  actually  “freeing”  it.  These  reactions  can  be  accounted  for 
through  the  use  of  “effective”  reactions  rates.  This  means  that  reaction  paths  such  as 


(ddju)d - >jU3He  +  n - >3He  +  /j  +n 


(A.256) 


are  considered  to  be  the  same  as 


(i ddju)d - >3He  +  jU  +n 


(A.257) 


with  the  rate  constants  being  adjusted  to  account  for  both  reaction  paths.  Rather  than 
expressing  a  reaction  rate  constant  for  each  reaction  path  leading  to  fusion,  it  is  more 
convenient  to  express  the  individual  reaction  rates  as  the  total  fusion  rate  of  a  given 
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reactant  times  the  probability  that  a  specific  path  will  be  taken.  The  following  reactions 
are  expressed  in  this  manner: 


( pd ju)  p — i2S3  >  Products 


( pd ju)d  — ^54  >  Products 


(pdju)t  — ^55  >  Products 


(ptju) p — 1256  >  Products 


(J ptju)d — Al51  >  Products 


(ptju)t — Al5S  >  Products 


(i ddju)p  — /'2"'>  >  Products 


(chip) d  — Al60  > Products 


(ddjLi)t  — ^61  >  Products 


(A.258) 


(A.259) 


(A.260) 


(A.261) 


(A.262) 


(A.263) 


(A. 264) 


(A.265) 


(A.266) 
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(dtjn) p — Al62  >  Products 


(A.267) 


(dtjLi)d  — 1263  > Products 


( dtju)t  — /"64  >  Produ cts 


( tip)  p  — ^6S  >  Products 


(ttju)d  — /"66  >  Products 


(A.268) 


(A.269) 


(A.270) 


(A.271) 


and 


(#//)£  — ^67  >  Produ  cts 


(A.272) 


Defining  effective  sticking  probability  (&> )  to  be  the  probability  that  a  muon  is 
“stuck”  to  a  nucleus  during  fusion  and  remains  “stuck”  until  it  decays, 


( ox  =The  effective  probability  of  p3He  sticking  during  ppd  fusion 
co2  =The  effective  probability  of  p4He  sticking  during  ppt  fusion 
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oo3  =The  effective  probability  of  t/j  sticking  during  dpd  fusion 
co4  =The  effective  probability  of  ju^He  sticking  during  d/ud  fusion 
m5  =  The  effective  probability  of  pA  He  sticking  during  dpt  fusion 
co6  =The  effective  probability  of  pAHe  sticking  during  tpt  fusion 

From  this  it  is  possible  to  express  the  fusion  rate  constants  in  terms  of  their  effective 
sticking  probability.  The  sticking  probability  is  not  a  constant;  rather  it  is  a  function  of 
the  temperature  and  pressure  in  the  reaction  chamber.  How  much  these  values  change 
with  temperature  and  pressure  is  an  ongoing  area  of  research. [2;  3;  4] 

Nuclear  fusion  occurs  via  the  following  reactions  (note:  although  electrons  are  not 
always  shown  in  the  reaction  equations  presented  in  this  section,  gammas  and  electron- 
positron  pairs  which  have  reasonably  high  formation  rates  are). 


(pdp)  p  ^530-^Xi-n)  >JU-  +  3 He  +  Hm  (A.273) 

(pdju)p — ^222^— > y  +  p^He  +  H *  (A.274) 

(pdjLi)p — ^53 ( )/l  >/  +  //+  3 He  +  H *  (A.275) 

(pdju)d-  >  M-  +  3He  +  D*  (A.276) 
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(pdju)d  >y  +  fi3He+D'  (A.277) 

(pdju)d  ^4(1~^)Xl  > y  +  p+ 3 He + D*  (A.278) 

( pdju)  t  jU~  + 3  He +  r  (A.279) 

( pd/u)t  — ^5sWl  > y  +  p 3 He  +  T*  (A.280) 

(pdju)t — — s-5  (  *7  A  -);/|  — >  /  +  p  +  3  He  +  r*  (A.281) 

(/??//) p  ^a-^Xi-^Xi— A)  >  ^  +  *He  +  H'  (A.282) 

(ptp)p — i256^2 (1~Pl }  >y  +  p 4 He  +  H'  (A.283) 

(ptp)p — /t256(1~^)^(1~pi)  >y  +  jLl  +  4He  +  H*  (A.284) 

(ptp)p — i256^2A  >  e+g~  +  p 4  He  +  //*  (A.285) 

( ptp ) d  ^a-^xi-AXi -A)  >  +  4 He +  D'  (A.286) 
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(ptjLi)d — hsi^hBl  >  y  +  pAHe+D'  (A.287) 

(ptju)d — »^  +  p+4He+D*  (A.288) 

(ptp)d — hsnPifh  >  e+e~  +  n  4  //g  +  £)*  (A.289) 

^5s(i-^Xi-r2)(i-A)  > /r  +  4//e  +  r  (A.290) 

(ptju)t — yy  +  p 4He  +  r*  (A.291) 

(/>*//)  t  hsdt^rdtPd  >  r  +  jU+4He  +  r  (A.292) 

(ptp)t — ^m^-^e+e~  +  p4He  +  T*  (A.293) 

(ddp)p — 'W1-^  )(!-£„)  >  t  +  p  +  //*  +  p~  (A.294) 

(i ddp)p — >tju  +  p  +  H'  (A.295) 

(< ddju)p  > 3 He  +  H'+  p  +  //“  (A.296) 
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[eld ft)  p  h^Bn  >  ju 3 He  +  Hm  +  nx  (A.297) 

(ddfi)d — >t  +  p  +  £)»+jU-  (A.298) 

( ddju)d — 17606)3 (1~B" )  >tju  +  p  +  D'  (A.299) 

(. dd ft) d  ^od-^4)A  >  3He  +  D'+  nx  +  fi  (A.300) 

(ddju)d  A76°c°4Bn  > // 3 He  +  £>*  +nl  (A.301) 

(ddju)t  ^id-^Xi ~Bn)  >t  +  p  +  r  +  jU~  (A. 302) 

( ddfl)t  ^O-A,)  >tjU  +  p  +  r  (A. 303) 

{ddju)t  }™(x~0)dBn  > 3 Re  +  r  +  ni+ju-  (A.304) 

(ddju)t  >n3He+r  +  nx  (A.305) 

(dtju)p — 1762  (1~Ws }  >  4  He  +  H'  +  n2  +  p~  (A.306) 


203 


(dtp)  p 


(A. 307) 


>p4He  +  H'  +  n2 

(dtju)d  ^3(^5)  >4He+Dm+n2+p-  (A.308) 

(dlp)d  >  p4  He +  D'  +  n2  (A.309) 

(dtp)t  ^4(!-a)  > 4 He  +  T‘  +n2+  p~  (A.310) 

(dtp)t — A°-mC°5  >  p 4 He  +  T *  +  n2  (A.311) 

(Up) p  >  *He  +  //•  +  2n2  +  p~  (A. 3 12) 

(»//) p  >  p4  He +Hm  +  2 n2  (A.313) 

(ttju)d- *  He  +  D'+2n2+p~  (A.314) 

(ttju)d  A2b6(°6  > p4 He  +  D’  +  2n2  (A.315) 

(, ttp)t  ^7(1^6)  >4He  +  r  +  2n3+p-  (A.316) 
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and 


(ttjLi)t — juAHe  +  T*  +2n^  (A.317) 

Reactions  A. 273  through  A.317  can  be  used  to  detennine  the  fusion  rate  of 
specific  paths.  yi  represents  the  probability  of  gamma  fonnation.  pi  is  the  probability  of 

an  electron  positron  pair  ( e  e+  )  forming,  d/id  fusion  can  follow  a  path  which  produces 
tritium  and  a  proton,  or  which  produces  He  and  a  neutron.  The  branching  ration  for 
these  two  paths  is  the  same  as  it  is  for  conventional  fusion.  The  probability  ( B '„)  of  the 

3 

He  and  neutron  path  being  followed  is  about  58%. [5]  The  rate  constants  shown  in  the 
following  equations  assume  the  non-reacting  hydrogen  to  have  negligible  impact  on  the 
sticking  probability,  the  probability  of  electron  positron  pair  ( e  e'  )  production  and  the 
probability  of  gamma  production. 

From  the  reaction  equations  presented  in  this  appendix,  it  is  possible  to  derive 
kinetic  rate  equations  for  the  various  products  and  intermediates: 

=  2,  [»•][>]  +  2,  [//,][//-]  +  2,  [HD][>- ]  +  2,  [OT][>]  -  \  [/>//] 

-  A3  [PfA  [D*]  “  A4  [P/A  [r  ]  -  As  [pp]  [HD]  -  ^16  [pp\ lHT ] 

-  \i  [pp]  [A  ]  -  As  [pp]  [DT ]  -  A 9  [pp]  [DT ]  -  Ao  [pp]  [Ti  ] 

-  Ai7  [pp]  [h2  (■ n )]  -  As  [pp]  [h2  (U)]  -  A,  [pp]  [HD]  (A.3 1 8) 

-  Ao  [pp]  [hd]  -  Ai  [pp]  [ d2  )]  -  A2  [pp]  [n 

-  Ab  [pp]  [HT]  -  A4  [pp]  [HT]  -  As  [pp]  [DT]  -  A6  [pp]  [DT] 

-  At  [pp]  [a  )]  -  As  [pp]  [Ti  (A i) 
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^4  =  ^[D-\^]+K[HD}[^}+A1[D1}[^}  +  \[DT}[^y^[d^ 

-  4  [d»][H'\- A„  [df,][r\-  A,.  [<///][«.]- 4,  [rf//][«D] 

- 125  [  J//]  [//r]  -  a26  [  J//]  [//r]  - 127  [dp]  [ dt ]  -  a28  [J//]  [r2  ] 

-  A49  [dp]  [tf2  (TT)]  -  250  [J//]  [i/2  (tl)]  - 25,  [<///]  [iZD]  (A.3 19) 

-  A52  [ dju\  [HD]  - A53  [J//] \d2  (TT)]  -  A54  [dp] \d2  (Tl)’ 

- 155  [dp]  [HT]  -  A56  [dp]  [ HT ]  -  A57  [dp]  [DT]  -  258  [dp]  [DT] 

- ^  [r  ][>"  ] + A,  [ht\[h-  ] + a, ,  [DT\[n-  ] + 4  [r2  ]  [V  ]  -  \  [dfi] 

i:]\j r \  [///][«_] 

- 4:  [</'][/m] -  a.,  [»/<][ A ] -  a„  [/,//][//r ] -  a„  [</i][fjr ] 

-l„[(//][//2(tt)]-l62[f//][//2(1'i)]-4[<//][HD]-iM[<//][HD]  (A. 320) 

-KA<A[DAn)]-xMtA[DAn)]-xM<A[HT]-^Vt‘IHT] 

-  4  [/./,]  [/_>/■]  -  [///]|  /->/■]  -2,,  [///]  [r2  (TT)]  -  A  [<//][r2  (U)‘ 

+  /l275  [(JJ//) /?]  +  /l279  [[ddju)d~^  +  / 1283  [( ddp)t^ 

\{ppm)p\  =  ^  [^2  (TT )]  +  238  [/>//]  [//2  (1M]]  -  A(l  \{ppfu)p] 

-  ^73  [(  ppm)  p]  [d’  ]  -  ^4  [(  ppm)  p]  [r  ]  -  ^75  [( PPM)  p]  [HD] 

~  K  [( PPM) p]  [HT ]  -  277  \[ppp)  p]  [D2  (TT )]  (A. 321) 

-K\{pPM)p\d2^1))[-{A19+\0)\(ppp)p][DT] 

“4si  \{pPM)p\t2  (tt)]-^2  \{ppp)p^T2  (Tl)] 
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d[{pPM)d] 

dt 


g[C ppm)*] 

dt 


d\[pd/u)p~] 

dt 


^9  [pm] [hd]  ~  K  [( ppp)d  ]  -  43  [(  ppm)  d][H'] 

-K\{ppM)d\r^-\5\{pp^)d^LH  2  (TT )] 
“4,  \_(pPP)d\H 2  (tl)]-^7  [(ppju)d][HD] 
-{\*+\v)[{pPP)d][HT]-K\(PPP)d][DT] 
“4  [_(PPM)d\T2  (tT)]-^  \[PPP)d\T2  (Tl)] 


4  M  ]  -  K  \_{PPP)  *]  ~  4  \_{PPP)  *]  [ H ’  ] 

-  4  [( PPM)*]  [O'  ]  -  4  [( PPM)*]  [H2  (tt )] 

“4  \_(pPP)t1^H2  (4)J-(4  +  A9S)[(ppp)t][HD] 
“4  [(ppfi)t][D2  (TT )]  -  Aoo  [(pPM)*][D2  (^)] 
~\m\{pPM)t][HTY\o2\{PPM)*][DT] 


4  [pm]  [HD]  +  4  [dp]  [H2  ) 


+  ^[dp][H2(U) 


-\[{pdp)p]-\^pdp)p][D’]-\QA\{pdp)p]\T’] 

~  4s  \[PdP)p]  [HD]  ~  4e  \{PdM) p]  [HT ] 

-  ^ 07  [( pdp)p]  [A  (ft )]  -  \08  [( pdp)p]  [d2  (t 

-^m+\m)\[PdP)p][DT]-^u\[PdP)p]\T2^) 
\_{PdM)p\T2  (ti)]-4s3  \{pdp)p] 


(A. 3 22) 


(A. 323) 


(A. 324) 
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d[{pdp)  d] 
dt 


d[{pdp)t] 

dt 


d[{ptp)p~] 

dt 


;  A  [/>//]  [a  (TT)]  +  A  [/V*]  [a  (T 'l')]  +  4 1  [dp]  [HD] 


~  ^0[{pd  p)d~]~ \n[{pd  p)d~][H’^~ \u  [{pd p)d~][r'\ 

-  A5  [(aAA][A  (tt)]-^16  (TT)] 

-^i7[(^//)j][JHD]-(ylii8+yl1i9)[(jpJ//)j][/fr] 

-  Ao  [( pdp)  d]  [DT ]  - \2l  [( pdp)d][r2  (TT )] 

~A2  [{Pdp)d1^T2fti)^-X2^[{pdp)d'] 


A  [pp] [DT ]  +  A  [dp] [HT ]  -  A0  [( pdp) t]  - [( pdp) t] [H’ ] 
-  A4  [( pdp)t\ [d’ ]  - \25  [( pdp) t] \jl2  (TT )] 

-\i6{{pdp)t\H2^l)^-(\21+\2%)\(pdp)t][HD] 

~  A9  \{pdp)t]  [A  (TT  )]  -  Ao  \[pdp)t\ [A  (T  ^)] 

-\n[{pdp)t~][HT]-\n[{pdp)t~][DT]-X25fi\{pdp)t~] 


AAA[pp][HT]  +  A6l[tp][H2ft  T) 


+A  [tp][H2ftlj 


~  A  {{Ptp)  P]  -  A3  [{PIP)  P]  ]  -  A4  [{PIP)  P]  [T’  ] 

-  As  [{ptp)  p]  [HD]  -  ^36  [{ptp)  p]  [HT ] 

-  \37  [{ptp)  p]  [d2  (TT )]  -  a8  [( ptp )  p]  [A  (T  T) 

"  ( A9  +  Ao  ) \{ptp)  p]  [DT ]  -  Al  [{ptp)  p]  [T2  (TT )_ 

-  A2  [{ptp)  p]  [t2  (A)]  -  A«  [{ptp)  p] 


(A. 325) 


(A. 326) 


(A. 327) 
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d\_{ptfu)d] 

dt 


d\_(pt/u)t~] 

dt 


d[(ddp)  j p] 
dt 


A6  [pm\  [DT ]  +  A63  [tju]  [HD]  -  A0  [( ptju)  d ]  - 


A43  [( ptp)d]  [_Hm  ]  -  Al44  [( ptp)d]  [r  ] 

-  A4S  [( ptp)  d]  [h2  (TT)]  -  A46  \[ptp)  d]  [h2  (Tl)] 

-  ^47  [( ptp)d] [HD]  ~(A4S  +  A49 )\_{ptp)d] [HT ] 
-AX5Q\{ptp)d][DT]-AXiX\{ptp)d]\T2[^)^ 

~  Asi  [( Ptp) d]  \_T2  ( tl)]  -  A?  [( ptp) d] 


A47  [ pp ]  [t2  (■ n )]  +48  [pp]  [t2  (U)]  + A  [tp]  [HT] 


-  A  [( ptp)t]  -  A53  [( ptp)t]  [h'  ]  -  A54  [( ptp)t]  [a  ] 
-Ass  [(ptp)t][H2  (TT)]- A6  [(ptp)t][H2  (TT)_ 

-  (^157  +  Ass  )  [(  ptp)  t]  [HD]  -  Al59  [(  ptp)t]  [a  (TT ) 

-  Aeo  [( ptp)t]  [■ £>2  (TT)]  -  Aei  [( ptp)t]  [HT ] 

-  A,62  [( ptp)t] [dt]-  A25S  [( ptp)t] 


:  Ai  [d  m][HD]  +  -Aq  [(ddp)  p]  —  A 163  [(ddp) 

-  AI 64  [( ddp)  p]  [r  ]  -  Al65  [(ddp)  p]  [HD] 

-  Al66  [(ddp)p]  [HT]  -  Al67  [(ddp)p]  [A  ( tt )] 

-  As  [(ddp)  p]  [a  (T  T)]  -  (Am  +  A70 )  [(ddp)  p]  [DT  ] 

-  Alll  [(ddp) p] [t2  (TT )]  —  \72  [(ddp)p]  [t2  (T i) 

-  As9  [( ddp)p ] 


(A. 328) 


(A. 329) 


(A. 330) 
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d[(dd//)d] 

dt 


d[(dd//)r] 

dt 


d\[dtfi)p\ 

dt 


/t53  [d//][d>2  (TT )J  +  /t54  [d//][d>2  (TT)J  -A0  [(dd//)  d] 

-  Aj73  [(dd//)  d][d/'  J  -  A[74  [(c/d//)  d][r*  ] 

-2^  [(dd//)d][/d2  (TT)] -\16  [(dd//)d][/Z2  (TT)] 

-\ni  [(dd//)d][d/D] -(2 118  +  Aj19)[(c/c///)c/][//r ] 

- \2Q  [(dd//)  d]  [DT]  -  421  [(dd//)  d]  [d2  (TT)] 

-^22[(dd//)d][r2(tl)]-A26o[(dd//)d] 


A57  [d//][l>r ]  -A0  [(dd//)r]-2 j83  [(dd//)r][d/‘] 

-  44  [(dd//). /]  [ZT  ]  -  4,5  [(dd//)  /]  [/d2  (TT)] 

-  4*6  [(dd//) r]  [h2  (TT)]  -  (^87  +  Agg )  [(dd//). r]  [dffi» ] 

-  49  [(dd//) r]  [a  (TT)]  -  ^90  [(dd//) r]  [a  (TT)] 

-/ll9l  [(dd//)r][d/T ]  —  \92  [(dd//)?][Z)r ]-/l26l  [(dd//)r] 


4  |/4]  ]  +  4  |/a]  [#°]  -  4  [idtM)p]  -  4s  [ ( d/ //)//]  [zr 


-  44  [(dr//)//]  [7”  ]  -  4s  [(*/*)  z]  [ZZD ]  -  Ai96  [(d?//)  //]  [AT ] 

-  47  [(dr//)  //]  [z»2  (TT)]— 4„  [(dr//)  /?]  [a  (T  T)] 


(^99  +  4o 


)  [(dtp)  p]  [DT ]  -  Am  [(dr//)  p ]  [r2  (TT ) 


-  4)2  [(dr//)  a][4  (TT)]  -  462  [(dr//)//] 


(A. 331) 


(A. 332) 


(A. 333) 
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d[(dtp)d~\i 

dt 


4  [d<u\ [DT]  +  A65  [tM][D2  (tt )]  +  A65  [tp] [D2  (TT)] 


-A0  [(dtp)d~^-A jgj  [(dtp) d~] ]//']-  A204  [(dtp)  j][7”] 
-a205  [(^//)flf][//2(TT)]-^06[(^//)j][//2(tl)] 


■  A207  [(dtp) d ]  [HD]  -  ( A208  +  A20g  ) [(dtp)d ]  [HT ] 

■  4o  [(dtp)  d]  [DT]  -  A2U  [(dfju)  d]  [t2  (TT)] 
-A212  [(dtp)d][T2  (TI)]-^63  [(dtp)d~\ 


4»M[-D7’]-^)  [(dtp)t\ 


-fi^J  =  A59[J//][r2(tt)]+260[jyu][r2(TI) 

-/I213  \{dtiu)t][H'^-A2U  [(dtp) t ] [z>* ]  -  A2l5  [(dtp)t][H2  (Tt) 

-  A216  [(dtp)t]  [h2  (TT)]  -  (47  +  A2I8 )  [(dtp)t\  [HD] 

-  A219  [(dtp)t] [  A  (TT)]  -  A220  [(dtp) t ] [A  (^)] 

-  A221  [(dtp)t\ [HT ]  - 1222  [(dtp) t ]  [DT ]  -  A264  [(dtp)t\ 


d[(ttp)p] 


dt 


=  4  M  [HT ]  ~  \  [(dp) p]  ~  43  [(dp)  p]  [D’  ] 

-A4  [(dp)p][TA~A 225  [(ttp)p][HD\ 

-  A226  [(dp)  p\ [HT  ]  -  A227  [(tip)  p]  [A  )] 

-  As  [(dp)p][D2  (TT)]  -  (A229  + A230  )[(ttp)  p]  [DT  ] 

-  4i  [(dp)p\ [t2  (TT )]  -  a232  [(ttp)  p]  [t2  (TT)] 
-A265  [(ttp)p] 


(A. 334) 


(A. 335) 


(A. 336) 
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^  -  =  ^70  M  [ DT ]  ~  4>  [(#/*)' d]  “  ^233  [(«/*)  <*][#’] 

-^34  [(«//)  ^][r  ]  -A235  [(«//)  </][#2  (TT )] 

-^\_{ttM)d\H2{^i)^-^1[(ttp)d'\[HD] 

-  (■ ^238  +  ^239  ) \{ttn  )< d  ]  [HT]  -  /L24o  [(* tt  p  )< d]  [DT] 

-  A241  [( ttju)d ]  \j2  (TT )]  -  A242  \{ttp)d~§T2  (T '!)] 
-^266  [(ttfi)d] 


^  ^  ^  =  K  M \t2  (TT)]  +  A72  [tp] [t2  (TT)]  - A0  [(ttp)t] 

-  ^43  [(ff/*  )  f  ]  [//*  ]  -  1244  \{ttp)t\  [D’  ] 

-A245  [(ttp)t]\_H2  (tt)]-^46  [(ttp) f][tf2  (TT) 

- (/L247  + 1248 ) \_{ttju)i] [HD]  -  A24g  [(«//)?] [d2  (TT ) 

-  A250  \{ttp)i]  [a  (TT)]  -  A25l  \{ttju)t]  [HT ] 

-A252  \{ttp)t~][DT]-Xm  [ \ttp)t ] 


d  [n]  ] 
dt 


^259 [{ddp)p]  +  A260Bn  ^(dd  p]d~]  +  A26lBn 


d[n2] 

dt 


A262  [(dtp)  p]  +  A263  [(dtp)d~]  +  A264  [{dtp)  t] 


<7[n3] 

dt 


2  {/L265  [{ttp)p]  +  A266  [{ttp) d]  +  A261  [{ttp)t^ 


(A. 337) 


(A. 338) 


(A. 339) 


(A. 340) 


(A. 341) 
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and 


-^a-iym+4  [//']+^  [Hi]+l>[HD]+X,lHT]+A,[iy~\ 

+  As [HD]  +  A,  [A ]H [DT]  + A,  [r  ]  + \t  [HT]  +  A,,  [DT]  +  4  [f2]} 

+ /L253  (i  -  )  [(^J//)  ^]  + /I254  (1  -  )  [(^J//)  J] 

+  ^255  (1  -®i)[(^//)?]  +  ^56(1-®2)(1-A)[(^)^] 

+  /L257  (1  -«2)(1-a)[(^)^]  +  ^258  (1  -®2)(1-a)[(^/0^]  (A. 3 42) 

+  {^59  (1  -  ®3  )  (1  -  Bn  )  +  ^59  (l  “  ®4  )  Bn  }  [(^//)  />] 

+  {^260  (l-®J)(l-^,)  +  4!o(1-®4)^}[(^)rf] 

+  {^6l(1-®3)(1-S„)  +  ;l26l(1-®4)5„}[(^)0 
+  ^262  (1  -  )  [(<*/*)  /?]  +  4,3  (!  -  )  [(dtfi)  d]  +  ^264  (*  ~  )  [(<*/*)  * ] 

+  A265  (l-co6)[(ttv)p]  +  A266  (l-co6)[(tt{j)d]  +  A267  (l-co6)[(tt{/)t] 


The  probability  density  of  free  muons  [//]  in  the  system  is  equal  to  the 
concentration  of  muons  input  in  the  chamber  minus  the  probability  the  muon  has  decayed 
or  is  bound  to  a  molecule  in  the  chamber.  This  is  written  as 


\.^=\_^^~[pA-[dA~[tA-\^ppti)p\-\ippti)d']-\ippti)t'\ 

-\[pdp)p]-\[pdp)d']-\{pdp)t\-\{ptp)p]-\[ptp)d~\ 

-[( ptju)t]  -  [( pt/Li)  p]-\{ptii)d\-  [( ptp)t\  - \{ddp)p]  (A. 343) 

-  [  3//e//]  -  [  ^Hep ] 


The  total  amount  of  muons  in  a  reaction  chamber  can  be  expressed  as 
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_^/]  =  [//o]e'V 


(A. 344) 


where  ju0  is  the  initial  concentration  of  negative  muons. 

Neglected  from  all  of  the  equations  in  this  appendix  are  the  effects  of  “muon 
scavenging.”  Muon  scavenging,  as  it  relates  to  muon-catalyzed  fusion,  refers  to  free 
muons  binding  to  elements  other  than  hydrogen,  and  muons  bound  to  hydrogen  being 
transferred  to  other  elements.  If  significant  quantities  of  helium  are  present  in  the 
reaction  chamber,  regardless  if  the  helium  is  a  product  of  fusion,  or  a  decay  product  of 
tritium,  then  muon  scavenging  needs  to  be  added  to  the  rate  equations.  Muon  scavenging 
also  needs  to  be  considered  if  there  is  a  significant  probability  of  muonic  molecular 
molecules  interacting  with  chamber  walls. 

Most  of  the  reaction  rate  constants  presented  in  this  appendix  have  been  looked  at 
computationally,  experimentally  or  both,  however,  the  accuracy  to  which  the  values  are 
known  varies  greatly.  Some  possible  reactions  which  are  not  believed  to  occur  in  any 
significant  amount  ( e.g .,  ddfi—>T'+pfi)  have  been  omitted.  Some  of  the  reactions 
listed  (e.g.,  t/u-^d/u)  occur  so  seldom  that  they  can  be  neglected.  The  reaction  rate 
constants  that  have  the  greatest  need  for  further  study  are  those  that  form  tri-atomic 
muonic  molecules.  It  has  been  shown  that  the  non-fusing  hydrogen  in  tri-atomic  muonic 
molecules  has  a  significant  impact  on  the  fonnation  rate  of  these  molecules. [6;  7]  The 
mass  of  the  non-fusing  atom  affects  the  separation  distance,  vibrations,  and  stability  of 
the  fusing  atoms,  and  as  a  result  affects  the  fusion  rate. 
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Appendix  B.  Equilibrium  Concentration  of  Hydrogen  Isotopes 


In  order  to  determine  the  relative  amounts  of  various  hydrogen  molecules  which 
will  be  present  in  a  reaction  chamber  at  a  given  temperature  and  pressure,  it  is  necessary 
to  know  the  equilibrium  constants  (Keq)  of  the  isotopes  involved.  Equilibrium  constants 
are  most  often  calculated  from  the  Gibb’s  free  energy  (AG)  by  the  equation:[l;  2] 

AG  =  —RT  In  Keq  (B.l) 

where  R  =  8.3 144  J/mole  K  and  T  is  the  temperature  in  Kelvin. 

Gibb’s  free  energy  can  be  detennined  from  partition  functions. [3]  For  non-tritium 
containing  hydrogen  molecules  these  values  have  been  published  in  the  Journal  of 
Physical  and  Chemical  Reference  Data. [4] 

Protium  (P),  deuterium  (D),  and  tritium  (7)  combine  to  form  a  variety  of 
hydrogen  molecules.  In  order  to  determine  the  equilibrium  concentration  of  these 
molecules  it  is  necessary  to  consider  the  reactions: 

HD - >\h2  +  4d2  (B.2) 

HT - >\h2+\t2  (B.3) 
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DT 


2 


2 


A 


>D 


and 


2 


Using  reaction  equations  B.l  to  B.6  it  is  possible  to  derive  the 
equilibrium  equations: 


K 


C,H  ~ 


K 


C,D  ~ 


K 


C,T  ~ 


(B.4) 


(B.5) 


(B.6) 


(B.7) 


following 


(B.8) 


(B.9) 


(B.10) 
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(B.ll) 


K 


C , HD 


c 

^  urt 


K, 


C,HT 


c 


(B-12) 


HT 


and 


K. 


C,DT 


r 

nr 


(B.  13) 


These  equations  do  not  include  all  possible  molecules  that  form;  however,  they  do 
include  all  of  the  molecules  that  occur  in  quantities  that  can  affect  nonnal  muon- 
catalyzed  fusion  reaction  conditions.[l;  5:20] 

From  the  initial  concentrations  of  H2,  D2,  and  T2  the  following  can  be  derived. 


C  = C  --C  --C  --C 

H2  i,H  2  H  2  HD  2  HT 

C  =C  --C  --C  --C 

^D2  ^ i,D  2  D  2  HD  2  DT 


(B!4) 


(B!5) 
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(B.  16) 


cT:=c,T-jcr-jCHr-jcDT 

Where  Ch,  Cd,  Ct,  Chd,  Cht,  and  Cdt  represent  equilibrium  concentrations  and  Q#,  Q/j, 
and  Cjj  represent  the  initial  molecular  concentrations  of  Hj,  D2,  and  77  respectively. 

Kc  is  the  molar  (n)  equilibrium  constant  and  can  be  calculated  for  the  equilibrium 
pressure  constant  Kp.[ 6: 169-179]  Assuming  ideal  gas  behavior,  Kp  =  Keq. 

Kc  =  P~MKeq  (B-17) 

At  temperatures  less  than  2000  K  only  diatomic  hydrogen  molecules  occur  in 
significant  quantities  and  many  of  the  tenns  in  the  above  equations  approach  zero. 
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